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DIRECTOR’S  OVERVIEW 


This  report  presents  a  summary  of  the  scientific  progress  and  accomplish¬ 
ments  on  research  projects  funded  by  the  Joint  Services  Electronics  Program 
(JSEP)  for  the  contract  period  from  l  January  1987  through  31  December 
1987.  It  does  not  contain  information  regarding  accomplishments  on  research 
projects  funded  in  other  ways. 

The  Joint  Services  Electronics  Program  at  Polytechnic  University  is  the 
core  of  interdisciplinary  research  In  electronics  encompassing  programs  in  the 
Departments  of  Electrical  Engineering,  Physics  and  Chemistry  under  the  aegis 
of  the  Weber  Research  Institute.  The  research  encompassed  by  this  program  is 
grouped  into  three  broad  categories:  Electromagnetics,  Solid  State  Electronics 
and  Information  Electronics.  The  detailed  projects  (research  units)  comprising 
the  complete  program  are  listed  In  the  Table  of  Contents. 

Following  our  previous  practice.  Section  2  of  every  research  unit  contains  a 
short  summary  of  the  recent  progress.  Further  details  regarding  that  progress 
are  contained  in  Section  3,  State  of  the  Art  and  Progress  Details. 

Each  of  the  research  units  described  in  this  report  is  designated,  for  exam¬ 
ple,  as  EM7-1,  SS7-1  or  IE7-1,  corresponding  to  the  category,  year  (1987),  and 
number  within  the  category.  These  numbers  follow  the  numbers  given  in  our 
annual  proposal  dated  December  31,  1988,  with  two  exceptions.  In  accord  with 
instructions  of  the  Technical  Coordinating  Committee  we  phased  out  the  Infor¬ 
mation  Electronics  unit  referenced  in  the  proposal  as  IE7-1  before  the  start  of 
the  period  covered  by  this  report.  Consequently,  we  have  renumbered  the 
other  two  Information  Electronics  units.  The  numbers  were  reduced  by  one: 
the  proposed  IE7-2  Is  reported  here  as  IE7-1;  and  the  proposed  IE7-3  is 
reported  here  as  IE7-2. 

Following  this  Overview,  we  present  the  Report  on  Significant  Scientific 
Accomplishments,  which  highlights  a  few  of  our  more  important  contributions 
this  year. 


REPORT  ON  SIGNIFICANT  ACCOMPLISHMENTS 
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(December  31,  1987)  ! 
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Focus  Wave  Mode  Causality  Dilemma  1 

I 

Professors  E.  Heyman  and  L.B.  Felsen 
Unit  EM 7-3 

Conductor-Backed  Lines:  Dangers  and  Full  Wave  Analysis 

Professor  AjV.  Oliner 
Unit  EM 7-1 

Optical  Multistability  of  Anistropic  Microparticles 

Professor  K.M.  Leung 
Unit  SS7-2 


1.  FOCUS  WAVE  MODE  CAUSALITY  DILEMMA 


Professors  E.  Heyman  and  L.  B.  Felsen 

Efficient  transfer  of  electromagnetic  energy  in  space  has  motivated  the 
study  of  time-dependent  wave  fields  that  propagate  in  missile-like  fashion 
without  diffusion  or  dispersion.  Among  new  wave  functions  addressing  these 
objectives,  the  focus  wave  modes  (FWM)  have  received  major  attention  [1,2]. 
These  wave  objects  are  exact  solutions  of  the  time-dependent  source-free  field 
equations,  and  they  have  infinite  energy.  They  travel  through  space  at  the  speed 
of  light  with  undistorted  profiles  that  match  the  focal  region  of  a  well- 
collimated  beam.  These  features  are  highly  desirable  for  directed  energy 
transfer. 

A  major  question  has  been  whether  and  how  the  FWM  can  be  excited  by 
causal  source  distributions  in  a  physical  aperture.  Through  several  alternative 
studies  involving  Green’s  function  and  spectral  techniques,  we  have  been  able 
to  show  that  backward  propagating  wave  phenomena  are  dominant  in  the  syn¬ 
thesis  of  a  forward  propagating  FWM  [3,4].  While  this  behavior  is  compatible 
with  the  source-free  regime,  it  can  not  be  established  by  causal  initial  condi¬ 
tions  in  an  aperture  plane.  This  Important  conclusion  has  led  to  a  reassessment 
of  the  role,  if  any,  to  be  played  by  the  FWM  in  physical  pulsed  beam  propaga¬ 
tion. 

REFERENCES 

1.  J.  N.  Brittingham,  "Focus  Wave  Modes  in  Homogeneous  Maxwell’s  Equa¬ 
tions:  Transverse  Electric  Mode,"  J.  Appl.  Phys.,  Vol.  54,  1983,  pp.  1179- 
1189. 

2  R.  W.  Ziolkowski,  "Exact  Solutions  of  the  Wave  Equation  with  Complex 
Source  Locations,"  J.  Math.  Phys.,  Vol.  26,  1985,  pp  881-883. 

3.  E.  Heyman,  B.  Z.  Steinberg  and  L.  B.  Felsen,  "Spectral  Analysis  of  Focus 
Wave  Modes,"  J.  Opt.  Soc.  Am.  A,  Vol.  4,  1987,  pp.  2081-2091. 

4  E.  Heyman,  "Focus  Wave  Modes:  A  Dilemma  with  Causality,"  submitted 
to  IEEE  Trans.  Antenna  &  Propagat 


2.  (X)NDUCTOR- BACKED  LINES:  DANGERS  AND  FULL  WAVE 
ANALYSIS 

Professor  Arthur  A.  Oliner 

In  millimeter-wave  integrated  circuits  and  in  monolithic  microwave  circuits 
(MMICs),  it  is  a  common  complaint  that  components  that  work  well  separately 
fail  when  combined  together.  The  most  common  cause  of  such  failure  is  that 
unexpected  and  unwanted  radiation  is  produced  that  results  in  cross  talk  and 
other  deleterious  effects.  These  ill  effects  occur  in  pronounced  fashion  when 
certain  transmission  lines,  such  as  slot  lines  and  coplanar  waveguides,  are 
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conductor  backed. 


It  is  very  tempting  to  introduce  such  conductor  backing  because  such  back¬ 
ing  has  many  advantages,  such  as  a  lower  Z0,  reduced  dispersion,  better 
mechanical  strength,  and  so  on.  Raytheon  is  planning  to  employ  conductor- 
backed  coplanar  waveguide  in  its  MMICs  because  of  easier  fabrication  and 
lower  cost,  in  addition  to  these  advantages.  However,  there  are  also  certain 
important  dangers,  or  serious  potential  problems,  that  are  introduced  by  the  use 
of  conductor  backing.  They  include  leakage  of  power  into  surface  waves  or  into 
the  dielectric  region  between  the  plates,  unexpected  cross  talk,  significant 
alteration  of  the  guide  wavelength,  and  unexpected  or  unwanted  coupling  to 
neighboring  lines. 

These  possible  unpleasant  consequences  are  generally  not  recognized  or 
discussed  in  the  literature,  except  for  a  single  recent  paper  which  provided  no 
quantitative  information.  The  present  study  has  developed  two  different  full- 
wave  analyses  for  the  class  of  conductor-backed  lines,  and  has  applied  them  to 
slot  lines  as  an  important  example.  A  set  of  measurements  was  also  taken  on  a 
particular  conductor-backed  slot  line  structure.  The  results  obtained  show  that 
the  numbers  obtained  from  the  two  different  theoretical  approaches  agree  well 
with  each  other  and  with  the  measurements  taken. 

Several  very  interesting  effects  were  found,  including  the  fact  that  the 
power  leakage  rate  is  rather  high.  The  results  verify  the  concern  expressed 
above  regarding  the  potential  dangers  in  the  use  of  conductor-backed  lines,  and 
the  study  itself  provides  the  first  quantitative  information  regarding  the  serious¬ 
ness  of  these  effects. 


3.  OPTICAL  MULTISTABILITY  OF  ANISTROPIC  MICROPARTICLES 

Professor  K.M.  Leung 


It  was  shown  in  our  previous  work  that  a  Rayleigh-sized  spherical  micro¬ 
particle  can  exhibit  optical  bistable  behavior  and  can  also  function  as  an  optical 
transistor.  In  order  to  enhance  the  functionality  and  capabilities  of  microparti¬ 
cle  photonics,  it  is  Important  to  investigate  more  complex  microparticle  com¬ 
ponents.  For  this  purpose,  we  have  recently  extended  our  work  from  spherical 
to  spheroidal  microparticles. 

In  comparison  with  a  spherical  particle,  a  spheroidal  particle  has  a  few  extra 
physical  characteristics,  such  as  the  particle  shape  and  the  particle  orientation 
with  respect  to  the  polarization  of  the  incident  field.  As  a  result,  it  is  capable  of 
exhibiting  a  richer  class  of  intrinsically  nonlinear  phenomena.  We  find  that 
besides  the  usual  bistable  behavior,  a  total  of  five  different  types  of  trtstable 
behaviors  are  possible,  depending  on  the  frequency  of  the  incident  light  and  the 
particle  shape  and  orientation. 
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SECTION  I:  ELECTROMAGNETICS 


A.  NEW  PHYSICAL  EFFECTS  INVOLVING  OPEN 
DIELECTRIC  STRUCTURES 


Professors  A.A.  Oliner  and  M.  Guglielmi 
Unit  EM7-1 


1.  OBJECTTVE(S) 


This  study  Is  concerned  with  guiding,  radiating  and  scattering  effects 
involving  open  dielectric  structures.  The  majority  of  electromagnetic 
waveguides  are  open  and  involve  dielectric  materials,  and  they  therefore  fall 
into  this  category.  Certainly,  most  waveguides  for  millimeter-wave  integrated 
circuits  and  all  waveguides  for  integrated  optics  are  in  this  category.  Moreover, 
problems  such  as  cross  talk,  leakage,  and  low  Q,  that  arise  In  microwave 
integrated  circuits,  ‘ and  coupling  effects  and  blind  spots  that  appear  in  many 
phased-array  antennas,  are  due  to  the  presence  of  dielectric  material  on  those 
open  structures.  It  should  be  clear,  therefore,  that  these  basic  studies  have 
important  practical  implications. 

Before  our  investigation,  almost  all  of  the  published  literature  on  these 
topics  involved  a  surface  wave  Incident  on  a  dielectric  Junction  or  grating  at 
normal  incidence,  where  the  resulting  boundary- value  problem  is  two- 
dimensional,  and  TE  and  TM  modes  do  not  couple  to  each  other.  In  many 
open  dielectric  waveguides  and  antennas  for  millimeter-wave  integrated  circuits 
and  for  integrated  optics,  however,  the  surface  wave  is  obliquely  incident.  The 
boundary- value  problems  then  change  from  the  scalar,  two-dimensional  ones  to 
vector  three-dimensional  ones,  and  the  TE  and  TM  modes  no  longer  remain 
independent  but  are  coupled  together.  The  vector  nature  of  these  problems 
Introduces  new  mathematical  challenges,  but,  more  interestingly,  the  mode 
coupling  produces  a  rich  variety  of  interesting  and  somettmes  unexpected  new 
physical  effects.  In  our  studies  so  far,  we  have  found  that  a  number  of  such 
effects  arise  in  potentially  important  guiding  and  radiating  structures,  either 
leading  to  possible  performance  problems  or  providing  new  opportunities  if 
those  effects  can  be  properly  utilized. 

We  are  continuing  to  explore  these  physical  effects  on  old  and  new  types 
of  dielectric  structure.  It  is,  of  course,  necessary  to  first  obtain  the  appropriate 
mathematical  solutions;  those  solutions,  or  the  approaches  used  to  obtain  them, 
are  often  of  interest  in  themselves.  We  then  examine  the  physical 
consequences  of  those  solutions,  and,  where  appropriate,  assess  their 
implications  for  device  performance  in  millimeter-wave  integrated  circuits  and 
antennas  and  in  integrated  optics. 
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SECTION  I:  ELECTROMAGNETICS 
2.  SUMMARY  OF  RECENT  PROGRESS 

Major  effort  during  the  past  year  was  placed  on  six  topics,  three  of  which 
are  related  to  aspects  of  an  important  basic  antenna  type  for  millimeter  waves. 
The  three  related  aspects  of  the  antenna  investigation  are  described  first,  and 
they  are  grouped  together  under  top'c  A  in  section  3.  The  antenna  is  the 
dielectric  image  guide  leaky-wave  antenna  loaded  by  periodic  metal  strips,  for 
which  many  measurements  have  been  presented  in  the  literature  but  for  which 
no  theory  was  previously  available  from  which  the  antenna  could  be  designed. 
The  three  related  aspects  under  study,  and  discussed  under  section  3,  are  first, 
a  detailed  parametric  analysis  of  the  antenna  behavior,  second,  a  novel 
optimization  procedure  for  tapering  the  antenna  to  control  the  sidelobes,  and 
third,  a  new  solution  for  scattering  by  a  multimode  grating  that  yields  arbitrary 
accuracy.  The  other  three  topics  are:  a  new  type  of  two-dimensional  scanning 
array  for  millimeter  waves,  a  leaky-wave  analysis  of  a  high-gain  printed  antenna 
configuration  that  employs  a  superstrate,  and  a  study  of  conductor-backed  slot 
line  for  MMIC  and  microwave  integrated  circuits  that  provides  full-wave 
analyses  for  the  first  time. 

For  the  first  topic,  that  of  a  detailed  parametric  analysis  of  the  behavior  of 
the  dielectric  image  guide  leaky-wave  antenna,  the  phase  and  leakage  constants 
of  the  leaky  structure  were  determined  as  a  function  of  frequency,  dielectric 
constant,  and  geometric  dimensions.  From  the  trends  in  behavior,  we  are  now 
able  to  design  antennas  to  meet  specified  performance  characteristics, 
something  certainly  not  available  previously.  As  an  example,  we  present  in 
section  3  the  curve  obtained  when  we  wish  to  optimize  the  bandwidth  over 
which  the  beam  width  remains  constant. 

The  second  topic  involves  a  new  method  for  optimizing  the  taper  in  the 
aperture  distribution  of  the  leaky-wave  antenna  studied  above  in  order  to 
control  the  sidelobe  level  in  a  specified  way.  As  a  challenging  example,  we 
chose  a  Taylor  distribution  corusponding  to  sidelobes  which  are  all  35  dB 
below  the  main  beam  maximum.  The  problem  Is  to  specify  the  taper  in  the 
dimensions  that  satisfied  the  Taylor  distribution  requirement,  but  the  difficulty 
is  that  only  the  width  and  location  of  the  metal  strips  are  accessible  parameters. 
The  novel  design  procedure  we  developed  changes  both  the  strip  widths  and 
the  strip  period  simultaneously  to  accomplish  this  purpose.  Since  the  strips  can 
be  deposited  lithographically,  or  etched  away,  the  new  procedure  allows  one  to 
design  a  mask  that  accomplishes  the  geometrical  arrangement  all  at  one  time. 
The  design  procedure  and  the  fabrication  process  are  therefore  practical  for  use 
at  millimeter  wavelengths. 

The  third  topic  is  only  ind.rectly  related  to  the  antenna  itself,  since  it 
involves  Just  a  constituent  of  the  antenna,  namely,  the  metal  grating  placed  on 
an  air-dielectric  interface.  As  a  part  of  the  analysis  of  the  antenna,  we  had 
derived  a  solution  for  plane-wave  scattering  by  such  a  grating.  We  are  currently 
in  the  process  of  deriving  an  alternative  solution,  which  involves  a  different 
series  expansion  technique  and  has  the  promise  of  yielding  arbitrary  accuracy. 
The  earlier  solution  utilized  a  small-argument  expansion. 
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SECTION  I:  ELECTROMAGNETICS 


The  next  topic  actually  makes  use  of  the  dual  to  the  solution  referred  to 
above,  that  is,  the  solution  previously  obtained  for  plane-wave  scattering  by  the 
metal-strip  grating  on  an  air-dielectric  Interface,  but  for  the  opposite  incident 
polarization.  This  solution  forms  part  of  the  overall  solution  for  a  new  type  of 
two-dimensional  scanning  antenna,  with  leaky-wave  scan  in  elevation  and 
phase-shift  scan  in  the  cross  plane.  The  antenna  possesses  several  key 
advantages:  no  grating  lobes,  no  blind  spots,  and  no  cross  polarization.  The  use 
of  a  periodic  structure  (grating)  permits  a  wider  scan  range,  as  discussed  in 
section  3. 

The  fifth  topic  relates  to  a  superstrate-substrate  structure  that  permits  a 
simple  printed-circuit  antenna  to  radiate  a  narrow  beam  that  scans  with 
frequency.  The  new  feature  here  is  that  we  show  that  this  performance  can  be 
elegantly  described  in  terms  of  leaky  waves.  It  is  first  shown  that  the  leaky- 
wave  approach  reproduces  all  the  asymptotic  expressions  and  numerical  values 
for  the  far  field  obtained  by  means  of  a  straightforward  field  approach,  but 
derived  in  a  simpler  way,  using  leaky  waves.  More  important,  however,  is  that 
the  leaky-wave  approach  furnishes  insight  into  the  physical  processes  involved, 
and  also  provides  new  practical  information,  as  explained  in  section  3. 

The  sixth  and  last  topic  applies  to  a  set  of  potentially  serious  problems  that 
arise  in  MMICs  and  microwave  integrated  circuits.  Conductor  backing  has 
been  Introduced  into  such  circuits  when  transmission  lines  such  as  slot  line  and 
co planar  waveguide  are  used,  because  of  the  many  advantages  that  are 
provided.  There  are,  on  the  other  hand,  possible  unpleasant  and  unexpected 
consequences  to  their  use  that  seem  not  to  be  fully  recognized  or  understood. 
In  addition,  there  are  no  full-wave  analyses  available  that  permit  us  to  analyze 
quantitatively  how  serious  these  effects  may  be.  Some  of  these  potentially 
serious  effects  are  described  In  section  3.  We  have  developed  two  full-wave 
theoretical  approaches,  one  purely  numerical  and  the  other  in  network  form 
and  analytical  in  nature,  that  agree  well  with  each  other  and  with  measurements 
in  a  special  case.  These  theoretical  approaches  are  currently  being  employed  to 
determine  quantitatively  the  physical  behavior  corresponding  to  some  of  these 
potential  surprises  in  specific  cases  in  which  either  the  slot  line  conductors  or 
the  backing  conductor  are  of  Infinite  or  finite  lateral  extent. 


3.  STATE  OF  THE  ART  AND  PROGRESS  DETAILS 

A.  Dielectric  Image  Guide  Leaky-Wave  Antennas 
Loaded  by  Periodic  Metal  Strips. 

In  last  year's  Annual  Report  we  presented  a  simple  and  accurate  transverse 
equivalent  network  for  the  class  of  leaky-wave  antennas  shown  in  Fig.  1. 
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Fig.  1.  Leaky-wave  antenna  consisting  of  a  periodic  array  of  metal  strips 

placed  on  a  dielectric  image  guide. 

The  theory  developed  allowed  us  to  understand  its  basic  behavior  and  to  exert 
a  great  degree  of  control  over  its  performance,  showing  for  the  first  time  how 
this  class  of  antennas  can  be  effectively  designed  to  achieve  desired 
characteristics. 

The  details  of  the  theory  developed  and  some  of  the  results  obtained  were 
presented  at  a  conference  [l],  at  which  the  paper  was  awarded  the  Second  Prize 
in  the  1987  U.R.S.I.  National  Student  Paper  Prize  Competition.  Additional 
results  and  some  applications  of  the  theory  developed  were  also  presented  at  a 
later  conference  [2]. 

In  the  past  year,  our  research  efforts  in  connection  with  this  class  of 
antennas  have  been  aimed  in  three  main  directions.  First,  we  have  performed 
a  detailed  parametric  analysis  of  the  antenna  behavior  when  the  grating 
discontinuity  is  uniform  along  the  z  direction  of  Fig.  1.  Second,  we  have 
studied  suitable  non-uniform  (tapered)  grating  geometries  to  achieve  control 
over  the  sidelobe  levels.  Finally,  we  have  begun  the  development  of  a  new 
solution  for  the  description  of  the  grating  discontinuity  that  will  allow  us  to 
compute  to  an  arbitrary  degree  of  accuracy  the  elements  involved  in  its  network 
description. 

(1)  Description  of  the  Parametric  Analysis 

The  first  step  in  the  parametric  analysis  of  this  class  of  antennas  is  the 
derivation  of  a  suitable  set  of  constraints  that  can  facilitate  the  initial  choice  of 
the  structural  parameters.  Keeping  in  mind  that  we  want  to  have  only  one 
main  radiation  lobe  In  the  air  region  in  Fig.  1,  we  have  derived  the  approximate 
expressions  shown  in  Table  1. 
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Single 
Radiating 
Space  Harmonic 

Single  Mode  in  the 
Dielectric  Layer 


Backward 

Quadrant 


er>l 


P  “  2y/tx-  1 


Forward 

Quadrant 


er>q 


t  ^  yfh-  1 


Table  1. 


The  expressions  shown  can  be  used  to  choose  a  set  of  parameters 
so  that  the  antenna  in  Fig.  1  has  only  one  main  radiation  lobe. 


Next,  choosing  a  set  of  parameters  appropriate  for  operation  in  the 
backward  quadrant,  we  have  obtained  the  antenna  dispersion  curve  shown  in 
Fig.  2.  We  then  perturbed  the  structural  parameters  one  at  a  time,  observing 
that  the  real  part  of  the  complex  propagation  constant  remains  essentially 
unchanged  while  the  imaginary  part,  the  leakage  rate  a,  can  be  strongly 
modified.  The  results  obtained  can  be  summarized  as  follows: 

a.  Increasing  the  obstacle  size  d/p  generally  Increases  the  leakage  rate  but  its 
frequency  dependence  remains  essentially  the  same. 

b.  Decreasing  the  thickness  d  of  the  dielectric  layer  produces  a  twofold  effect: 
the  curve  for  a,  the  leakage  rate,  "flattens"  and  its  average  value  increases. 

c.  Decreasing  the  value  of  the  dielectric  constant  er  results  in  a  generally 
"milder"  leakage  behavior. 

As  a  result  of  the  accurate  parametric  analysis  performed,  we  are  now  able 
to  design  antennas  that  meet  specified  radiation  characteristics.  As  an  example, 
we  show  in  Fig.  3  the  results  obtained  when  we  optimize  the  antenna  behavior 
for  a  constant  half-power  beamwidth  over  a  maximized  scan  range. 

(2)  Optimized,  Aperture  Distribution  for  Sidelobe  Control 

The  antennas  discussed  in  the  previous  section  are  uniform  in  the 
longitudinal  direction  (z  direction  in  Fig.  I).  As  a  result  of  this  uniformity  the 
antenna  aperture  distribution  is  exponentially  decaying,  and  therefore  not 
appropriate  for  sidelobe  control.  To  achieve  a  specified  sidelobe  characteristic 
one  requires  a  tapered  aperture  distribution.  This  means  that  the  obstacle  size 
d/p  must  be  varied  along  the  z  direction  in  Fig.  1. 

The  theory  that  we  have  developed  for  this  class  of  antennas  allows  us  to 
design  such  tapered  geometries  using  the  assumption  that  the  structure  can  be 
considered  locally  uniform.  Th.s  assumption,  however,  can  be  introduced  only 
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If  the  rate  of  change  of  d/p  along  the  z  direction  is  small.  In  addition, 
changing  the  obstacle  size  d/p  not  only  affects  the  leakage  rate  a  but  also  the 
phase  constant  (30  of  the  structure.  As  a  result,  different  sections  of  the 
antenna  would  have  different  phase  constants  and  would  radiate  in  slightly 
different  directions,  thereby  degrading  the  overall  antenna  performance. 

It  is  therefore  evident  that  changing  d/p  alone  is  not  sufficient  for  this  type 
of  optimization.  A  much  better  result  is  obtained  if  one  varies  at  the  same 
time  both  the  obstacle  size  and  the  distance  between  adjacent  metal  strips, 
thereby  obtaining  the  desired  antenna  in  terms  of  a  locally  uniform  and  periodic 
structure. 

As  an  example  of  this  procedure,  we  have  synthesized  the  required 
antenna  geometry  corresponding  to  the  Taylor  aperture  distribution  appropriate 
to  a  35  dB  sidelobe  level.  Figure  4  presents  the  radiation  pattern  thus  obtained. 

One  feature  that  deserves  more  attention  is  the  actual  practical 
implementation  of  the  required  taper.  Once  a  given  aperture  distribution  is 
specified,  the  theory  developed  yields  the  width  of  each  metal  strip  and  the 
distance  between  adjacent  strips  for  the  whole  structure.  A  mask  can  then  be 
realized  so  that  the  whole  system  of  radiating  metal  strips  can  be  readily 
deposited  on  the  top  surface  of  the  basic  guiding  structure.  This  feature  is 
particularly  attractive  at  millimeter  wavelengths  where,  due  to  the  reduced 
dimensions,  mechanical  simplicity  is  an  important  issue. 

The  work  in  connection  with  this  optimization  procedure  has  been  done  in 
collaboration  with  Professor  J.  Encinar,  who  visited  Polytechnic  University  as  a 
NATO  Postdoctoral  Fellow  from  Spain  during  1087. 

(3)  Scattering  by  a  Multimode  Grating  at  a  Dielectric  Interface: 

A  New  Solution  that  Yields  Arbitrary  Accuracy 

In  a  recent  Ph.D.  thesis  [3],  we  developed  a  network  description  for  a 
II2I!  multimode  metal-strip  grating  at  a  dielectric  Interface  that  is  valid  in  the 
small-argument  limit.  This  solution  was  obtained  by  solving  a  Cauchy  singular 
integral  equation  of  the  type  shown  in  (1). 


.  ,  z»\  1  l  Gn(a^)  J(., 

r”(“«  -  —  /  -^r-ae 


where  fn(a£)  is  known  and  Gn(a£')  is  the  unknown.  To  obtain  an  explicit 
expression  for  Gn(a£')  we  expanded  the  known  function  fa(a£)  in  a  Taylor 
scries  with  respect  to  a  relevant  "small"  parameter  a,  obtaining 

fn(aO=  a^n)  -i-a1(n)a^  +  a,(n)(a02  +•*'  (2) 

Equation  (2)  was  then  used  to  obtain  the  analytical  solution  of  (1).  This 
procedure  led  in  a  few  simple  steps  to  an  approximate  closed  form  expression 
for  the  unknown  Gn(o:£)  but  had  a  limitation.  When  the  number  of  terms 
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Included  In  (2)  was  increased,  we  ran  into  a  serious  problem  in  evaluating  the 
proper  expansion  coefficients. 

In  order  to  solve  this  problem  we  have  now  adopted  a  different  series 
expansion  technique.  Using  Legendre  Polynomials,  we  can  write 

fn(a£)  =  S  bk(n)Pk(ae)  (3) 

k  =0 

where  the  coefficients  bk(n)  are  formally  given  by  [4] 

i 

bk(n)  =  (k+“)  /  fn(a0  Pk(af)d  af  (4) 

z  -l 

The  integral  in  [4]  can  now  be  obtained  in  closed  form  for  arbitrary  values 
of  k,  thus  permitting  the  evaluation  of  the  function  Gn(a£)  to  an  arbitrary 
degree  of  accuracy. 

At  the  present  time  we  have  completed  all  of  the  analytical  work  in 
connection  with  this  new  procedure,  and  we  are  in  the  process  of  developing  a 
code  to  verify  the  accuracy  of  the  results  obtained.  This  work  has  been  done  in 
collaboration  with  Professor  H.  Hochstadt  of  the  Mathematics  Department  of 
Polytechnic  University. 

B.  A  Novel  Scanning  Array  for  Millimeter  Wavelengths 

The  scanning  array  that  was  devised  and  analyzed  during  the  past  year  is  a 
modification  of  the  class  of  scanning  arrays  studied  by  us  earlier  [5,8].  Those 
arrays  consisted  of  a  linear  array  of  leaky-wave  line  sources,  where  scanning  in 
elevation  is  achieved  by  frequency  scan  or  electronic  scan  of  the  leaky-wave 
line  sources,  and  scanning  in  the  cross  plane  is  obtained  by  placing  phase 
shifters  at  the  feed  ends  of  the  successive  line  sources  in  the  array.  The  angle 
in  the  cross  plane  then  depends  on  the  phase  shift  introduced  between  the 
successive  line  sources.  The  radiation  in  pencil  beam  form  will  scan  in  both 
elevation  and  azimuth  in  a  conical  scan  manner.  The  leaky-wave  line  sources 
were  made  to  be  uniform  along  the  length  of  the  line  sources  in  order  to 
provide  simplicity  in  their  configuration.  As  a  result,  however,  the  scan  in 
elevation  is  restricted  to  the  forward  quadrant  only;  the  two-dimensional  scan 
provided  by  the  whole  array  then  covers  only  half  of  space  at  best. 

The  leaky-wave  line  source  in  the  array  discussed  here  is  periodic  instead  of 
uniform  longitudinally,  consisting  of  a  periodic  array  of  metal  grids  or  slots,  as 
shown  in  Fig.  5.  The  periodic  grid  is  in  fact  still  simple  in  configuration,  and  it 
can  be  fabricated  at  one  time  by  lithographic  means.  The  mask  for  such 
fabrication  can  incorporate  the  tapering  of  the  aperture  distribution  for  sideiobe 
control,  following  our  improved  procedure  in  Sec.  A,  (2). 

The  periodic  nature  of  the  leaky-wave  line  source  permits  greater  flexibility 


SECTION  I:  ELECTROMAGNETICS 


in  scan  coverage.  If  the  period  of  the  periodic  strips  were  sufficiently  large 
(somewhat  greater  than  a  half  guide  wavelength  along  the  line  source  axis),  the 
n  =  -1  space  harmonic  of  the  periodic  array  would  become  fast  and  would  leak 
power  away  at  some  angle.  The  angular  coverage  in  elevation  would  then 
comprise  the  whole  background  quadrant  and  some  or  all  of  the  forward 
quadrant,  depending  on  the  parameter  values. 

It  is  easy  to  dream  up  a  variety  of  array  structures  that  should  work  well, 
but  we  are  also  restricted  by  the  fact  that  we  must  be  able  to  analyze  the  array 
accurately.  In  the  antenna  in  Fig.  5,  the  solution  for  the  array  of  metal  strips 
placed  on  a  dielectric  layer  comes  from  the  same  general  integral  equation 
solution  as  the  one  we  used  for  the  antenna  described  under  A  above,  with  the 
important  difference  that  the  polarization  is  opposite.  The  form  of  the 
equivalent  network  is  similar,  however,  as  are  many  of  the  basic  features. 

In  the  array  In  Fig.  5,  the  individual  line  sources  each  are  fed  by  a 
dielectric-filled  rectangular  waveguide  placed  on  its  side,  so  as  to  create  a 
horizontal  electric  field  polarization.  The  electric  field  lines  are  then 
perpendicular  to  the  metal  baffles,  so  that  the  spacing  between  them  can  be 
whatever  we  wish;  in  particular,  we  can  make  them  less  than  about  X0/2  so  that 
grativj  lobes  are  automatically  avoided  .  We  can  also  view  the  structure  as  a 
grounded  dielectric  layer  with  a  metal  grating  on  it.  The  periodic  structure 
produces  an  infinity  of  space  harmonics,  all  of  which  will  be  slow  if  the  grating 
period  is  too  small.  For  each  of  these  space  harmonics,  including  the  n  =  0 
one,  the  field  will  then  be  above  cutoff  in  the  vertical  direction  in  the  dielectric 
regions  but  evanescent  in  the  air  regions.  If  the  period  is  increased  sufficiently, 
or  alternatively  If  the  frequency  is  increased  sufficiently,  one  of  these  space 
harmonics  (the  n  =  -l  one)  will  become  fast  and  will  become  propagating  in  the 
air  regions  in  the  vertical  direction.  As  a  result,  a  mode  akin  to  a  TEM  mode 
in  each  parallel  plate  region  will  propagate  at  an  angle  toward  the  open  top  and 
will  leak  power  away. 


a 

Fig.  5.  Configuration  of  the  novel  two-dimensional  scanning  array 

composed  of  a  linear  array  of  leaky-wave  line  sources,  for  which 
there  are  no  grating  lobes,  no  blind  spots,  and  no  cross 
polarization. 
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Tlie  analysis  of  this  array  and  numerical  results  for  it  show  that  when  the 
ratio  of  strip  width  to  period  is  varied  the  value  of  a  changes  greatly  but  the 
value  of  0  is  modified  only  very  slightly.  This  behavior  is  Just  what  we  seek  so 
that  the  beam  width  and  the  angle  of  the  beam  maximum  can  be  specified 
essentially  independently,  and  so  that  the  antenna  aperture  can  be  tapered  to 
control  the  sidelobe  level. 

For  the  dielectric  constant  value  selected  (er  =  4.0),  this  array  can  scan  in 
elevation  essentially  only  in  the  backward  quadrant,  with  only  a  small  excursion 
into  the  forward  quadrant.  The  angular  coverage  can  be  increased  when  er  is 
made  larger,  but  we  must  take  care  that  the  second  mode  remains  below  cutoff. 
However,  the  angular  coverage  is  already  much  greater  than  that  obtained  with 
the  uniform  apertures.  Also,  when  the  strips  are  singly  periodic  as  indicated, 
the  broadside  angle  corresponds  to  a  stop  band,  and  scanning  in  that  narrow 
angular  region  must  be  omitted.  By  pairing  the  strips  a  quarter-guide- 
wavelength  apart  at  the  broadside  angle,  cancellation  can  occur  and  the  array 
can  be  scanned  through  broadside  as  well,  but  we  have  not  explored  that  option 
in  thfc  3tudy. 

An  important  advantage  of  this  array  is  that  the  radiation  into  space  has 
pure  horizontal  electric  field  polarization.  Another  important  advantage  is  that, 
with  the  baffles  present  (of  height  c  in  Fig.  5),  no  blind  spots  are  found  in  the 
scan  angle  behavior.  Without  any  baffles,  experience  with  other  arrays 
indicates  that  blind  spots  are  usually  present  when  dielectric  layers  are  part  of 
the  array  configuration.  The  short  baffles  appear  to  prevent  such  an  occurrence. 
In  addition,  they  separate  the  discontinuities  at  the  air-dielectric  interface  and  at 
the  open  upper  end,  making  the  analytical  result  simpler. 

The  attributes  of  this  novel  array  are,  therefore,  that  there  are  no  grating 
lobes,  no  cross  polarization,  and  no  blind  spots.  These  significant  advantages 
are  combined  with  simplicity  of  structure,  because  the  array  is  intended  for 
application  at  millimeter  wavelengths.  The  simplicity  requires  the  use  of  leaky- 
wave  line  sources,  which  may  not  be  suitable  for  certain  scanning  applications. 
The  important  new  feature  in  this  array,  as  compared  with  those  having 
uniform  apertures  [5,6],  is  that  the  appropriate  use  of  a  periodic  aperture 
permits  a  larger  and  more  flexible  scan  angle  coverage. 

C.  A  Leaky-Wave  Analysis  of  the  High-Gam 
Antenna  Configuration 

A  method  for  increasing  the  gain  of  a  printed-circuit  antenna  was 
presented  recently  [7]  which  involves  the  use  of  a  superstate  layer.  The 
antenna  element,  which  for  simplicity  was  taken  as  a  horizontal  electric  dipole, 
is  embedded  within  a  two-layer  geometry  consisting  of  a  bottom  (substrate) 
layer  of  dielectric  constant  €j  placed  on  a  metal  ground  plane  and  a  top 
(superstate)  layer  of  substantially  higher  dielectric  constant,  €2.  For  layer 
thicknesses  for  which  the  bottom  layer  is  effectively  a  half  "wavelength”  thick  in 
terms  of  the  transverse  (vertical)  wavenumber  in  that  medium,  and  effectively 
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a  quarter  "wavelength"  thick  in  a  similar  fashion  in  the  upper  layer,  a  narrow 
radiated  beam  about  any  desired  angle  0p  is  produced.  The  above-mentioned 
condition  was  called  a  "resonance  condition." 

In  [7]  it  was  shown  that  as  the  frequency  is  lowered  the  beam  scans  from 
endflre  to  broadside.  Asymptotic  formulas  were  derived  for  the  gain  and  the 
beamwidth  as  e2,  the  dielectric  constant  of  the  superstrate,  is  made  large.  In 
addition,  radiation  patterns  were  computed  numerically. 

The  present  study  involves  the  recognition  that  this  high-gain  effect  can  be 
interpreted  in  terms  of  leaky  waves.  More  precisely,  the  "resonance  condition” 
found  in  [7]  is  actually  a  good  approximation  to  the  transverse  resonance 
condition  for  the  existence  of  a  leaky  mode  on  the  structure.  This  leaky  mode 
travels  along  the  surface  radially  away  from  the  exciting  antenna,  and  it  leaks 
power  away  as  it  travels.  The  wavenumber  /?-  Ja  of  the  leaky  wave  along  the 
surface  defines  the  aperture  distribution  of  the  antenna,  and  determines  the 
shape  of  the  radiation  pattern.  The  leaky-wave  approach  furnishes  an 
alternative  interpretation  of  the  physical  effect. 

It  Is  not  necessary  to  employ  the  leaky-wave  approach  to  obtain  all  the 
numbers  for  the  far-fleld  performance,  but  this  approach  furnishes  insight  into 
the  physical  processes  involved,  provides  new  practical  information,  and  permits 
one  to  derive  the  same  results  as  those  reported  in  [7)  in  a  much  simpler  way. 
The  new  practical  information  includes  the  knowledge  of  the  field  behavior 
along  the  surface  in  a  simple  convergent  fashion.  That  information  allows  us  to 
determine  readily  the  size  of  the  supporting  structure  necessary  to  produce  the 
narrow-beam  radiation  patterns,  a  result  that  would  be  difficult  to  obtain 
otherwise. 

An  investigation  of  the  leaky  waves  excited  by  the  printed  antenna  shows 
that  two  leaky  waves  become  dominant,  a  TE-mode  leaky  wave  that  determines 
the  H-plane  radiation  pattern  and  a  TM-mode  leaky  wave  that  determines  the 
E-plane  pattern.  These  patterns  can  be  obtained  from  the  leaky-wave 
properties  in  a  simple  and  direct  way.  A  comparison  is  presented  in  Fig.  6  of 
the  H-plane  patterns  computed  directly  from  the  single  TE-mode  leaky  wave 
and  computed  exactly  numerically.  For  this  case,  the  beam  maximum  is  at  45°, 
and  the  dielectric  constants  of  the  substrate  and  the  superstrate  are  2.1  and 
10.0,  respectively.  It  is  seen  that  the  agreement  between  the  two  curves  is  very 
good,  verifying  that  the  field  along  the  surface  in  that  direction  is  indeed 
dominated  by  that  single  leaky  wave.  Contributions  from  the  other  leaky 
waves,  if  they  were  taken  into  account,  would  make  up  the  difference  between 
the  two  curves.  We  should  also  recognize  that  the  leaky  wave  result  is  not  only 
much  simpler  to  obtain  than  the  exact  one,  but  it  also  readily  exhibits  the 
physics  underlying  the  behavior. 

This  study  has  been  conducted  in  collaboration  with  Professor  D.R. 
Jackson  of  the  University  of  Houston.  The  results  obtained  from  the  first 
phase  have  been  presented  at  symposia  and  accepted  for  publication  [8-10]. 
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Comparison  of  calculations  for  the  H-plane  radiation  pattern  of  a 
high-gain  printed-circuit  dipole  antenna,  where  the  beam 
maximum  is  at  45°,  and  where  the  substrate  and  superstate 
dielectric  constants  are  2.1  and  10.0,  respectively.  The  solid  and 
dashed  lines  represent  the  exact  numerical  values  and  the  result 
obtained  from  a  single  leaky  wave  only. 
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D.  Conductor-Backed  Slot  Line: 

Dangers  and  Full-  Wave  Analyses 

It  is  very  tempting  in  microwave  integrated  circuits  to  introduce  a 
conductor  backing  for  slot  line  or  coplanar  waveguide  because  such  conductor 
backing  has  many  advantages.  Among  these  advantages,  as  quoted  in  the 
literature  [e.g.,  11,12],  are:  lowers  Z0,  lessens  dispersion,  improves  mechanical 
strength,  allows  easier  implementation  of  mixed  coplanar-microstrip  or  slot 
line-microstrip  circuits,  helps  in  grounding  floating  regions,  convenient  for  dc 
biasing,  and  provides  a  convenient  heat  sink.  There  are  also  certain  important 
dangers,  or  serious  potential  problems,  that  are  introduced  by  the  use  of 
conductor  backing.  These  potential  problems  include  leakage  of  power  into 
surface  waves  or  into  the  dielectric  region  between  the  plates,  unexpected  cross 
talk,  significant  alteration  of  the  guide  wavelength,  and  unexpected  or 
unwanted  coupling  to  neighboring  lines.  These  possible  unpleasant 
consequences  seem  not  to  have  been  discussed  in  the  literature. 

The  type  of  potential  problem  outlined  above  depends  on  the  lateral  extent, 
or  width  ,  of  either  the  conductor  backing  or  the  plates  comprising  the  slot  line 
or  the  coplanar  waveguide.  Several  examples  of  conductor-backed  slot  line, 
with  various  conductors  of  infinite  and  finite  lateral  extent,  are  presented  in 
Fig.  7;  the  upper  shielding  plate  or  cover,  which  is  customarily  there  but  plays  a 
minor  role,  is  omitted  from  Fig.  7  for  simplicity  (and  clarity). 

Structure  (a),  and  its  equivalent  in  coplanar  waveguide,  will  always  leak 
power  into  the  dielectric-filled  parallel  plate  regions,  and  the  leakage  rate  is 
rather  high,  as  we  find  quantitatively.  Structure  (b)  represents  the  other 
extreme,  a  conductor-backed  two-strip  line;  however,  it  is  also  a  pair  of  coupled 
microstrip  lines  in  its  odd  mode.  That  structure  and  those  in  (c),  (d),  (e)  and 
(h)  will  not  leak  at  lower  frequencies  but  may  leak  into  a  surface  wave  in  the 
outer  regions  at  higher  frequencies.  Ail  structures  with  finite  metal  lengths  that 
are  not  short,  such  as  (c),  (e),  (f)  and  (h),  face  the  danger  that  when  no 
leakage  occurs  the  guide  wavelength  can  be  seriously  altered  from  its  expected 
value  by  the  loading  produced  by  the  finite  transverse  line  lengths.  Structure 
(d),  which  involves  only  a  short  conductor  back,  and  which  is  utilized  in  the  de 
Ronde  coupler  [13],  will  not  face  this  danger.  Finally,  in  (g)  and  (h)  coupling 
(desired  or  not)  can  occur  between  parallel  slot  lines  not  near  each  other  on  the 
same  side  or  on  opposite  sides. 

The  only  full-wave  analyses  that  appear  in  the  literature  apply  to  the 
infinitely  wide  structure  in  (a)  [ ll]  and  the  other  extreme  in  (b).  The  analysis 
that  applies  to  (a)  was  actually  performed  for  coplanar  waveguide,  but  it  could 
Just  as  easily  have  been  done  for  slot  line.  However,  that  calculation  was  made 
several  years  ago,  before  it  was  recognized  that  leakage  could  occur;  that 
analysis  is  therefore  incomplete  since  it  misses  the  leakage  effect  altogether. 
Many  full-wave  analyses  for  the  structure  in  (b)  appear  in  the  literature  (see 
[14]  for  a  very  good  summary),  in  the  context  of  the  odd  mode  of  a  pair  of 
coupled  microstrip  lines.  None  of  those  treatments,  however,  inquired  as  to 
whether  or  not  leakage  into  a  surface  wave  could  occur  at  the  higher 
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frequences.  On  the  other  hand,  an  approximate  analysis  [  15]  was  performed 
for  the  two  strips  without  conductor  backing  which  does  ask  this  question,  and 
Indeed  answers  it  in  the  affirmative,  with  measurements  to  back  up  the 
conclusion.  To  our  knowledge,  there  are  no  full-wave  analyses  for  any  of  the 
other  situations. 

We  have  recently  developed  two  methods  of  analysis  that  are  applicable  to 
all  of  the  structures  in  Fig.  7,  thereby  providing  for  the  first  time  a  systematic 
full-wave  analytical  procedure  for  determining  the  effect  of  finite  plate  widths 
on  the  modal  propagation  characteristics.  One  method  is  based  on  mode 
matching  and  the  other  employs  a  new  transverse  equivalent  network  where  the 
elements  are  expressed  in  simple  closed  form.  Both  methods  were  applied  to 
the  infinitely  wide  structure  in  Fig.  7(a),  and  the  results  from  the  two  methods 
agreed  rather  well  with  each  other.  Measurements  were  also  made  for  both  the 
guidance  and  leakage  properties.  In  addition  to  verifying  that  the  conductor 
backing  leads  to  a  high  leakage  rate,  the  comparisons  with  measurement  and 
those  between  the  theories  show  that  the  theoretical  methods  are  reliable. 
These  studies  are  being  conducted  in  collaboration  with  Professor  H.  Shigesawa 
and  Dr.  M.  Tsuji  of  Doshisha  University,  Kyoto,  Japan. 

The  several  unexpected  dangers  that  can  arise  when  the  plates  are  finite  in 
lateral  extent  are  continuing  to  be  investigated  quantitatively.  Some  of  these 
effects  may  explain  why  designers  of  MMIC  or  microwave  integrated  circuits 
often  find  that  components  that  work  well  separately  fail  when  combined 
together. 
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1.  OBJECTTVE(S) 

This  investigation  aims  to  provide  a  new  and  systematic  approach  to  the 
propagation,  scattering,  guiding  and  coupling  of  beam  waves  in  open  layered 
structures,  particularly  those  of  the  thin-fllm  variety.  In  the  past,  problems 
involving  wave  progression  through  different  media  have  been  treated  mostly 
in  terms  of  plane  waves  of  infinite  extent,  or  in  terms  of  strongly  localized  rays 
that  connect  a  source  to  an  observation  point.  However,  an  increasing  number 
of  application  areas  utilize  fields  that  propagate  by  means  of  well-bounded 
beams,  such  as  those  provided  by  lasers.  Because  they  occupy  an  intermediate 
position  between  unbounded  plane  waves  and  rays  of  infinitesimal  cross- 
section,  beams  exhibit  propagation  characteristics  that  are  different  in  many 
respects  from  either  plane  waves  or  rays.  For  example,  they  are  subject  to 
nongeometric  lateral  displacements,  peculiar  focal  phase  shifts,  angular 
deviations,  anomalous  absorption  effects,  and  '  other  recently  recognized 
phenomena  that  have  not  as  yet  been  investigated  except  in  a  few  special  cases. 
It  is  thus  desirable  to  study  these  phenomena  and  other  related  aspects  in  terms 
of  beam  fields,  which  are  realistic  and  meaningful  representations  of  many 
practical  electromagnetic  (and  other)  fields. 

To  achieve  results  that  are  relevant  to  a  wide  spectrum  of  applications 
areas,  this  study  addresses  basic  aspects  of  beam  propagation  through,  and 
guidance  along,  open  layered  configurations  that  are  typical  of  a  very  broad 
variety  of  situations.  For  this  purpose,  we  consider  beams  of  the  Gaussian 
form  because  they  are  more  easily  tractable  from  mathematical  points  of  view 
and,  fortunately,  they  are  also  of  greatest  practical  application. 

The  configurations  of  interest  are  first  expressed  in  terms  of  canonic 
models  consisting  of  a  single  layer  placed  between  two  open  media  having 
different  electromagnetic  properties;  structures  with  multiple  layers  can 
thereafter  be  treated  as  composite  configurations  made  up  of  several  individual 
canonic  models.  While  the  investigation  has  focused  on  lossless  passive  media, 
the  projected  work  has  also  considered  media  with  arbitrary  absorption  (or 
other)  losses.  Some  of  the  specific  aspects  to  be  explored  include  beam-profile 
distortions,  the  interaction  of  beam  waves  with  other  (e.g.,  surface,  leaky) 
waves,  and  the  conversion  of  beam  fields  into  other  characteristic  field  types. 
Although  phrased  in  electromagnetic  terms,  the  results  of  this  basic 
investigation  are  applicable  also  to  other  areas,  such  as  acoustic  waves  along 
interfaces  or  layers,  elastic  waves  along  geophysical  strata,  plasmon  waves  in 
metal  films,  and  other  configurations  involving  waves  and  beams  through  or 


-  19  - 


SECTION  I:  ELECTROMAGNETICS 


along  layered  media. 


2.  SUMMARY  OF  RECENT  PROGRESS 

In  the  preceding  yearly  progress  report,  we  mentioned  that  our  studies  had 
included  three  phases,  namely  (a)  nonspecular  phenomena  in  beams  incident 
on  multilayered  media,  which  appear  as  peculiar  beam-shifting  and  other 
effects;  (b)  new  types  of  leaky  waves  supported  by  guiding  structures;  and  (c) 
surface  and  leaky  waves  along  general  open  configurations,  such  as  thin-fllm 
structures  for  integrated  optics,  microwave  acoustics  and  other  areas  of  current 
interest.  While  all  of  these  topics  were  pursued  during  the  past  year,  most  of 
the  work  during  that  period  has  focused  on  surface  plasmon  waves  along  thin 
metallic  films. 

In  our  previous  yearly  report,  we  have  discussed  that  our  work  on  leaky 
waves  along  canonic  structures  has  determined  all  the  possible  types  of  such 
waves  and  has  provided  a  detailed  catalog  of  their  properties.  We  found  that  a 
total  of  eight  leaky-wave  varieties  are  possible,  of  which  half  are  forward¬ 
traveling  and  the  other  are  backward-traveling.  However,  every  type  of  layered 
configuration  can  support  only  some  of  these  waves.  Because  the  (rather 
novel)  type  of  leaky  waves  having  backward-traveling  characteristics  can  occur 
only  in  the  presence  of  metallic  layers,  we  have  examined  such  (plasmon- type) 
structures  In  detail. 

In  particular,  we  have  explored  thin  metal  film  geometries  of  the  type  that 
permit  the  propagation  of  long  range  surface  plasmons  (LRSP)  waves.  In  this 
context,  a  closer  analysis  of  the  various  leaky  waves  that  may  be  supported  by 
such  geometries  has  revealed  novel  plasmon  waves  having  properties  that  are 
different  from  those  reported  so  far.  Specifically,  by  introducing  a  dielectric 
layer  in  parallel  to  the  metal  film,  we  have  shown  that  plasmon  waves  can  be 
generated  having  propagation  lengths  that  are  one  or  two  orders  of  magnitude 
larger  than  those  of  the  LRSP  variety.  We  have  therefore  denoted  this  new 
variety  by  the  name  of  extended  range  surface  plasmons  (ERSP)  waves.  In 
addition,  we  have  determined  the  optimal  conditions  for  exciting  such  ERSP 
waves  by  means  of  incident  beams  having  Gaussian  or  rectangular  profiles. 
The  results  obtained  show  that  ERSP  waves  can  provide  fields  with  greatly 
enhanced  intensities  over  very  long  distances.  These  properties  are  very  useful 
for  applications  that  require  strong  field  Interactions  with  nonlinear  media. 

Papers  covering  the  above  topics  were  written  and  presentations  were 
made  at  several  technical  meetings. 
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3.  STATE  OF  THE  ART  AND  PROGRESS  DETAILS 

A.  Background 

The  interest  in  beam  fields  was  stimulated  by  applications  in  laser 
technology,  fiber  optics,  ultrasonic  diagnostics  and  other  areas  which  have 
generated  a  number  of  fundamental  studies  of  beam  propagation  through  a 
variety  of  media.  However,  most  of  the  studies  have  dealt  with  beams 
propagating  over  long  distances  through  uniform  or  "slowly  varying"  media. 
Thus,  comparatively  few  investigations  have  been  carried  out  on  beam  fields 
interacting  with  abrupt  transitions  and/or  multiple  boundaries  separated  by 
distances  of  the  order  of  a  wave-length,  such  as  occur  in  thin  film 
configurations. 

Because  lasers  generate  beam  fields,  a  considerable  number  of 
investigations  on  beams  inside  laser  cavities  were  undertaken  in  the  early 
sixties.  These  studies  viewed  the  beams  primarily  in  terms  of  ray  optics  and 
considered  their  interaction  with  mirror  reflectors  or  thin  lenses  [ l ,2] .  The 
fields  of  beams  that  are  emitted  by  lasers  into  free  space  have  also  been 
studied,  mostly  in  terms  of  their  Hermite-Gaussian  modes  [1-3],  and  interest 
into  the  fine  structure  of  those  propagating  fields  is  continuing  to  this  day  [4,5]. 
However,  all  of  those  studies  are  concerned  primarily  with  a  beam  that 
propagates  basically  through  a  single  uniform  medium. 

When  the  medium  inside  which  the  beam  propagates  is  not  uniform,  the 
usual  assumption  has  been  that  the  medium  properties  change  slowly  over  a 
distance  of  a  wavelength,  i.e.,  one  deals  with  a  slowly  varying  medium.  The 
treatment  has  then  involved  a  variety  of  techniques,  which  include  complex 
rays  [6,7],  or  the  propagating- be  am  method  [8],  which  relies  heavily  on 
numerical  computations. 

In  contrast  to  the  above,  our  study  is  concerned  with  beam  fields  that 
propagate  through  media  that  can  change  rapidly  and  are  therefore 
characterized  best  by  boundaries  that  separate  media  having  different 
electromagnetic  properties.  In  this  context,  most  recent  studies  have  dealt 
almost  exclusively  with  the  lateral  Goos-Haenchen  shift  of  beams  incident 
under  total  reflection  conditions  at  a  single  interface  between  different  media 
[9-15].  Only  a  very  few  investigations  have  examined  incidence  of  beams  on 
layered  configurations  having  two  or  more  interfaces  [16-19] .  Of  these,  a 
considerable  amount  of  basic  work  was  carried  out  by  us  [  17,19];  in  this 
context,  we  have  also  shown  that  electromagnetic  beam  phenomena  on  layered 
configurations  have  analogous  counterparts  in  acoustics  [20,21],  thus 
stimulating  investigations  in  that  area  which  are  still  continuing  [22-34]. 

In  addition  to  the  (Goos-Haenchen)  lateral  displacement  of  the  beam  axis, 
recent  studies  have  shown  that  the  beam  may  also  exhibit  an  angular  deviation 
[11,26,27],  and/or  a  focal  shift  [24,25].  However,  those  effects  were  examined 
only  for  beams  incident  from  a  denser  medium  onto  a  single  interface  to  a  rarer 
medium  under  total- re  (lection  conditions.  By  extending  an  analytical  pole-zero 
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approach  recently  developed  by  us  [33],  we  have  obtained  qualitative  and 
quantitative  results  for  all  of  these  effects  for  arbitrary  incidence  on 
multilayered  structures  [34].  This  approach  has  revealed  a  hitherto  unknown 
phenomenon,  which  manifests  itself  as  a  reduction  or  increase  in  the  effective 
beam  width.  The  unified  analysis  of  the  three  (lateral,  focal  and  angular) 
beam-shifting  effects  and  the  fourth  (beam-waist  modification)  phenomenon 
were  described  by  us  in  previous  progress  reports  and  subsequently  published 
[34-38]. 

B.  Beam  Excitation  of  Plasmon  Waves  Along  Metal  Films 

As  shown  by  us  in  the  past  [31-38],  the  nonspecular  effects  described 
above  are  intimately  related  to  the  excitation  of  leaky  waves  which  can  be 
supported  by  the  pertinent  layered  structures.  As  part  of  the  effort  during  the 
preceding  year,  we  have  therefore  devoted  a  large  portion  of  our  studies  to 
explore  and  classify  the  various  varieties  of  leaky  waves  that  can  be  guided  by 
general  layered  media  of  the  type  described  in  Fig.  1.  The  results  of  that  study 
were  reported  in  the  previous  annual  report,  and  most  of  the  relevant  material 
has  also  been  published  [39,40,41,43].  While  we  have  continued  these  general 
studies  of  leaky-wave  fields,  we  have  concentrated  during  the  past  year  on  the 
canonical  configurations  described  in  Fig.  1  for  the  specific  case  where  the  layer 
shown  therein  by  e(z)  contains  multiple  lossy  or  lossless  dielectric  and/or  metal 
layers. 

Because  strong  interest  has  recently  been  shown  [42,43]  in  plasmon  waves 
guided  by  thin  metal  films,  we  have  specifically  focused  on  the  particular  leaky 
waves  that  are  responsible  for  the  plasmon  modes.  The  general  situation  we 
examined  is  then  shown  in  Fig.  2,  wherein  a  beam  is  assumed  incident  in  a 
prism  placed  upon  a  multilayered  geometry  which  includes  a-metal  film.  The 
(leaky-wave)  plasmon  modes  excited  by  the  incident  beam  are  then 
characteristic  solutions  of  the  pertinent  boundary- value  problem,  with  fields  in 
the  form 


F  =  Fjexp  ik(«x  +  t jz)  -  iwt 

•  a 


where  J  denotes  the  specific  layer  under  consideration  and  k=27r  /X.  Here  k  and 
t j  are  normalized  wavenumbers  which  satisfy 

« 2  +  r2  =  €j  ,  (2) 

and  Cj  refers  to  the  dielectric  constant  of  the  J1*  layer.  The  wavenumbers  k  and 
t j  are  found  from  Eq.  (2)  in  conjunction  with  a  dispersion  (secular)  relation  in 
the  form 

D(k)  =  f(K:r,)  =  0  .  (3) 
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Z 


CANONIC  CONFIGURATIONS  FOR  e(z)  =  eh  =  const. 


Duct: 


«s  >  «h  >  eu  >  0 


Transition:  es  >  «h  >  *u  >  0 


Gap: 
Plasma : 


>  e„  >  eK  >  0 


fis  >  «u  >  0  >  eh 


Open  layered  structure  supporting  leaky  waves  and  the  canonic 
configurations  that  characterize  typical  situations. 


Long-Range 

Surface-Plasmon 

(LRSP) 

Geometry 


Dielectric  layer  for 
Extended-Range 
Surface  Plasmon  (ERSP) 


Geometry  for  exciting  plasmon  waves  by  means  of  an  incident  beam 
of  width  2w.  In  the  absence  of  the  dielectric  layer  shown  by  dashed 
lines,  the  configuration  corresponds  to  Sarid's  structure  [42]  for  LRSP 
waves.  The  presence  of  the  dielectric  layer  indicated  by  dashed  lines 
permits  the  excitation  of  ERSP  waves  having  much  longer  propagation 
ranges. 
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By  solving  Eq.  (3)  under  certain  restricted  conditions,  Sarid  has  shown 
[42]  the  possibility  of  exciting  plasmon  waves  having  long  propagation  lengths, 
which  were  subsequently  labelled  long  range  surface  plasmons  (LRSP)  waves. 
We  have  verified  the  presence  of  such  LRSP  waves  and  have  shown  [44]  that 
they  can  be  more  strongly  excited  if  the  angle  of  the  incident  beam  is  offset 
from  the  angle  that  satisfied  a  reflection  dip.  (This  reflect, ion-dip  condition  had 
in  the  past  been  assumed  to  provide  the  best  excitation  parameters.)  More 
interestingly,  however,  we  have  solved  the  dispersion  equation  (3)  under  more 
general  conditions  and  have  found  that  plasmon  waves  can  be  generated  with 
propagation  lengths  that  are  considerably  longer  than  those  of  the  previously 
known  LRSP  waves. 

To  clarify  these  novel  plasmon  wave  types,  consider  Fig.  2  in  which  the 
conventional  LRSP  is  indicated  by  solid  lines.  This  consists  of  the  prism,  a  gap 
of  thickness  tlf  a  metal  film  of  thickness  tm  and  a  substrate  having  a  dielectric 
constant  ea.  To  obtain  the  larger  propagation  ranges,  an  additional  dielectric 
layer  having  a  dielectric  constant  es>ed  is  introduced,  as  shown  by  the  dashed 
lines  in  Fig.  2.  The  function  of  the  additional  layer  is  to  serve  as  a  planar 
waveguide  whose  field  can  couple  to  the  surface  plasmon  guided  by  the  metal 
film. 

We  have  thus  shown  [45]  that  this  coupling  effect  provides  a  mechanism 
whereby  the  additional  layer  acts  as  an  energy  reservoir  which  increases  the 
power  storage  of  the  system.  As  a  result,  the  net  propagating  plasmon  wave 
can  progress  over  a  distance  Le  which  is  more  than  one  order  of  magnitude 
larger  than  the  propagation  range  L0  of  the  (ordinary)  LRSP  waves  known  so 
far  [42,43].  To  distinguish  between  the  two  plasmon  waves,  we  label  the 
plasmon  having  the  larger  (extended)  range  as  an  extended  range  surface 
plasmon  (ERSP)  wave. 

The  importance  of  the  ERSP  wave  is  illustrated  in  Fig.'  3(a),  where  the 
ratio  Le/L0  between  the  propagation  ranges  mentioned  above  is  shown  by  the 
solid  line.  As  a  function  of  the  separation  between  the  metal  film  and  the 
auxiliary  dielectric  layer,  it  is  clearly  possible  to  obtain  ratios  Le/Lo=l00  and 
larger.  In  Fig.  3(b),  the  power  density  ratio  Pxe/Pxo  between  the  ERSP  and 
LRSP  is  shown  by  the  solid  line.  This  ratio  is  important  for  field  interactions 
with  non-linear  media  because  the  power  density  Px0  of  the  LRSP  wave  is  very 
large  in  the  metal  film  compared  to  the  power  flux  of  the  incident  beam  that 
excites  the  plasmon  field.  As  can  be  seen  from  Fig.  3(b),  the  ratio  Pxe/Pxo  is 
unity  or  larger,  so  that  the  excitation  of  the  ERSP  wave  can  provide  a 
considerably  longer  propagation  range  at  no  reduction  in  field  enhancement. 

To  complete  the  picture,  we  have  investigated  the  presence  of  small  losses 
in  the  dielectric  media.  The  results  are  shown  by  the  dotted  lines  in  Fig.  3. 
which  show  that  the  propagation  range  Le  can  be  about  100  times  larger  than 
L0  at  Ucsil.5^m.  Under  those  conditions,  Pxe/Pxo  —  0.85,  which  implies  that 
the  field  intensity  enhancement  achieved  by  the  ERSP  wave  is  only  lb%  lower 
than  that  of  the  LRSP. 


SECTION  I:  ELECTROMAGNETICS 


As  a  result,  we  And  that  the  plasmon  waves  of  the  ERSP  variety  in  realistic 
geometries  can  have  propagation  ranges  of  the  order  of  1  cm  or  larger  with  an 
attendant  field  enhancement  that  is  only  slightly  smaller  than  that  of  the  much 
shorter-propagating  LRSP  waves.  The  excitation  of  ERSP  waves  therefore 
produces  remarkably  improved  conditions  over  long  distances  in  non-linear 
regimes. 

Publications  concerning  the  above  effects  have  been  already  made  [44-46] 
and  a  Doctoral  thesis  dealing  with  the  detailed  work  is  presently  being  written 
[47]. 
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C.  MIXED  SPECTRAL  TECHNIQUES  FOR  WAVE  PROPAGATION 
AND  DIFFRACTION 

Professor  L.B.  Felsen 

Unit  EM7-3 


1.  OBJECTTVE(S) 

High  frequency  or  transient  propagation  in,  or  transmission  through,  lay¬ 
ered  media,  and  high  frequency  or  transient  scattering  by  impenetrable  and 
penetrable  targets  usually  requires  synthesis  in  terms  of  a  large  number  of  basic 
wave  processes.  For  the  guiding  or  ducting  problem,  these  wave  processes  are 
either  normal  (discrete  and  continuous)  modes  or  ray-optical  fields.  For  the 
transmission  problem,  the  basic  wave  processes  are  traveling  waves  which 
undergo  multiple  internal  reflection  at  the  layer  boundaries.  For  the  transient 
scattering  problem,  the  wave  processes  involve  multiple  wavefront  and  reso¬ 
nance  fields.  Because  descriptions  by  multiple  propagation  events  are  often 
poorly  convergent  and  do  not  provide  physical  interpretation  in  compact  form, 
it  is  desirable  to  seek  collective  descriptions  of  multiple  phenomena. 

Thus,  the  objective  of  this  fundamental  study  is  the  construction  of  a  new 
theory  of  propagation,  transmission  and  scattering  that  has  broad  implications 
for  a  general  class  of  electromagnetic  and  other  wave  problems.  The  approach 
is  to  seek  suitable  spectral  representations  in  the  spatial  and  temporal  domains 
and  to  adapt  these  by  arguments  of  asymptotic  localization  to  successively  more 
complicated  environments.  The  basic  spectral  building  blocks  are  time- 
harmonic  and  transient  local  plane  waves,  and  the  spectra  may  be  real  or  com¬ 
plex.  Localized  spectra  can  represent  ray  fields,  true  and  local  mode  fields, 
hybrid  ray-mode  forms,  wavefront-resonance  phenomena,  and,  in  more 
fleshed-out  form,  also  transitional  effects  where  these  simple  compact  spectral 
formulations  may  fail.  The  tools  involve  the  theory  of  spectral  representations 
and  asymptotic  treatment  of  integrals  and  partial  differential  equations. 


2.  SUMMARY  OF  RECENT  PROGRESS 

This  section  presents  a  brief  summary  of  recent  progress;  more  detailed 
descriptions  of  selected  portions  are  contained  in  the  next  section. 
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A.  Time-Harmonic  Fields 

(1)  Wave  Coupling  into  Large  Enclosures 

One  phase  of  our  investigations  aimed  at  understanding  propagation  in 
rather  general  environments,  especially  at  high  frequencies,  has  been  concerned 
with  wave  coupling  into  enclosed  guiding  regions  via  openings  in  the  enclosure. 
Our  goal  here,  as  in  all  of  our  model  studies,  is  to  arrive  at  a  "good”  parameter¬ 
ization  of  complicated  propagation  and  diffraction  phenomena.  Good  parame- 
terizations  are  stable  under  perturbation  and  can  therefore  be  extended  from  an 
idealized  prototype  to  more  general  configurations  belonging  to  the  same 
phenomenological  class.  From  this  perspective,  we  have  examined  in  detail  the 
coupling  of  a  plane  wave  into  a  large  open-ended  parallel  plane  waveguide. 
Although  this  problem  is  classical,  we  have  asked  new  questions,  especially 
within  the  framework  of  our  self-consistent  hybrid  ray-mode  format,  which  is  a 
crucial  constituent  in  the  good  parameterization  scheme.  Selected  results  are 
summarized  in  Section  3.  A  similar  study  of  plane  wave  coupling  into  an 
open-ended  circular  waveguide  is  nearing  completion. 

(2)  Complex  Rays  and  Beams 

The  virtues  and  limitations  of  complex  ray  and  paraxial  Gaussian  beam 
(GB)  theory  continue  to  be  explored  for  pertinent  model  problems,  with 


emphasis  on  questions  relating  to  algorithm  implementation,  and  also  wiih 
respect  to  a  novel  usage  of  GB’s  as  basis  elements  in  rigorous  field  representa¬ 
tions  spanning  a  discretized  (configuration)^ wavenumber)  phase  space  [  1,2] .  A 
detailed  prototype  study  of  the  complex  ray  method  for  beam  reflection  near 
critical  incidence  has  been  completed  and  published  [3]. 


B.  Transient  Fields 


Our  general  spectral  theory  of  transients  (STT)  [4j,  which  represents  tran¬ 
sient  wave  phenomena  as  superpositions,  over  spatial  wavenumber  spectra,  of 
transient  plane  waves  Instead  of  the  usual  decomposition  into,  and  reconstruction 
via,  time-harmonic  wave  fields,  has  been  applied  to  pulsed  input  signals,  espe¬ 
cially  those  generating  pulsed  focused  beams.  The  analysis  for  one  such  input, 
the  complex  source  pulsed  beam,  is  described  in  more  detail  in  Section  3.  The 
analysis  for  another,  the  focus  wave-mode,  is  nearing  completion. 
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C.  Effects  of  Random  Media 

This  new  phase,  which  was  introduced  in  the  previous  progress  report  and 
which  has  been  supported  by  the  Rome  Air  Development  Center  under  a  spe¬ 
cial  fund  (but  administered  through  JSEP),  addresses  the  modification  of  high 
frequency  propagation  and  diffraction  phenomena  when  the  ambient  propaga¬ 
tion  medium  undergoes  weak  large-scale  random  fluctuations.  Such  effects  are 
important  for  millimeter  and  optical  wave  propagation  in  the  atmosphere,  and 
in  other  environments  with  properties  so  complicated  as  to  require  statistical 
treatment.  Our  aim  has  been  to  construct  a  new  stochastic  geometrical  theory 
of  diffraction  (SGTD)  which  attempts  to  merge  the  constructs  of  deterministic 
GTD  with  the  stochastic  modifications  introduced  by  the  fluctuating  medium. 

Formulation  of  the  theory  has  been  completed,  and  we  have  shown  how  it 
would  be  applied  to  reflection  from,  and  refraction  through,  curved  interfaces, 
to  edge  diffraction,  and  to  diffraction  by  small  scatterers  embedded  within  the 
fluctuating  environment  [5-7].  A  comprehensive  manuscript  is  in  preparation. 
The  theory,  which  embodies  approximations,  must  now  be  tested  by  numerical 
experiments  and  comparisons  with  other  solutions  (if  they  exist)  or  with  exper¬ 
imental  data.  This  work  is  presently  in  a  (hopefully  temporary)  state  of  suspen¬ 
sion  because  Dr.  R.  Mazar,  a  visiting  post-doctral  scientist  who  provided  the 
expertise  in  random  propagation  theory,  has  now  returned  to  his  home  country 
(Israel). 

3.  STATE  OF  THE  ART  AND  PROGRESS  DETAILS 

A.  Background 

Many  electromagnetic  propagation  environments,  whether  natural  or  man¬ 
made,  are  so  complicated  that  direct  solution  of  the  field  equations  to  deter¬ 
mine  signal  characteristics  is  beyond  the  scope  of  present  analytical  and  compu¬ 
tational  capabilities.  At  high  frequencies,  propagation  can  be  localized  and 
approximated  as  ray  fields  which  undergo  reflection,  refraction  and/or 
diffraction  on  their  path  from  a  source  at  S  to  an  observer  at  P.  While  ray 
theory  provides  a  fundamental  view  of  the  propagation  process  by  tracking  local 
plane  wave  fields  emanating  at  the  source,  such  tracking  becomes  cumbersome 
when  many  ray  paths  exist  between  S  and  P.  It  would  therefore  be  desirable  to 
deal  with  multipath  effects  collectively.  In  guided  propagation  along  a  refracting 
channel,  rays  may  form  caustics  (convergence  or  focusing  zones  of  enhanced 
field  strength)  where  simple  ray  theory  fails.  When  these  caustics  are 
sufficiently  distinct,  one  may  correct  ray  theory  by  uniform  asymptotic  transi¬ 
tion  functions,  but  situations  arise  for  rays  with  many  reflections  where  an 
accumulation  of  caustics  makes  such  corrections  impractical  and  even  impossi¬ 
ble.  Here,  again,  a  collective  alternative  to  multiple  ray  reflections,  for  exam¬ 
ple,  by  employing  guided  modes,  is  desirable  or  actually  necessary.  However, 
the  modal  approach  to  ducted  propagation  is  beset  with  similar  difficulties  when 
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the  required  number  of  modes  is  large.  It  would  then  be  useful  to  express  the 
Interference  properties  of  clusters  of  modes  collectively  in  terms  of  simpler 
events. 

The  preceding  discussion  makes  evident  the  importance  of  collective  treat¬ 
ment  of  mode  or  ray  fields  when  many  of  these  are  required  to  synthesize  the 
signal  in  a  particular  transmission  or  guided  propagation  channel,  or  when 
failures  in  approximate  mode  or  ray  theory  make  these  descriptions  inapplica¬ 
ble.  Substantial  progress  in  this  direction  has  been  made  through  the  new 
hybrid  ray-mode  theory  developed  by  us.  It  has  been  shown  how  clusters  of  ray 
fields  excited  by  a  localized  source  can  rigorously  be  converted  into  clusters  of 
guided  mode  fields  plus  a  (usually  small)  remainder,  and  vice  versa.  The 
theory  has  been  applied  to  a  series  of  "canonical"  problems  involving  guided 
electromagnetic  propagation  along  concave  surfaces  (here,  the  guiding  mechan¬ 
ism  is  provided  by  "whispering  gallery”  effects),  in  tropospheric  ducts,  in  plane 
parallel  homogeneously  filled  waveguides,  and  in  graded  index  waveguides. 
The  theory  has  also  had  impact  on  other  fields  such  as  underwater  acoustic  pro¬ 
pagation  and,  with  generalization  to  time-dependent  signals,  the  modeling  of 
seismic  events.  Concern  in  these  applications  has  been  with  the  greater  compu¬ 
tational  efficiency  of  the  hybrid  formulation,  with  the  avoidance  of  singular 
regions  in  ray  fields  or  mode  fields  by  filling  these  regions  with  mode  fields  or 
ray  fields,  respectively,  and  also  with  the  penetrating  physical  insight  of  the  pro¬ 
pagation  mechanism  provided  by  the  hybrid  method.  The  collective  approach 
linking  ray  fields  and  mode  fields  has  been  applied  also  to  the  complex  spectra 
emitted  by  a  source  at  a  complex  location.  Such  a  source  generates  in  physical 
space  a  Gaussian  beam,  and  therefore  makes  the  important  class  of  beam  pro¬ 
pagation  and  diffraction  problems  amenable  to  the  hybrid  format.  When  a 
guiding  environment  exhibits  weak  variations  along  the  guiding  direction,  the 
constructs  above  are  generalized,  by  localization,  to  accommodate  adiabaticaliy 
adaptive  wave  phenomena.  This  has  been  formalized  in  our  intrinsic  mode  and 
scaled  spectrum  theory.  All  of  these  accomplishments  have  been  documented 
in  a  comprehensive  series  of  publications,  cited  in  previous  reports. 

B.  Selected  Results 

(1)  Rays,  Modes  and  Beams  for  Plane  Wave  Coupling  into  a  Wide 
Open-Ended  Parallel- Plane  Waveguide 

A  plane  wave  incident  obliquely  onto  the  open  end  of  a  wide  parallel-plane 
waveguide  generates  an  interior  field  that  resembles  a  geometric  optical  sheet 
beam  whose  initially  well-defined  shape  gets  diffused  after  undergoing  multiple 
reflection.  This  behavior  can  be  confirmed,  but  not  physically  well  explained, 
from  numerical  results  obtained  by  summing  over  the  aperture-excited  guided 
modes.  A  phenomenologically  incisive  representation  (a  good  parameteriza¬ 
tion)  is  obtained  by  subjecting  the  guided  mode  sum  to  partial  Poisson  summa¬ 
tion  that  yields  an  exact  hybrid  field  comprising  the  multiple  reflected  geometric 
optical  beam,  some  multiple  reflected  edge-diffracted  rays,  some  guided  modes. 
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plus  remainder  terms.  Asymptotic  approximations  reduce  these  exact  forms  to 
the  wave  constituents  of  the  nonuniform  and,  when  required,  uniformized 
geometrical  theory  of  diffraction.  These  basic  features  are  illustrated  by  numer¬ 
ical  r°s"!ts  presented  for  various  parameter  regimes.  They  delineate  the  ranges 
of  applicability  of  the  asymptotic  forms  and  also  the  relative  importance  of  the 
individual  wave  contributions.  In  the  parameter  range  where  geometrical 
features  of  the  field  can  still  be  resolved,  the  hybrid  form  is  not  only  physically 
appealing  but  also  numerically  more  efficient  than  alternatives  involving  only 
rays  or  only  modes. 

The  physical  configuration  is  shown  in  Fig.  1,  which  also  schematizes  the 
relevant  ray  phenomena  (Fig.  la),  the  physical  optics  (PO)  method  employed 
for  calculation  of  the  reference  solution  by  modal  summation,  and  the  decom¬ 
position  of  a  typical  guided  mode  into  four  modal  ray  congruences.  The  edge 
diffracted  ray  fields  are  decomposed  Into  four  ray  species  in  a  similar  manner. 
The  hybrid  ray-mode  decomposition  is  schematized  in  Fig.  2.  It  shows  the 
multiply  reflected  geometric  optical  (GO)  beam  and  also  the  multiple  reflected 
edge  diffracted  ray  fields,  which  can  be  taken  to  arrive  from  an  array  of  image 
sources  (Fig.  2a);  the  total  field  calculated  in  this  manner  can  be  converted  into 
the  complete  set  of  guided  modes.  For  the  hybrid  form,  only  the  relevant  ray 
fields  are  converted  (Fig.  2b);  their  departure  angles  lie  in  the  shaded  sections. 
Figures  3  and  4  show  typical  results  obtained  by  the  strategies.  Details  may  be  ■ 
found  in  reference  8. 

(2)  Complex-Source  Pulsed  Beams 

By  assigning  complex  values  to  the  source  coordinates  and  pulse-initiation 
time  of  the  time-dependent  Green's  function  in  free  space,  one  may  generate  a 
new  exact  field  solution  that  behaves  like  a  propagating  pulsed  beam.  Although 
the  conventional  pulsed  line  source  response  Is  known  in  closed  form,  the  com¬ 
plex  extension  cannot  be  performed  directly  thereon  because  of  the  nonanalytic 
behavior  of  the  causal  field.  The  analytic  continuation  has  been  carried  out  by 
spectral  analysis  and  synthesis,  utilizing  the  recently  formulated  spectral  theory 
of  transients  (STT)  [4].  This  approach  not  only  guarantees  uniqueness  but  also 
elucidates  the  spectral  content  of  the  resulting  waveform,  which  is  composed  of 
contributions  from  singularities  in  the  complex  spectral  wave-number  plane. 
The  new  pulse  solution  is  of  interest  not  only  in  free  space  but  also,  by  similar 
analytic  extension  of  time-dependent  Green’s  functions  in  more  complicated 
environments,  for  the  direct  construction  of  the  transient  field  produced  in 
these  environments  by  the  incident  pulsed  beam. 

We  exemplify  construction  of  the  solution  for  the  two-dimensional  case  in 
a  g  —  (x,z)  coordinate  system.  The  conventional  time-dependent  Green's  func¬ 
tion  for  real  (g',t'),  where  g'=(x',z')  is  the  source  location  and  t'  is  the  source 
initiation  time,  is  known  to  be  given  in  closed  form  by 

G(g.g,.t,t/)  =  H  ^t-t')  -  s/v  j  ,  A=  [(t- t')‘ -  (s/v)-  .  (i) 
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Edge  diffraction 
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Plane  wave  excitation  of  open-ended  parallel-plate  waveguide. 

(a)  Physical  configuration  and  GTD  schemattzatton  of  waveguide 
mode  excitation. 

(b)  PO  formulation  in  terms  of  equivalent  magnetic  currents  M  on 
a  perfect  conductor  closing  the  waveguide  aperture  at  0<x<a,  z=o. 

(c)  Modal  congruences  J=i . 4  and  their  shadow  boundaries  SB 
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GO  beam  due  to  all  species,  and  image  construction  of  multiple 
reflected  edge-diffracted  rays  for  ray  species  J=  1.  SBj  denotes  the 
shadow  edge  of  the  multiple  reflected  truncated  plane  wave;  the  GO 
ray  contribution  at  P  due  to  this  wave  is  shown  as - . 

(a)  Full  conversion  of  modes  to  rays  through  Poisson  sum  formula; 
all  multiple  reflected  rays  are  required. 

(b)  Partial  conversion  of  modes  um,  Mx <m <M2,  to  rays;  only  those 
rays  whose  departure  angles  0{l*  lie  In  the  shaded  interval 

are  necessary. 


kz=*204 


Field  magnitude  along  cross  section  kx=204  in  a  waveguide  with 
normalized  width  ka=l007r,  where  k  is  the  wavenumber.  100  modes 
can  propagate  here.  Plane  wave  incidence  angle  0j=3O°.  The  refer¬ 
ence  solution  (solid  curve)  is  generated  by  mode  summation.  The 
rectangular  trace  outlines  the  GO  beam. 
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Detailed 
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fields  over  portions  of  the  cross  section  in 
are  indicated.  Solid  curve:  100  mode : 
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Here  H(a)=l  or  H(a)=0  for  a>0  or  a<0,  respectively,  v  is  the  wave  propaga¬ 
tion  speed  in  the  (nondispersive)  medium,  and 


s 


1  1/2 

(x-x')2  +  (z  -  z')2 


(2) 


is  the  distance  from  £'  to  the  observer  at  £ . 

When  the  source  location  £'  is  made  complex,  direct  substitution  into  the 
closed-form  expression  for  the  time-harmonic  free-space  Green’s  function 
G(£,£/,o;)  is  known  to  generate  at  real  observation  points  £  a  closed-form  field 
that  behaves  like  a  Gaussian  beam.  The  proper  analytic  continuation  of  the 
complex  distance  s  in  Eq.  (2)  from  real  to  complex  values  is  easily  established 
in  this  case.  However,  for  the  causal  transient  field  in  Eq.  (1),  the  interpreta¬ 
tion  of  complex  £'  and  (or)  t'  substitution  is  less  obvious.  The  problem  can  be 
resolved  by  the  spectral  analysis  and  synthesis  route  of  STT.  We  choose 

x'  =  0,  z'  =  ib,  b  >  0,  (3) 

which,  under  time-harmonic  conditions,  is  known  to  generate  a  beam  field  with 
maximum  along  the  z  axis,  waist  at  z=  0,  and  beam  width  determined  by  b 
(Fig.  5).  For  the  time-dependeflt  case,  we  also  assign  the  complex  value 

i  =  ib/v  (4) 


Fig.  5  The  z-directed  beam  generated  by  a  time-harmonic  line  source  at  the 
complex  location  £ '=(x',z')=(0,ib).  The  complex  source  location  is 
schematized  by  a  dashed  line.  The  parameter  b  is  related  to  the 
beam  width  at  the  waist  by  W(0)=(b/k)l/2.  At  any  location  z.  the 
beam  width  is  W(z)=W(o)[i+(z/b)2] l/2. 


to  the  source-initiation  time. 
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By  STT,  the  source-excited  transient  held  is  synthesized  as  a  spectrum  of 
source-excited  transient  plane-waves  A(£)/(t-  t'-  r(£)],  where  f  is  the  spatial 
wavenumber  corresponding  to  the  x-coordinate.  Explicitly  [9] 

G  =  Re  G  +  (5) 


where 


G+(t) 


s  - 

l(t) 


JL  Mil 

V  r'(f) 


and 


00 


A (O 


A(D 


t-t'-r(£)  t-t'-r(f) 


dC  - 


r(£,£/;^)  *  —  [(x-x')£±  (z-zOVl-  £2  ]  ,  z  <  z' , 


(6) 


(7a) 


A(£,g';f)  =  -  1  ■  .  (7b) 

2Vl-e 

The  integral  in  (6)  may  be  evaluated  in  terms  of  the  singularities  of  the 
Integrand  in  the  complex  £  plane.  These  involve  the  singularities  of  A  and  r, 
and  the  time-dependent  poles  |(t) 

r[|(t)]=t-t'  (8) 


that  describe  the  transient  local  plane  waves  arriving  at  the  observer  at  time  t, 
The  two  roots  of  Eq.  (8)  are 

«*  =  [(> 


where 


and 


with 


;x-x')(t-t')  ±  i(z-z')A 


v/s2 


1/2 


A  =  ^t2-  |— |  -  2l~  jt- —  costfj  J 

-  7r  <  argA  <  0 


s  =  (/?'-  b“-  2ibp  cos#) 


Re  s  >  0 


1/2 


(9) 


(10) 

(11) 

(12) 


Here,  p=(x2+z2)1^2  and  0=sin"  l(x/p).  The  integral  can  now  be  reduced  to 
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the  closed  form 


G(g  .g'jt.t')  =  Re 


1 

27T  A 


1  +  —  log 

7Ti 


tr-  t^+ A  | 
s/v  )  _ 


(13) 


with 


-37T/2  <  Imlog  <  7T/2  .  (13a) 

This  result  constitutes  the  unique  analytic  continuation  of  (l)  from  real  to  com¬ 
plex  source  coordinates. 

The  waveform  generated  in  this  manner  has  the  characteristics  of  a  pulsed 
beam  in  the  z>0  region,  which  is  of  concern  here  (see  Fig.  Q).  The  field  peaks 
at  approximately  t=p/v,  preceded  by  a  weak  precursor.  The  strongest  peaks 
occur  near  the  z-axis  and  reach  a  maximum  on  axis.  Additional  details  pertain¬ 
ing  to  the  physical  behavior  of  this  pulse,  as  well  as  pulse  shapes  generated  by 
other  choices  of  parameters  and  alternative  analytic  continuation  of  the  com¬ 
plex  distances,  may  be  found  in  reference  9.  The  three-dimensional  version  of 
the  complex  space-time  source  coordinate  substitution  has  been  explored  as 
well. 


Pulsed-beam  waveform  at  various  observation  points  z/b=0.l,  near 
the  initial  plane,  and  in  the  far  zone  p/ b=2.  The  waveform  is  nor¬ 
malized  to  (b/v)G.  The  normalized  time  (v/b)t  and  space  coordi¬ 
nates  (x/b,  z/b)  are  drawn  to  the  same  scale. 


Fig.  0 


SECTION  I:  ELECTROMAGNETICS 


The  pulsed  beams  described  above  are  unique  in  that  they  represent  ana¬ 
lytic  continuations  of  Green’s  functions.  It  is  therefore  suggestive  that  they 
may  be  useful  for  synthesis  of  more  general  source  distributions  designed  to 
create  well-focused  fields.  This  is  presently  under  investigation.  As  already 
mentioned,  another  important  attribute  is  the  ability  to  generate  directly,  again 
by  extension  to  complex  (g'.t')  coordinates,  the  response  to  the  pulsed  beam  in 
any  environment  from  the  impulse  response  in  that  environment. 
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A.  X-RAY  COUPLED  WAVE  INTERACTIONS  AT  CRYSTAL 
SURFACES 

Professors  H.J.  Juretschke  andB.  Post 
Unit  SS7-1 


1.  OBJECTTVE(S) 

To  use  the  multiple  interaction  of  x-rays  in  a  Bragg  geometry  in  order  to 
develop  simple  and  compact  methods  for  obtaining  direct  phase  information 
about  the  crystal  scattering  factors,  and  for  characterizing  the  mode  structure; 
to  extend  the  interactions  to  include  coupling  to  other  waves  that  can  modify 
the  mode  structure,  and  that  can  be  used  to  explore  nonlinear  interactions  of 
x-rays;  and  to  understand  the  local  x-ray  fields  in  the  immediate  surface  region 
of  the  solid,  as  well  as  the  effect  of  stringent  boundary  conditions  on  all 
modified  x-ray  waves,  especially  those  originating  in  the  interior.  High  resolu¬ 
tion  x-ray  diffraction  experiments  will  be  supported  by  theoretical  studies  of  the 
predictions  of  n-beam  dynamical  theory,  and  by  rigorous  extensions  of  the 
theory  to  include  the  coupling  with  other  waves. 

2.  SUMMARY  OF  RECENT  PROGRESS 

The  progress  made  by  the  two  investigators  In  this  unit  is  summarized 
separately  below. 

The  main  research  on  the  dynamical  theory  of  coupled  wave  phenomena 
for  x-rays  (H.J.  Juretschke)  has  concentrated  on  probing  the  validity  of  the 
approximate  n-beam  formulation  developed  earlier  [l]  in  the  presence  of  devia¬ 
tions  from  the  perfect  periodicity  which  forms  the  usual  underlying  starting 
condition  for  most  dynamical  effects.  This  work  establishes  the  framework  for 
applying  the  theory  to  real  crystals,  and  at  finite  temperatures.  In  addition,  that 
theory  has  been  expressed  explicitly  for  applying  corrections  to  high  precision 
structure  factor  determinations  due  to  the  neighboring  n-beam  interaction 
regions.  High  precision  experimentation  has  now  reached  the  point  where  such 
corrections  cannot  be  ignored.  As  to  earlier  work,  the  two  papers  discussed  in 
last  year’s  report  have  not  been  published  (2,3). 


Two  papers  dealing  with  the  experimental  determination  of  x-ray  reflection 
phases  have  been  published  in  1987.  These  dealt  with  the  determination  of  the 
phases  of  ultra  weak  reflections  [4]  and  with  the  optimization  of  experimental 
conditions  for  enhancing  the  visibility  of  phase  indications  [5].  Two  additional 
papers  dealing  with  the  determination  of  the  phases  of  reflections  from  acentric 
organic  crystal  (benzil),  are  being  prepared  for  publication. 
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3.  STATE  OF  THE  ART  AND  PROGRESS  DETAILS 
A.  Theory  (H.J.  Juretschke) 

(l)  The  Effect  of  Thermal  Phonons  on  Diffraction  Near  an  n-Beam  Point 

Two  of  the  lines  of  investigation  we  have  been  developing  previously  have 
been  combined  in  this  study.  These  are:  the  theory  of  diffraction  near  an  n- 
beam  point;  and  the  effect  of  phonons  on  the  dynamical  diffraction  of  x-rays. 
The  main  result  of  the  first  line  of  research  has  been  that  the  first  order  correc¬ 
tions  due  to  the  presence  of  the  n-beam  point  can  be  taken  into  account  fully 
within  the  traditional  two-beam  formulation  merely  by  modifying  both  the 
structure  factors  and  the  origin  of  the  dispersion  surface  (Lorentz  point)  of  the 
pure  two-beam  case.  The  major  result  of  the  second  line  of  research  is  the 
incorporation  of  phonons  within  a  self-consistent  dynamical  approach,  such  that 
phonon-excited  x-ray  fields  are  also  fully  dynamical. 

The  influence  of  phonons  on  diffraction  near  an  n-beam  point  is  of  interest 
for  a  number  of  reasons.  Since  phonons  destroy  the  translational  periodicity  of 
a  perfect  crystal,  such  a  study  is  a  test  of  how  many  of  the  dynamical  features 
predicted  by  theory  for  a  perfect  crystal  survive  when  the  perfection  is  reduced, 
either  by  phonons  or  by  the  practically  equivalent  deviations  from  periodicity 
caused  by  weak  mosaic  structure.  In  addition,  such  a  study  can  make  quantita¬ 
tive  predictions  about  the  influence  of  temperature  on  diffraction  when  more 
than  one  structure  factor  is  involved.  It  is  common  practice  in  that  case  to 
assign  a  Debye-Waller  factor  correction  separately  to  each  structure  factor, 
although  the  validity  of  this  practice  has  apparently  never  been  investigated  in 
detail.  It  is  most  likely  based  on  the  argument  that  for  incoherent  thermal  pho¬ 
nons  all  multiple  scattering  effects  between  different  diffracting  planes  average 
out  to  zero.  Finally,  such  a  study  also  sheds  light  on  how  the  modified  two- 
beam  formulation  of  diffraction  near  an  n-beam  point  is  affected  by  tempera¬ 
ture,  and  in  what  manner  it  may  need  to  be  renormalized  at  each  temperature. 

It  has  now  been  shown  In  our  study  that  phonons  can  be  incorporated  in 
the  modified  two-beam  formulation  without  any  complications,  and  that  the 
main  results  of  their  presence  are: 

a)  The  intrinsic  modified  two-beam  parameters  of  diffraction  'structure  fac¬ 
tors,  absorption,  etc.)  remain  completely  the  same  as  in  the  absence  of 
phonons.  Just  as  in  a  pure  two-beam  case.  The  only  change  is  that,  Just  as 
in  that  case,  there  appear  Debye-Waller  factors  and  thermal  diffuse  scatter¬ 
ing. 

b)  The  Debye- Waller  factor  is  entirely  that  of  the  primary  two- be  am  interac¬ 
tion,  and  is  a  prefactor  of  the  modified  structure  factor.  The  other  two- 
beam  structure  factors  contributing  to  this  modified  structure  factor  show 
no  additional  temperature  softening.  Hence  none  of  the  fine  structure  in 
diffraction  associated  with  the  modified  structure  factor  disappears  in  the 
presence  of  phonons,  but  is  merely  broadened  in  the  manner  characteristic 
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of  all  temperature  effects  In  two-beam  diffraction. 

c)  Since  the  above  is  true  for  deviations  from  periodicity  caused  by  weak 
mosaic  structure  (thermal  phonons  are  practically  static  on  the  time  scale 
of  x-rays).  Hence,  certainly  mildly  imperfect  crystals  are  expected  also  to 
show  all  the  diffraction  fine  structure  of  perfect  crystals,  suitably  broadened 
by  characteristic  averaging.  Of  course,  since  this  averaging  may  involve  a 
distribution  of  asymmetric  line  features,  some  of  this  asymmetry  may 
become  obscured. 

The  underlying  reason  for  the  above  results  is  that  phonon-excited  x-rays 
near  an  n-beam  situation  are  also  subject  to  n-beam  interactions,  and  that  the 
many  more  coherent  coupling  channels  existing  here  undergo  cancellations 
leading  to  a  net  effect  due  to  the  primary  interaction  alone. 

It  is  concluded  from  this  study  that  both  thermal  and  mosaic  effects  can  be 
handled  in  the  neighborhood  of  an  n-beam  point  entirely  in  the  traditional 
two-beam  fashion,  but  now  with  respect  to  the  modified  two-beam  parameters. 
Therefore,  the  dynamical  theory  of  these  effects  is  not  confined  to  perfect  cry¬ 
stals,  as  already  established  experimentally  (e.g.  see  Ref.  6).  This  work  has 
been  accepted  for  publication  in  Acta  Crystallographies-  A  preliminary  report 
of  it  was  presented  at  the  14th  International  Congress  of  Crystallography  in 
Perth,  Australia,  in  August  1S87. 

(2)  Corrections  Due  to  Multple  Beam  Interactions  in  High  Accuracy 
Two-Beam  Structure  Factor  Determinations 

Recently  there  has  been  increased  activity  in  the  determination  of  two- 
beam  structure  factors  to  high  accuracy,  as  a  means  of  studying  electron  distri¬ 
butions  in  chemical  bonds  of  valence  crystals  such  as  Si.  At  the  0.1%  level  of 
accuracy,  it  becomes  necessary  to  take  into  account  a  number  of  corrections, 
among  which  is  the  contribution  from  neighboring  multiple  interaction  points. 

Our  previously  developed  theory  for  the  asymptotic  effects  of  such  interac¬ 
tions  gives  a  full  prescription  for  the  needed  corrections,  which  is  quantitati\ », 
and  which  can  be  easily  carried  out  once  the  relative  location  and  identification 
of  the  neighboring  Interaction  points  are  known.  Hence  it  was  considered 
worthwhile  to  formulate  such  a  correction  procedure  as  an  explicit  application 
of  the  theory.  As  applied  to  recent  measurements  of  Sl(l2  12  0)  with  AgKa 
radiation  (6],  it  predicts,  under  the  most  favorable  circumstances,  corrections  at 
least  of  the  order  of  0.1%.  However,  since  in  the  tradition  of  most  such  meas¬ 
urements,  no  information  is  supplied  about  the  exact  orientation  of  any  two- 
beam  experiment  in  reciprocal  space,  it  is  difficult  to  apply  the  method  to  a 
specific  situation,  in  which  it  is  likely  that  the  corrections  would  in  fact  be 
larger.  But  the  sample  computation  makes  clear  that  such  corrections  cannot  be 
ignored  at  the  0.1%  level;  In  order  to  be  applied,  of  course,  the  experimental 
information  must  include  a  Tull  indexing  of  all  neighboring  multiple  interaction 
points.  It  is  hoped  that  as  a  result  of  this  study,  such  information  will  be  pro¬ 
vided  routinely  in  future  experimental  work. 
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This  work  was  presented  at  the  "Accuracy  in  Structure  Factor  Determina¬ 
tion”  Symposium  in  Warburton,  Victoria,  Australia,  in  August  1987.  It  will  be 
published  in  the  Australian  Journal  of  Physics. 

B.  Experimental  (B.  Post) 

During  the  past  year  we  have  continued  our  investigation  of  relevant  vari¬ 
ables  involved  in  the  experimental  determination  of  x-ray  reflections  phases. 
The  sensitivity  of  our  procedures  has  been  demonstrated  by  determining  the 
phases  of  the  very  weak  "forbidden"  reflections  of  germanium,  using  CuKal 
radiation.  Invariant  triplets  of  germanium  reflections  were  generated  in  n-beam 
patterns.  Each  triplet  Included  two  strong  reflections  whose  phases  were 
known,  in  addition  to  the  one  "forbidden"  reflection.  Triplet  phases  were 
clearly  displayed  in  the  asymmetric  line  profiles  of  the  n-beam  interactions;  sub¬ 
traction  of  the  pairs  of  known  phases  yielded  the  phases  of  the  corresponding 
forbidden  terms  [4]. 

Experimental  conditions  for  optimizing  the  phase  Indications  shown  in  n- 
beam  patterns  have  also  been  studied. 

Large  numbers  of  phases  of  acentric  organic  crystal  (benzil)  have  been 
determined  using  CuKal  and  CrKal  radiations.  The  experimental  findings 
correspond  satisfactorily  with  phases  calculated  from  published  data.  They  indi¬ 
cate  that  experimental  values  of  phase  angles  of  individual  reflections  can  be 
determined  to  ±  45°,  or  and  possibly  to  ±  30°.  Providing  that  information  to 
crystallographers  should  make  possible  the  determination  of  many  crystal  struc¬ 
tures  which  have  resisted  solution  up  to  the  present. 
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B.  MICROSTRUCTURE  PHOTOPHYSICS 


Professor  S.  Arnold  and  K.M.  Leung 


Unit  SS7-2 


1.  OBJECITVE(S) 

Our  objective  is  to  develop  a  quantitative  understanding  of  the  detailed 
interaction  of  electromagnetic  radiation  with  microstructures  of  dimensions  (a) 
comparable  to  and  much  less  than  the  wavelength  (X)  of  the  incident  radiation. 
In  some  cases  the  broad  features  of  this  Interaction,  and  especially  the  enor¬ 
mous  resonant  enhancements,  may  be  understood  qualitatively  by  using  a  clas¬ 
sical  model  for  the  local  electromagnetic  field  (e.g.,  surface  enhanced  Raman 
scattering)  [l].  In  others  this  approach  has  not  proven  to  be  fruitful  (e.g., 
enhanced  photoemission  from  small  Ag  particles)  [2].  The  lack  of  a  cohesive 
framework  must  in  part  be  due  to  the  limitations  of  the  current  classical  model 
in  the  face  of  the  properties  of  real  materials  and  structures.  Unfortunately 
careful  attention  has  not  been  paid  to  these  limitations.  This  neglect  has  usually 
been  dictated  by  experimental  bounds  or  by  theoretical  limitations.  But,  at  this 
time,  a  fuller  understanding  of  the  limitations,  and,  by  implication,  a  more  sys¬ 
tematic  exploitation  of  their  remarkable  consequences/ rests  precisely  on  more 
attention  to  detail  and  on  quantitative  comparison  of  experiment  and  theory. 

Here  we  propose  such  a  study,  relying  on  tight  interplay  between  theory 
and  experiment,  that  emphasizes  both  far  field  and  near  field  features  such  as 
elastic  scattering  and  fluorescence,  on  single  structures  of  fully  characterized 
geometry,  and  extending  into  a  limit  a/X<<i  where  particle  properties  are 
expected  to  deviate  from  their  macroscopic  behavior. 

In  the  coming  period  special  emphasis  will  be  placed  on  nonlinear  processes. 
This  is  an  especially  challenging  area  since  nonlinear  theory  for  the  interaction 
of  radiation  with  a  particle  of  a  generalized  size  is  virtually  nonexistent.  How¬ 
ever,  in  the  particular  case  of  spherical  particles,  and  most  probably  for  particles 
of  other  regular  morphologies,  resonances  having  unprecedentedly  high  Q's 
have  been  detected.  Thus  the  field  strength  within  a  particle  even  in  the  pres¬ 
ence  of  a  laser  of  modest  power  easily  produces  a  sizable  nonlinear  effect.  For 
example.  Just  this  last  period  we  have  demonstrated  Optical  Bistability  with  a 
threshold  intensity  in  the  W/cm2  range  [3],  More  will  be  said  about  this  effect 
in  the  summary.  All  of  our  work  thus  far  implies  that  this  area  is  fertile  for 
both  basic  and  practical  results.  On  the  practical  side  we  see  a  new  technology 
evolving.  Microparticle  Photonics. 
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2.  SUMMARY  OF  RECENT  PROGRESS 


This  section  presents  a  brief  summary  of  recent  progress;  more  detailed 
descriptions  are  contained  in  the  next  section  la  conjunction  with  the  state  of 
the  art  so  that  the  nature  of  the  contributions  can  be  understood  more  clearly. 

The  program  in  Microstructure  Photophysics  is  composed  of  two  interac¬ 
tive  parts.  Experiments  are  carried  out  at  the  Microparticle  Photophysics 
Laboratory  (MP3L)  under  the  direction  of  S.  Arnold,  while  theoretical  work  on 
associated  phenomena  is  carried  out  by  KM.  Leung.  We  have  been  fortunate 
during  the  past  year  in  as  much  as  the  group  has  grown  by  one  member  (L.M. 
Folan)  and  the  laboratory  has  been  enlarged  to  Include  a  dynamics  component 
with  measurement  capabilities  down  to  a  few  picoseconds.  This  past  year  has 
been  particularly  productive  in  two  major  areas  which  will  be  discussed  below 
(Secs.  2A  and  2B).  In  addition  the  past  year  has  brought  further  Independent 
funding  for  the  experimental  and  theoretical  studies  from  NSF  (to  S.  Arnold 
and  LM.  Folan)  and  ONR  (to  KM.  Leung),  respectively.  . 

A.  Rayleigh  Particle  Optical  Bistability  (OB) [4] 

The  primary  report  predicting  Rayleigh  particle  OB  in  semiconductors  is 
now  published.  Over  this  last  year  these  calculations  have  been  extended  to  a 
nonsphericai  particle  where  Optical  Tristability  is  found.  Although  our  work  is 
seminal  In  this  area,  a  number  of  workers  are  currently  Jumping  into  the  fray. 
A  review  of  the  past  year’s  theoretical  progress  in  different  labs  will  be  dis¬ 
cussed  In  Section  3. 

B.  Mie  Particle  Optical  Bistability  [3] 

In  A.  the  work  is  theoretical.  The  Insights  gained  from  these  models  have 
been  used  in  the  MP®L  to  provide  the  first  demonstration  of  microparticle  OB. 
This  work  has  recently  been  published,  and  a  patent  disclosure  is  in  progress. 
Work  over  the  last  year  has  taken  our  Initial  single  resonance  model  for  the 
optical  bistability  of  a  cloud  droplet  and  extended  It  to  multiple  resonances 
through  a  model  which  can  be  considered  nearly  comprehensive. 

A  number  of  other  research  projects  which  are  relevant  to  the  above  have 
been  investigated  during  the  year  and  have  led  to  additional  publications.  In  ail 
the  publications  already  credited  to  JSEP  for  the  current  contract  are: 

1.  K.M.  Leung,  "Propagation  of  Nonlinear  Surface  Polaritons,"  Phys.  Rev.A, 

31r  1180  ( 1985) . 

2.  S.  Arnold,  E.K.  Murphy,  and  G.  Sageev,  "Aerosol  Particle  Molecular  Spec¬ 
troscopy,"  Appl.  Opt.,  24,  1048  (1985). 
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3.  K.M.  Leung,  "Aerosols  of  Anisotropic  Metallic  Microparticles  as  Artificial 
Kerr  Media,"  Opt.  Lett.,  IQ,  347  (1985). 

4.  A.B.  Plucliino  and  S.  Arnold,  "A  Comprehensive  Model  of  the  Photo- 
phoretic  Force  on  a  Spherical  Microparticle,"  Opt.  Lett.,  261  (1985). 

5.  L.M.  Folan,  S.  Arnold  and  S.D.  Druger,  "Enhanced  Energy  Transfer  within 
a  Microparticle,”  Chem.  Phys.  Lett.,  118.  322  (1985). 

6.  S.  Arnold  and  N.  Hessel,  "Photoemission  from  Single  Electrodynamically 
Levitated  Microparticles,"  Rev.  Sci.  Inst.  1086  (1985). 

7.  K.M.  Leung,  "P-Polarized  Nonlinear  Surface  Polaritons  in  Materials  with 
Intensity  Dependent  Dielectric  Functions,"  Phys.  Rev.  B,  22L*  5093  (1985). 

8.  K.M.  Leung,  "Optical  Bistability  in  the  Scattering  and  Absorption  of  Light 
from  Nonlinear  Microparticles,"  Phys.  Rev.  A  2£,  2461  (1986). 

9.  SA.rnold  and  L.M.  Folan,  "  A  Fluorescence  Spectrometer  for  a  Single 
Electrodynamically  Levitated  Microparticle,"  Rev.  Sci.  Inst.,  ££,  2250 
(1986). 

10.  S.  Arnold,  K.M.  Leung  and  A.B.  Pluchino,  "The  Optical  Bistability  of  an 
Aerosol  Particle,"  Opt.  Lett.  U,  800  (1988). 

11.  K.M.  Leung  and  He  Ming-Guo,  "Tristabie  Optical  Behavior  of  Anisotropic 
Rayleigh-sized  Microparticles,"  Submitted  to  J.  Wave-Material  Interaction. 

12.  S.D.  Druger,  S.  Arnold,  and  L.M.  Folan,  "Theory  of  Enhanced  Energy 
Transfer  between  Molecules  Embedded  in  Spherical  Dielectric  Microparti¬ 
cles,"  J.  Chem.  Phys.  27,  2849  (1987). 

13.  S.  Arnold,  and  L.M.  Folan,  "A  Spherical  Void  Electrodynamic  Levitator,” 
Rev.  Sci.  Inst.  2J£,  1732  (1987). 

14.  L.M.  Folan  and  S.  Arnold,  "Determination  of  Molecular  Orientation  at  the 
Surface  of  an  Aerosol  Particle  by  Morphology  Dependent  Photoselection," 
Opt.  Lett.  (Jan.-88). 

15.  S.  Arnold,  "Spectroscopy  of  Single  Levitated  Micron  Sized  Particles,"  in 
Optical  Effects  Associated  with  Small  Particles  Eds.  P.W.  Barber  and  R.K. 
Chang  (World  Scientific  Publishing,  1988). 
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3.  STATE  OF  THE  ART  AND  PROGRESS  DETAILS 

The  area  of  nonlinear  physics  is  concerned  with  physical  phenomena  that 
are  based  on  a  nonlinear  response  of  a  medium  to  incident  electromagnetic 
radiation.  The  nonlinearities  responsible  for  most  of  these  phenomena  in  the 
past  are  often  quite  weak  and  therefore  the  processes  involved  are  actually 
quasi-linear.  Consequently,  theoretical  threatments  have  been  almost 
exclusively  perturbative  in  nature,  and  experimentally,  no  minimum  threshold 
powers  are  needed  to  see  these  effects. 

Recently  there  are  growing  interests  in  intrinsically  (or  fully)  nonlinear 
interactions  of  electromagnetic  radiation  with  matter.  The  nonlinearities  appear¬ 
ing  in  the  dielectric  function  are  rather  strong  and  Maxwell’s  equations  become 
a  set  of  coupled  nonlinear  partial  differential  equations.  Calculations  of  these 
effects  must  be  carried  out  in  a  fully  non-pcrturbativc  fashion.  The  arsenal  of 
mathematical  techniques  such  as  Green’s  functions,  Fourier  transforms  and 
integral  equation  methods,  which  are  so  successful  in  dealing  with  linear  prob¬ 
lems,  are  totally  incapable  of  treating  these  intrinsically  nonlinear  physical 
phenomena.  The  associated  effects  are  not  merely  small  extensions  of  the  linear 
or  quasi-linear  ones,  but  must  be  regarded  as  separate  entities  with  characteris¬ 
tic  behavior  such  as  hysteresis  and  switching,  having  no  counterpart  what-so- 
ever  in  the  linear  limit.  With  the  incident  field  as  an  extra  physically  controll¬ 
able  parameter,  which  actively  and  quite  spontaneously  changes  the  effective 
value  of  the  dielectric  function  of  the  medium,  a  host  of  entirely  new 
phenomena  are  possible. 

The  distinct  nature  of  nonlinear  effects  is  especially  true  in  the  case  of 
interactions  of  intense  electromagnetic  radiation  with  microstructures  where 
such  effects  are  expected  to  be  enhanced  owing  to  two  separate  mechanisms. 
The  first  is  the  so-called  dielectric  confinement  or  local  field  effect  which 
increases  the  local  field  intensity  in  the  neighborhood  of  microstructures  when 
the  frequency  is  close  to  the  surface  plasmon  or  morphological  resonances.  This 
mechanism  is  responsible  for  the  sizable  enhancement  of  a  variety  of  linear 
processes,  such  as  surface  enhanced  Raman  scattering  [1],  as  well  as  nonlinear 
processes,  such  as  the  lowering  of  thresholds  for  lasing  [5]  and  optical  switching 
[4].  It  is  advantageous  to  have  high  local  Intensity  especially  for  nonlinear 
processes  since  the  magnitudes  of  the  effects  go  up  faster  than  linearly  with 
intensity.  The  second  mechanism  has  to  do  with  the  fact  that  noniinearities  in 
microcrystals  and  microparticles  can  actually  be  much  larger  than  those  found 
in  bulk  materials.  This  is  due  to  electronic  localization  or  quantum  size  effects, 
which  are  responsible  for  a  host  of  unusual  phenomena,  including  enhanced 
nonlinearities,  even  at  room  temperature,  in  semiconductor  heterojunctions  and 
superlattice  structures  [6].  The  electrons  and  holes  are  squeezed  into  a  two 
dimensional  space.  In  the  case  of  microparticles  electrons  are  confined  on  all 
sides  by  the  boundaries  and  therefore  they  are  expected  to  behave  more  like 
zero-dimensional  objects  or  quantum  dots,  with  enhanced  exciton  lifetime. 

Our  research  in  this  exciting  new  area  has  resulted  In  some  very  encourag¬ 
ing  results  which  will  be  discussed  below.  Our  work  has  been  followed  up  by 
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others,  especially  the  group  at  Bell  Laboratories  [7]. 

A.  Rayleigh  Particle  Optical  Bistability  (OB)  and  Optical  Tristability  (OT) 

For  a  given  optically  nonlinear  material  the  OB  behavior  of  a  spherical 
Rayleigh  particle  can  only  be  controlled  with  the  intensity  and  the  frequency  of 
the  incident  light.  In  order  to  increase  the  capability  of  microparticle  photonics, 
we  want  to  investigate  more  complex  microparticle  devices.  For  this  purpose  we 
have  recently  extended  our  work  on  the  OB  of  spherical  Rayleigh  particles  to 
nonspherical  (spheroidal)  particles  (8j. 

We  And  that  because  of  their  additional  physical  characteristics,  such  as  the 
particle  shape  and  the  orientation  of  the  particle  with  respect  to  the  incident 
light,  spheroidal  particles  are  capable  of  exhibiting  a  much  richer  class  of  intrin¬ 
sically  nonlinear  behavior.  Depending  on  the  frequency  of  the  incident  light  and 
the  particle  shape  and  orientation,  we  And  that,  besides  OB  behavior,  a  total  of 
Ave  different  types  of  optical  tristability  (OT)  behaviors  are  possible.  The 
regions  in  parameter  space  where  these  Ave  (Iv)  types  of  OT  behaviors  occur 
are  shown  in  Fig.  1  for  Axed  particle  shape.  Region  (vi)  where  only  OB  can 
exist,  and  region  (vti)  where  the  internal  intensity  is  a  single-valued  function  of 
the  incident  intensity  are  also  shown. 

We  should  also  point  out  that  nonspherical  particles  do  exist,  and  can  be 
found  in  the  earth’s  atmospheres,  interplanetary  space,  and  In  interstellar 
media.  In  addition  such  particles  can  be  fabricated  in  the  laboratory.  In  fact  a 
two  dimensional  array  of  silver  particles  each  shaped  as  an  ellipsoid  of  revolu¬ 
tion  has  been  made  a  number  of  years  ago  for  the  study  of  surface  enhanced 
Raman  scattering  [9],  In  addition  Ashkin  has  demonstrated  that  such  particles 
can  be  constructed  through  optical  levitation  [lOj. 

B.  Mie  Particle  Optical  Bistability  [3] 

Unlike  the  Rayleigh  particle  we  require  no  frequency  dependence  in  the 
dielectric  function  to  insure  Aeld  enhancement  in  the  Mie  case.  Resonances  are 
assured  by  the  phase  matching  of  waves  which  circumnavigate  the  perimeter  of 
the  particle.  Indeed  these  resonances  may  have  Q’s  greater  than  10s  [lij.  Con¬ 
sequently,  the  Aeld  enhancements  within  the  particles  are  projected  to  be  enor¬ 
mous.  One  recent  article  calculates  the  Aeld  at  the  surface  of  a  particle  several 
microns  in  size  to  be  9  orders  of  magnitude  larger  than  the  Aeld  of  the  wave 
Incident  on  the  particle(  under  resonant  conditions)  (12).  Thus  a  slight  change 
in  the  dielectric  function,  size  or  shape  of  the  particle  with  intensity  can  lead  to 
catastrophic  effects.  This  is  best  demonstrated  by  our  work  on  a  simulated 
cloud  droplet. 

A  cloud  droplet’s  size  is  controlled  by  transport  through  its  phase  boun¬ 
dary.  If  the  droplet  is  placed  in  an  atmosphere  at  the  same  vapor  pressure  as  its 
surface,  it  does  not  evaporate.  Optical  absorption  by  the  droplet  causes  its 
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Fig.  1  Diagram  showing  the  parameter  regions  which  exhibit  different 
types  of  nonlinear  optical  response  of  a  spheroidal  particle  with  a 
self-focusing  Kerr  nonlinearity.  Here  q  Is  the  angle  between  the 
polarization  of  the  Incident  light  and  the  long  symmetry  axis  of  the 
particle,  and  d  Is  the  frequency  detuning. 

surface  vpnor  pressure  to  increase,  and  Its  size  to  decrease.  If  the  absorbed 
energy  is  from  a  narrow  band  laser,  one  observes  bistable  behavior  in  light 
scattering  from  the  droplet  by  scanning  wavelength  [3]  or  intensity  [13].  In  Ref. 
[3]  we  developed  a  heurestic  model  for  this  behavior  based  on  an  optical 
response  having  a  simple  Lorentzlan  shape.  Over  this  last  period,  we  have 
taken  a  different  approach  in  order  to  provide  a  more  comprehensive  model. 
Our  approach  is  a  graphical  one  which  tries  to  predict  "equilibrium"  states 
through  the  graphical  solution  of  two  equations.  The  first  is  an  equation  which 
expresses  the  fact  that  the  optical  size  X  (circumference/wavelength)  decreases 
in  the  presence  of  light  of  intensity  I  [14,15], 


X  =  X0  -  rQa(X)I  . 


(1) 


The  second  is  an  expression  for  the  absorption  efficiency  Qa(X)  as  obtained 
from  Mie  theory. 


co 


Qa(x)  S  ( 2q+! )  j  fq(X.mX)  +  gq(X,mX) 


Xa 


(2) 
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where  m  is  the  refractive  index  and,  fq  and  gq  are  the  absorption  contributions 
associated  with  internal  fields  with  angular  mode  number  q  and  are  of  TE  and 
TM  type  respectively.  Since  the  cloud  droplet  interaction  occurs  at  an  intensity 
less  than  100W/cm2  [3],  refractive  index  variations  are  not  anticipated  and  the 
right  hand  side  of  Eq.  (2)  is  simply  a  function  of  X.  Unfortunately,  even  with 
this  simplification  an  analytical  solution  of  Eqs.  (1)  and  (2)  is  not  possible, 
however  a  graphical  approach  provides  an  accurate  means  for  obtaining  results 
from  which  physical  implications  can  be  drawn.  Figure  2  shows  the  basic 
scheme.  The  straight  line  is  a  plot  of  Eq.  (1)  while  the  curved  line  represents  a 
calculation  based  on  Eq.  (2)  over  a  small  region  in  optical  size.  The  meaning  of 
X10  Is  that  the  calculation  is  carried  out  at  a  particular  wavelength  Xj  for  which 
the  optical  size  at  an  infinitesimally  small  Intensity  is  X10  (and  the  correspond¬ 
ing  radius  Is  a^).  Since  intensity  controls  the  slope  of  the  straight  line,  absorp¬ 
tion  by  the  particle  is  particularly  sensitive  to  intensity.  By  the  same  token  an 
increase  in  wavelength  for  the  same  intensity  is  represented  by  shifting  the 
straight  line  to  the  right  without  changing  its  slope.  Such  a  wavelength  scan  also 
produces  a  significant  effect  on  the  particle  absorption  efficiency.  Both  of  these 
effects  are  most  pronounced  in  the  region  of  a  morphological  resonance.  In 
such  a  region  either  fq  or  gq  is  resonantly  enhanced  and  bistable  effects  are  antici¬ 
pated. 


Q  * 
a 


tan (tl) 


Xjo  X|(I) 

Fig.  2  Graphical  scheme  for  simultaneous  solution  of  Eqs.  (l)  and  (2). 

Ftgure  3  shows  the  results  of  an  actual  calculation  on  a  particle  of  refrac¬ 
tive  Index  1.38+10"  5l  and  for  which  the  unperturbed  radius  is  4.828/i.  The 
apparent  Lorentzlan  response  In  Fig.  3a  Is  due  to  the  resonant  stimulation  of 
the  gjjQ  coefficient  in  Eq.  (2)  (a  TM  mode).  The  simulation  was  carried  out  for 
an  intensity  of  50W/cm2  for  which  the  particle  size  vs.  absorption  is  depicted  by 
the  straight  dashed  lines.  As  one  changes  the  wavelength  of  the  laser  by  a  few 
angstroms  from  2jt Sq/x^  to  27ra0l/x20  we  see  that  the  number  of  solutions  for 
absorption  changes  from  3  to  l.  The  multivaluedness  of  the  absorption  is  most 
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apparent  when  the  solutions  are  plotted  against  wavelength  (Fig.  3b).  The  tran¬ 
sitions  shown  in  Fig.  3b  correspond  to  the  points  of  infinite  slope.  At  these 
points  catastrophes  will  take  place.  The  abrupt  changes  seen  in  light  scattering 
[3]  may  be  calculated  by  mapping  the  various  solutions  at  a  given  wavelength 
onto  the  scattering  response  (Fig.  3c).  The  overall  light  scattering  response  vs. 
wavelength  is  shown  in  Fig.  3d.  With  this  scheme  we  have  predicted  a  host  of 
new  phenomena  including  the  observation  of  hidden  resonances  [16]  and 
multistability  [17] .  In  some  sense  one  may  consider  the  scheme  depicted  in  Fig. 
3  to  be  exact.  However,  It  should  be  emphasized  that  no  dynamics  are  included 
in  this  diagram  and  the  particular  problem  is  perhaps  the  simplest  Mie  bistabil¬ 
ity  problems  which  one  can  Imagine. 


0.0033T 


0.0022 L 


b  a 


TV  7 

n  / 

I  \j  / 

n  / 


d  c 


588  591 

WAVELENGTH  (nm) 


2583  25.73  25.63 

OPTICAL  SIZE,  X 


Fig.  3  The  results  of  calculation  on  a  cloud  droplet  particle  of  refractive 
index  1.38+10" 5I  and  for  which  the  unperturbed  radius  is  4.828/i. 

The  next  step  is  clearly  to  take  a  careful  look  at  the  Mie  problem  in  the  pres¬ 
ence  of  a  Kerr  effect.  Both  experimental  and  theoretical  work  on  this  problem 
is  currently  underway. 
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1.  OBJECTIVE(S) 

To  develop  a  theoretical  and  experimental  understanding  of  some  specific 
Interactions  of  surfaces  with  electrons,  phonons,  atoms,  and  with  each  other, 
that  govern  the  dynamical  properties  of  surfaces  and  interfaces  -  the  emerging 
crucial  "missing  link"  between  current  surface  science  studies  and  the  practical 
behavior  of  surfaces.  Since  no  unified  framework  exists  to  attach  the  myriad  of 
problems  of  probable  importance,  we  will  study  three  situations,  representative 
of  the  range  of  questions  that  must  be  answered,  that  offer  the  greatest  hope 
for  fruitful  analysis  and  potential  impact  In  this  new  area.  The  general  theme 
of  phonon  interactions,  from  microscopic  to  macroscopic  time  and  length 
scales,  couples  the  three  target  problem  areas  within  the  field  of  dynamical  and 
nonequilibrium  phenomena. 

Under  this  common  theme,  the  specific  theoretical  problems  chosen 
address  the  relaxation  mechanisms  of  hot  electrons  in  very  small  structures  and 
the  dynamics  of  surface  reactions  with  incoming  atoms  or  molecules.  Accom¬ 
panying  experiments  will  concentrate  on  developing  second  harmonic  light  gen¬ 
eration  In  Interface  regions  as  a  tool  for  studying  time-dependent  strains  in 
those  regions. 

The  emphasis  of  most  current  surface  studies  on  the  gross  aspects  of  sur¬ 
faces,  such  as  their  electronic  properties  or  their  average  structural  relaxation, 
has  been  very  successful  In  understanding  many  of  the  major  new  features  seen 
or  expected  there.  However,  at  the  next  stage  of  detailed  quantitative  descrip¬ 
tions  of  particular  phenomena,  evidences  of  new  complexities  are  accumulating 
that  point  to  the  need  for  re-examining  critically  the  underlying  theory  in  order 
to  understand  more  fully  the  restrictions  and  new  conditions,  as  well  as  the  new 
degrees  of  freedom,  imposed  by  the  existence  of  a  surface,  and  by  its  dynami¬ 
cal  response.  For  example,  electrons  participating  in  current  transport  in  a  nar¬ 
row  surface  region  may  encounter  closely  spaced  barriers  before  being  able  to 
shed  the  excess  energy  gained  In  local  regions  of  high  field,  and  thus  continue 
to  remain  "hot."  Similarly,  atoms  colliding  with  a  surface  often  do  so  in  times 
so  short  that  they  do  not  allow  for  a  full  exchange  of  energy  to  reach  thermaii- 
zatlon.  The  common  features  of  phenomena  such  as  these  are,  first  of  all,  that 
they  deal  with  processes  that  start  out  far  from,  and  may  never  reach,  steady 
states  close  to  thermal  equilibrium;  and,  secondly,  that  because  of  the  existence 
of  the  surface,  or  an  interface,  the  Importance  of  various  relaxation  mechan¬ 
isms  may  differ  markedly  from  their  relative  strengths  in  bulk  material. 
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2.  SUMMARY  OF  RECENT  PROGRESS 

A.  Electron-Phonon  Scattering  by  Hot  Electrons  in  Quantum  Well 
Microstructu  *es 

The  research  has  focused  on  the  mechanisms  whereby  optically  excited 
electrons  in  the  immediate  vicinity  of  an  interface  relax  towards  lower  energy 
states.  We  have  developed  a  formalism  that  is  capable  of  describing  the  scatter¬ 
ing  of  electrons  by  interface  optical  phonons.  Unlike  previous  formulations 
[14,15]  the  dynamical  properties  of  the  materials  on  both  sides  of  the  interface 
play  important  roles  in  the  scattering  processes.  The  equilibration  of  optically 
excited  electrons  proceeds  mainly  by  intersubband  scattering,  since  this  dom¬ 
inates  the  rate  of  energy  loss. 

The  rates  of  decay  of  electronic  states  have  been  calculated  in  the  Born 
approximation.  However,  since  the  corresponding  wave  function  renormaliza¬ 
tions  are  unusually  large,  these  results  are  not  conclusive.  The  inclusion  of 
higher  order  processes  results  in  a  significant  change  of  the  decay  rates.  There¬ 
fore,  a  strong  electron-phonon  coupling  theory  has  been  constructed  to  describe 
the  electron  (polaron)  state  formed  at  the  interface. 


3.  STATE  OF  THE  ART  AND  RESEARCH  DETAILS 

A.  Electron-Phonon  Scattering  by  Hot  Electrons  in  Quantum  Well 
Microstructures 

The  size  of  electronic  components  has  decreased  rapidly  over  the  last 
decade,  to  such  an  extent  that  many  of  the  classical  concepts  developed  for 
bulk  materials  have  become  completely  Inappropriate  for  their  description.  An 
important  consequence  of  the  reduction  of  the  size  of  electronic  devices  is  that 
the  internal  electric  fields  become  very  large.  Fields  of  the  order  of  l03V/cm 
are  quite  common  in  present  day  integrated  circuits,  while  field  strengths  of  the 
order  of  lOsV/cm  can  be  found  between  the  neighboring  gates  of  charge  cou¬ 
pled  devices  and  in  quantum  well  layered  devices.  At  such  large  fields,  the 
electrons  are  rapidly  accelerated  and  acquire  large  kinetic  energies.  The  net 
result  is  that  the  conducting  electrons  cannot  be  regarded  as  being  in  equili¬ 
brium  with  the  rest  of  the  microstructure,  and  the  conduction  process  becomes 
extremely  non-ohmic  at  room  temperatures  [2].  This  regime  is  characterized 
by  extremely  nonequilibrium  thermodynamics. 

The  manner  in  which  the  electron  gas  exchanges  energy  with  the  rest  of 
the  system  becomes  extremely  important.  As  we  shall  outline  later,  there  is 
considerable  evidence  that  the  dissipative  process  for  these  high  energy  elec¬ 
trons  is  completely  different  from  the  processes  that  usually  occur  for  the  elec¬ 
trons  in  thermal  equilibrium.  This  is  even  more  true  for  quantum  well  systems 
such  as  MOSFETs  [3]  or  quantum  well  heterojunctions  [4]  in  which  the 
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saturation  velocities  of  10s  cm/sec  are  considerably  lower  than  the  bulk  values 
of  107  cm/sec  [5].  The  presence  of  the  interfaces  may  be  of  crucial  importance 
in  the  description  of  the  nonequilibrium  effects. 

At  room  temperature,  it  is  generally  believed  that  the  mobility  is  governed 
and  limited  by  phonon  scattering,  and  that  most  of  the  excess  energies  are 
transferred  to  the  lattice  via  optical  phonons.  But,  as  Ando  et  al.  [l]  have 
pointed  out  in  their  recent  review  article,  the  theory  of  phonon  scattering  for 
the  electronic  mobility  is  quite  unsatisfactory.  The  calculated  mobility  [6]  is 
much  larger  than  the  experimental  value.  Even  more  distressing  is  that  the 
one-phonon  process  leads  to  a  mobility  that  is  almost  independent  of  the  den¬ 
sity  of  electrons  in  the  Inversion  layer.  The  experiments  show  a  Ns1/3  depen¬ 
dence  where  N,  is  the  electron  density. 

At  this  time,  therefore,  a  more  general  investigation  of  the  processes 
through  which  an  electron  can  emit  or  absorb  phonons  and  equilibrate  with  the 
lattice  is  called  for.  In  particular,  it  should  take  into  account  the  specific  nature 
of  the  phonon  spectra  near  the  interfaces,  as  well  as  the  possibility  that  the 
energy  exchange  between  the  electrons  and  phonons  may  become  a  resonant 
process.  This  possibility  is  offered  since  the  separations  of  the  various  sub¬ 
bands  nearly  match  the  bulk  phonon  energies  in  some  systems  [7] .  In  addition, 
the  role  of  the  effects  of  higher  order  phonon  processes  must  be  clarified.  In 
bulk  materials,  such  effects  are  usually  negligible,  as  may  be  implied  by 
Migdal’s  theorem  [8],  However,  in  quantum  well  systems  the  theorem  no 
longer  applies,  so  the  multiphonon  processes  could  be  of  importance.  The 
presence  of  these  multiphonon  events  may  strongly  influence  the  hot  electron 
relaxation  rates,  but  may  also  be  responsible  for  the  discrepancy  between  the 
experimental  data  on  mobilities  In  the  ohmic  regime  and  theoretical  calculations 
that  only  consider  single- phonon  scattering  events  [ l] .  Other  authors  [l]  have 
speculated  that  these  discrepancies  may  be  due  to  the  effects  of  electron  corre¬ 
lations  in  the  final  states,  but  apparently  no  one  has  investigated  the  possibility' 
of  the  multiphonon  processes  that  we  suggest. 

The  optical  properties  of  inversion  layers  show  related  anomalies,  presum¬ 
ably  also  associated  with  the  electron  phonon  scattering  mechanisms. 
Kneschaurek  and  Koch  [9]  have  studied  the  temperature  dependence  of  the 
intersuoband  optical  transitions.  These  and  related  experiments  (10,11]  show 
that  the  widths  of  the  resonance  lines  present  a  problem  at  high  temperatures 
In  that  the  widths  do  not  increase  concomitantly  with  the  increase  in  the 
phonon- induced  relaxation  rate  inferred  from  transport  properties.  This  raises 
the  question  of  how  well  the  transport  scattering  rate  is  correlated  with  the 
corresponding  optical  rate  in  the  presence  of  multiphonon  events.  Theie  are 
two  obvious  differences  that  come  to  mind.  One  is  the  weight  of  the  factor 
(I-COS0)  that  usually  enters  into  the  moment  turn  scattering  rate  [12],  and  the 
other  is  the  frequency  dependence  of  the  scattering  rate  [13J.  The  presence  of 
important  multiphonon  processes  [17]  may  be  such  that  these  relaxation  rates 
become  completely  different  from  each  other. 

The  research  focuses  on  the  properties  of  an  electron  gas  at  the  interface 
between  a  semiconductor  and  another  material  where  band  bending  effects 
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Induced  by  the  external  field  drive  the  conduction  band  below  the  fermi  level 
in  the  vicinity  of  the  surface.  The  electronic  motion  normal  to  the  surface  is 
ccnflried  by  the  electrostatic  potential.  The  confinement  within  this  narrow 
potential  well  causes  the  electronic  energy  levels  to  be  quantized.  The  electrons 
remain  free  to  move  parallel  to  the  surface.  Thus,  we  consider  the  sequence  of 
subbands  with  energy  levels  E(n,k|j)  given  by 


E(n,k||)  =  An  + 


where  ky  is  a  two-dimensional  wave  vector  parallel  to  the  surface,  m*  is  the 
electronic  effective  mass  and  An  is  the  energy  of  the  n111  subband  at  ky  =  0.  The 
subband  energies,  An  and  the  z-dependence  of  the  wavefunctions  i/sn(z)  are 
found  as  the  solution  of  the  eigenvalue  equation 


E2  32 

7 - TT  +  eEoz  ^n(z)  =  An^n(z) 

2m  dz2 


which  has  the  solutions 


0n(z)  =  A  i 


z~  zn 


depicted  in  Figure  I.  In  (3), 


2m*eE„ 


is  a  measure  of  the  thickness  of  the  depletion  layer,  and  zn  —  A/eE0.  The  sub¬ 
band  energy  An  is  fixed  by  the  boundary,  condition  at  z  =  0,  which  requires  that 
Ai(-zn/d)  vanishes.  Thus,  we  have 


An  =  eE0  xn 


where  xt  =  2.338,  x2  =  4.087,  which  are  the  zeros  of  the  Airy  function. 

We  have  explored  the  influence  of  electron-phonon  interactions  on  the 
properties  of  the  electrons  in  these  inversion  layers,  in  particular  the  multipho¬ 
non  processes. 

We  have  generalized  the  electron-phonon  interaction  formulated  by  Evans 
and  Mills  [3  4]  to  the  case  of  an  interface.  We  have  considered  the  scattering 
from  both  the  bulk  and  interface  phonons,  since  it  has  been  noted  [15,18]  that 


Si 


& 


-  82  - 


SECTION  II:  SOLID  STATE 


V(z) 


Figure  1  The  z-dependence  of  Che  wave  functions  of  Che  first  and  second 
lowest  subbands. 


at  interfaces  with  large  dielectric  mismatch  the  optical  phonons  may  have  large 
energies  (670  K). 

Thus,  the  two  frequencies  of  the  interface /optical  phonons  are  given  by 
the  condition 


(4) 


in  which  and  e2(w)  are  the  frequency-dependent  dielectric  functions  for 

the  materials  on  either  side  of  the  interface.  The  frequencies  of  the  interface 
phonons  are  generally  Intermediate  between  the  LO  and  TO  phonon  frequen¬ 
cies  of  each  medium.  The  electrostatic  potential  set  up  by  such  phonons 
satisfies  Laplace’s  equation,  if  one  neglects  the  effects  of  retardation.  Hence, 
the  potential  for  an  interface  phonon  with  a  wave  vector  qy  parallel  to  the  sur¬ 
face  has  the  spatial  dependence  of 


exp 


1  cf||  ’  *11  -  Qj| 


(5) 


where  z  is  the  distance  measured  normal  to  the  interface.  Because  of  this  z- 
dependence,  we  find  that  electrons  trapped  on  one  side  of  the  interface  should 
be  appreciably  scattered  by  the  potential  set  up  by  both  interface  modes. 
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In  contrast,  the  bulk  phonons  decrease  rapidly  in  amplitude  at  the  interface 
boundary,  as  can  be  seen  by  the  following  elementary  argument.  The  bulk 
phonons  of  region  l  satisfy  the  condition  v  •  D  =  0,  by  virtue  of  €i(wLOl)  =  °- 
In  region  2,  if  7^  °>  we  nrnst  have  v2  <t>2  =  Thus  for  z  <  0,  the 

potential  is 


<£2(z)  =  <f>$  exp  £i  •  T||  +  q||Z  j 


Continuity  of  <f>  at  z  =  0  requires  <f>x  =  <t>°.  Continuity  of  D  then  gives 
<^(0)  =  02(O)  =  0,  so  that  the  potential  <f>x  has  the  form 


^i(z)  =  exp  [i^||  •  ?||  J  sin  qzz 


which  vanishes  as  z  approaches  zero.  Thus,  the  potential  due  to  the  high  fre¬ 
quency  bulk  LO  phonons  will  not  penetrate  across  the  interface  to  where  the 
electrons  are  located,  while,  as  argued  above,  the  interface  phonons  do. 

The  frequency  of  the  interface  optical  phonons  is  given  by  u;s  where 


2  A 
u;  =  — 

3  A 


ei(0)+e2(°°) 

C1(0O)+e2(00) 


^(CQ+e^oo) 

e2(oo)+e1(oo) 


f  €i(0)+€2(oo)  |  2 

|  e2(0)+el(oo)  |  2  \ 

\  ei(oo)+e2(oo)  ) 

{  €j(oo)+€2(oo)  J  U 2  j 

f  ‘.(OWoo)  )  I 

f  C2(0)-C2(°0)  ) 

\  e^ooJ-M^oo)  j  1 

i  ex(cx:)+e2{(x>)  1  _ 

in  which  ux  and  w2  are  the  transverse  optical  frequencies  of  the  materials  on 
either  side  of  the  interface.  Likewise,  e^O),  e2(0)  and  e^cc),  e2(oo)  are  the 
zero  frequency  and  high  frequency  limits  of  the  dielectric  constants.  Thus, 
electrons  in  material  1  can  scatter  with  Interface  phonons,  and  the  energy 
transfer  can  involve  large  frequencies,  of  the  order  of  w2.  By  contrast,  the  bulk 
phonon  scattering  rate  may  only  involve  energies  of  the  order  Hwj,  which  could 
be  significantly  smaller,  as  may  be  seen  from  Figure  2. 

The  lifetime  of  an  electronic  state  can  be  analyzed  from  a  study  of  the 
on; -electron  Green's  function,  defined  by 
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Figure  2  The  graphical  solution  of  the  equation  e,(cj)  =  -e2(o>)  for  the  fre 
quency  of  the  Interface  phonons. 


G„VV>  =  ' ««([%»«•  ^(0)1)  <«> 

where  a£((n  and  a^n  are  the  operators  which  annihilate  and  create  electrons  of 
wave  vector  kn,  in  the  subband  n. 
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The  Fourier  transform,  with  respect  to  time,  obeys  the  Dyson  equation 


a;  -  An  - 


Gnn'  (kll>w) 


~  2  2  (^||»**0  Ga"n'  (^||»^)  —  ^nn' 


which  defines  the  self-energy  The  lowest  order  contribution  to  the 


self-energy  is  given  by 


2  ~  2  ^nn"  (^f||)  Mn«n/  ("3||) 

nn'  %  n" 


w  l+NC^jJ-f^kip^n)  N(u>s)+fn«(k|p  qt||) 

✓x  ^  *4“  " 

(j-  An« - -(k|r^,)2-a;s-  ir;  u-  An - —(kip-fr,)2-^--  iq 

2m2  2m2 


in  which  N(x)  and  f(x)  are  the  Bose-Einstein  and  Fermi-Dirac  distribution 
functions,  respectively. 

The  imaginary  part  of  the  diagonal  self-energy  gives  rise  to  the  lifetime  of 
the  electronic  state,  as  calculated  in  the  Born  approximation.  .The  energy  loss 
associated  with  single-phonon  events,  which  scatter  electrons  from  the  state 
(n,k||),  has  been  tabulated  in  [18].  Although  these  calculations  are  in  reason¬ 
able  quantitative  agreement  with  other,  similar  theoretical  work,  they  are  not  to 
be  relied  upon.  As  Kawamoto,  Kalla  and  Quinn  [19]  pointed  out,  the  real  part 
of  the  electronic  self-energy  represents  a  large  fraction  of  the  interband  sub- 
spacings  and  it  is  strongly  k-dependent.  We  find  that  it  is  also  rapidly  varying 
with  u,  which  strongly  renormalizes  the  lifetime,  through  the  wave  function 
renormalization.  For  example,  on  approximating  the  self-energy  by 

eff  _ 

2  (kn-w) 


2(^(1’^)  2 


2(k||,w)  2(k||,u;) 

n.n'  n'n 


w  ^n'  ,  2 

2m  mm 


£ 
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the  lifetime  r  of  the  state  (k||,n)  Is  given,  by 


eff 

Im  £  (kjj.cu) 

nn _ 

d  eff  - 

—  Re  £(k„,a/) 


(12a) 


which  is  to  be  evaluated  at  the  frequency  u,  given  by  the  solution  of  the  equa¬ 
tion 


n^ki?  «ff  _ 

W-  An-  —  =  £(k||.w) 

2m  nn 


(12b) 


Tne  rapid  variation  of  Re  £(kj|,a;)  may  result  in  an  enhancement  of  1/2 

nn 

by  up  to  an  order  of  magnitude  greater  than  the  Born  approximation  value. 
Inclusion  of  higher  order  processes  significantly  changes  these  results.  As  a 
consequence  of  these  strong  manybody  effects,  we  have  adopted  a  multiphonon 
approach  to  the  calculation  of  the  lifetimes. 

The  interface  phonons  are  treated  by  performing  two  canonical  transforma¬ 
tions 


H'  =  Dj+D+HDjDj 


where 


£(V  s’iibib5 

l  <iii 


D2  =  exp  £ 
<iii 


n  <■»  +  <■ 


2mo;, 


and  X  is  a  coupling  constant. 

This  transformation  complicates  the  kinetic  energy  terms  of  the  Hamil¬ 
tonian,  but  does  yield  a  z-dependent  Franck-Condon-like  shift.  This  yieids  part 
of  the  image  potential  at  the  interface.  One  finds  that  if  e00,>e00l  the  electrons 
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are  attracted  to  the  interface  while,  if  £^>6^,  the  electrons  are  repelled. 

The  image  potential  does  result  in  a  change  in  the  z  dependence  of  the 
electronic  wavefunctions,  but  oniy  the  lowest  states  are  appreciably  effected. 
The  states  with  higher  principle  quantum  numbers,  n,  when  evaluated  with  the 
W.K.B.  approximation  do  not  significantly  deviate  from  the  corresponding 
W.K.B.  approximation  for  the  Airy  functions.  As  a  result  of  these  changes,  the 
matrix  changes  involved  in  the  thermalization  process  of  higher  energy  elec¬ 
trons,  do  not  change  by  an  order  of  magnitude,  one  might  infer  from  the 
exponential  dependence  of  the  interface  phonon  amplitudes.  The  variation  is 
limited  to  a  factor  of  2  or  3,  for  reasonable  values  of  the  parameters.  This 
results  in  a  change  of  the  electronic  lifetimes  by  a  factor  of  9.  This  is  of  the 
right  order  of  magnitude  needed  to  reconcile  the  theory  and  experimental  data. 
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1.  OBJECTIVE(S) 

Standard  CVD  methods  for  the  growth  of  silicon  films  require  that  the 
substrate  be  heated  to  temperatures  that  cause  the  decomposition  of  silane 
(Ts  >  600°C).  Unfortunately,  these  high  temperatures  limit  the  incorporation 
of  hydrogen  in  the  film.  This  is  detrimental  to  the  electronic  properties  of  the 
amorphous  films.  A  high  concentration  of  unpaired  spins,  a  consequence  of  low 
[H],  pins  the  Fermi  level.  However,  with  incorporation  of  sufficient  hydrogen 
the  defect  concentration  can  be  lowered  enough  to  allow  the  Fermi  level  to  be 
adjusted  by  substitutional  doping.  Therefore,  the  uncoupling  of  the  gas 
temperature  from  substrate  temperature,  to  allow  for  a  lowering  of  the 
nucleation  temperature,  has  been  of  great  interest  in  the  growth  of  silicon.  In 
the  growth  of  III-V  and  II- VI  compound  semiconductors  there  are  temperature 
problems  associated  with  substrate  decomposition  (InSb,  InP),  film  breakdown 
(HgCdTe)  and  dopant  migration. 

Two  methods  that  have  been  widely  used  to  separate  growth  temperature 
from  gas  phase  chemistry  have  been  plasma  assisted  CVD  and  UV  photolysis 
CVD.  Both  techniques  have  Inherent  drawbacks;  for  example  a  plasma  can 
cause  sputtering  damage  to  the  film,  and  UV  sources  exhibit  both  high  cost  and 
low  efficiencies.  However,  cw  C02  lasers  offer  an  inexpensive  and  efficient 
source  of  photons  which  could  be  used  to  lower  the  substrate  temperatures  .  It 
had  been  reported  in  the  literature  that  continuous  vibrational  excitation  of 
silane  gas  above  a  heated  substrate  would  lower  the  temperature  necessary  for 
amorphous  silicon  deposition  [1-4].  We  felt  that  if  it  were  possible  to  lower  the 
temperature  of  deposition  by  vlbrationally  exciting  gas  phase  molecules,  that  we 
would  be  able  to  isolate  the  substrate  requirements  for  many  reactions  that 
have  temperature  limitations.  In  addition,  we  felt  that  it  would  be  important  to 
have  a  better  understanding  of  what  role  the  vibrational  excitation  played  in  the 
chemical  reactions  occurring  in  the  deposition  process. 

We  thought  that  the  best  system  to  begin  our  studies  was  on  the  simple 
decomposition  of  silane  to  silicon.  This  system  has  been  widely  studied  by 
classical  CVD  and  is  an  important  technological  area.  In  addition,  since  silane 
does  not  undergo  an  electronic  transition  until  the  vacuum  ultravtolet,  the  C02 
laser  appears  to  be  the  best  choice  for  excitation  of  the  gas  species. 
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2.  SUMMARY  OF  RECENT  PROGRESS 

We  have  studied  in  detail  the  irradiation  of  silane  with  energy  supplied  by 
a  cw  C02  laser.  We  have  determined  that  vibrationally  excited  silane  rapidly 
transfers  its  energy  to  the  surrounding  molecules  and  heats  the  gas  phase. 
Therefore,  contrary  to  previous  results,  we  have  shown  that  vibrationally 
excited  silane  molecules  do  not  deposit  films  at  temperatures  lower  than  those 
necessary  for  classical  CVD.  In  addition,  we  have  shown  that  using  proper 
conditions,  the  C02  laser  can  increase  the  temperature  of  the  gas  to 
temperatures  where  homogeneous  gas  phase  chemistry  occurs  with  silane.  The 
most  likely  Initial  step  is  the  formation  of  S1H2  which  reacts  readily  with  the 
surrounding  gas  molecules.  In  a  system  containing  a  high  concentration  of 
silane,  polysilanes  are  formed  in  the  gas  phase.  These  polysilanes  can  travel 
long  distances  within  the  reactor  ( >  l  cm)  and  deposit  amorphous  silicon  films 
at  temperatures  as  low  as  200°C.  In  systems  containing  other  reactive  species, 
for  example  ammonia  or  nitrous  oxide,  the  S1H2  can  form  gaseous 
intermediates  capable  of  forming  silicon  nitride  or  oxide  films  at  similar  low 
temperature. 


3.  STATE  OF  THE  ART  AND  RESEARCH  DETAILS 

All  experiments  were  performed  in  a  CVD  reactor  equipped  with  ZnSe 
windows  for  the  introduction  and  monitoring  of  the  ir  laser  light. 

In  initial  experiments  we.  attempted  to  decompose  silane  using  an 
unfocused  laser  beam  (the  laser  operated  In  multi-line  TEMqq  mode).  The 
laser  beam  traversed  the  chamber  1  cm  above  the  substrate.  The  substrate  was 
held  to  a  temperature  below  700°  K.  Neither  film  growth  nor  powder 
generation  was  observed  for  laser  intensities  up  to  100  watts  cm-2  and  silane 
pressures  from  1  to  6  torr.  The  addition  of  5  torr  of  argon  gas  to  reduce  the 
thermal  conductivity  of  the  reactor  gas  caused  no  noticeable  change.  We 
proceeded  to  assess  the  effects  of  Increased  power  densities  on  the 
decomposition  process.  The  laser  beam  was  focused  with  a  20  cm  F.L.  front 
surface  spherical  mirror  before  its  entry  into  the  reaction  zone.  The  beam  waist 
within  this  zone  was  calculated  to  be  « 0.3  mm  In  radius.  With  power  densities 
up  to  17  kw  cm"2  no  decomposition  or  deposition  was  observed.  The  beam 
path  was  then  altered  so  that  it  impinged  on  the  leading  edge  of  the  substrate 
Fig.  1  to  produce  a  small  incandescent  spot.  The  substrate  temperature  thus 
decreased  with  circular  symmetry  away  from  this  spot. 
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The  resulting  growth  rings  Fig.  2  were  semicircular  with  no  perturbation  to  the 
ring  uniformity  in  the  regions  Just  below  the  passing  laser  beam.  This  result 
indicated  that  a  grazing  incidence  laser  beam  neither  promoted  growth  where 
the  substrate  was  too  cold  nor  altered  the  growth  rate  where  the  substrate 
temperature  was  sufficient  for  normal  cvd.  Although  there  had  been  reports  in 
the  literature  that  vibrationally  excited  silane  could  decompose  at  a  low 
substrate  temperature,  probably  the  surface  substrate  temperature  had  been 
raised  beyond  that  indicated  [1-4].  Under  the  conditions  within  the  reactor, 
similar  to  those  previously  reported,  the  mean  free  path  will  be  very  short  (less 
than  l  micron)  and  collisions  with  gas  molecules  will  cause  energy  transfer  and 
general  heating  of  the  gas.  The  hot  gas  can  then  exchange  heat  very  rapidly 
upon  collision  with  the  substrate  causing  a  rise  in  rhe  .  urface  temperature. 
These  initial  experiments  led  us  to  believe  that  in  order  for  pure  gas  phase 
modifications  to  occur  we  would  have  to  induce  homogeneous  chemical 
reactions  in  the  gas  phase,  and  that  the  new  gas  phase  species  could  react  at 
different  substrate  temperatures.  We  performed  a  series  of  experiments  to 
determine  if,  and  under  what  conditions  of  gas  composition,  pressure,  and  laser 
intensity  silane  homogeneous  reactions  would  occur. 

Visually  detectable  decomposition,  noted  by  the  formation  of  a  light  yellow 
film  on  a  silica  substrate,  occured  with  a  focused  50  w  laser  beam  at  10.5  torr  of 
silane  and  with  a  focused  40  w  laser  beam  at  7  torr  of  silane  and  20  torr  of 
argon.  The  addition  of  the  argon  served  to  lower  the  thermal  conductivity  of 
the  gas  within  the  reactor  and  In  turn  lead  to  the  lowering  of  the  power 
required  for  the  threshold  of  detectable  deposition.  We  proceeded  to  perform 
another  series  of  runs  with  argon  in  excess  at  a  pressure  of  100  torr,  varying 
the  silane  pressure  and  incident  laser  power  until  film  formation  was  detected. 


.-■rv-js--:-  mssSS 


Figure  2.  Growth  rings  showing  regular  deposition 

These  experiments  were  performed  to  aid  In  the  calculation  of  the  gas  phase 
decomposition  temperature.  To  perform  this  calculation  a  boundary  condition 
was  required.  This  was  supplied  by  a  silica  substrate  (2.5  cm  i.d.  0  cm  length) 
held  coaxial  with  the  laser  beam.  It  is  important  to  note  the  Importance  of  the 
position  of  the  beam  focus  along  the  line  connecting  the  input  and  output  ZnSe 
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windows  of  the  reactor.  If  the  focus  fell  towards  the  rear  of  the  reactor  the 
silane  gas  attenuated  the  incident  laser  beam  too  strongly  for  the  power 
densities  required  for  decomposition  to  be  reached.  In  our  configuration  the 
focus  fell  «3cm  from  the  front  window.  The  silica  tube  was  set  in  a  cradle  and 
the  cradle  was  adjusted  so  that  its  midpoint  coincided  with  the  position  of  the 
focus.  The  temperature  of  the  silica  tube  was  measured  by  a  thermocouple 
attached  to  a  silica  hook  formed  on  the  tube  and  held  in  intimate  thermal 
contact  by  gallium  metal  Fig.  3.  It  was  found  that  noticeable  film  production 
took  place  when  the  substrate  reached  a  temperature  ~  100°C  above  reactor 
ambient.  This  value  was  used  as  the  boundary  condition  for  the  subsequent  gas 
decomposition  temperature  calculations.  The  following  are  the  set  of 
conditions  for  silane  pressure  and  laser  power  that  lead  to  a  noticeable  rate  of 
film  growth  at  the  previously  determined  boundary  condition. 

Silane  Pressure  (torr)  Laser  Power  (watts) 


9.3 

3 

8.7 

8.2 


18 

20 

28 

34 


At  these  conditions  we  immediately  observed  disilane  and  trisilane,  in  the 
gas  phase,  by  quadrupole  mass  spectrometry.  If  the  laser  power  was  increased, 
the  higher  resultant  gas  temperatures  lead  to  the  rapid  formation  of  powder 
deposits  on  the  substrate  and  on  the  reactor  walls.  Infrared  spectroscopy  of 
these  deposits  showed  strong  SIH2  and  SIH3  peaks  indicative  of  polysilane. 
When  the  laser  power  was  slightly  reduced,  and  the  substrate  temperature 
raised,  it  was  possible  to  grow  films  at  a  substrate  temperature  as  low  as  200°C. 
These  films  could  be  grown  even  when  the  laser  beam  was  more  than  a 
centimeter  from  the  substrate.  These  distamces  are  far  greater  than  those 
possible  for  the  transfer  of  energetically  excited  species,  but  rather  correspond 
to  the  diffusion  of  a  reasonably  stable  chemical  species  created  at  the  laser 
focus. 

4.  CALCULATION: 

We  chose  to  approximate  the  system  as  an  infinite  cylinder  ignoring  the  z 
dependence  in  the  heat  equation,  since  the  horizontal  extent  of  the  substrate 
was  several  times  as  great  as  its  radius.  This  is  similar  to  the  case  of  an  infinite 
cylinder  with  a  central  uniformly  heated  core.  The  equation  governing  our 
model  then  reduces  to  that  of  the  steady  state  heat  equation  in  r: 

d2T/dr2  +  l/r(dT/dr)+q/K  =  0 

q  =  heat  source  strength  per  unit  volume  per  unit  time 

K  =  thermal  conductivity 

A  Pascal  program  was  written  which  numerically  integrated  the  heat 
equation  subject  to  the  condition  that  at  r  =  R,  the  radius  of  the  tube,  the 
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temperature  converges  to  that  given  by  the  boundary  condition.  A  total  of  500 
integration  steps  were  performed  along  the  radius  extending  from  the 
theoretical  beam  waist  to  r  =  R.  The  attenuation  coefficient,  which  was 
required  for  the  determination  of  power  loss  within  a  unit  length  of  the  beam 
path,  was  taken  from  curves  previously  published  for  the  pressure  dependent 
absorption  coefficient  for  silane  [3].  Its  value  was  estimated  to  be  ~  4.9x10"  4 
/cm_1torr-1  a  moderately  strong  absorption.  Since  we  were  working  in  dilute 
solution  with  argon  gas  we  could  use  the  thermal  conductivity  of  argon 
(1.75X10-4  watts /( cm2) (°C/cm).  Corrections  were  performed  for  the  affects  of 
temperature  on  the  values  of  the  thermal  conductivity  and  attenuation 
coefficient.  Further,  local  temperature  changes  along  1  cm  increments  of  the 
tube's  major  axis  were  used  to  correct  each  elements  thermal  conductivity  and 
in  turn  its  attenuation  coefficient.  This  lead  to  a  more  realistic  value  for  the 
actual  power  loss  in  the  element  coinciding  with  the  position  of  the  focus.  The 
correction  made  to  the  thermal  conductivity  Is  given  by: 

K  =  K^T^/To)05 

where  =  the  thermal  conductivity  at  20°C.  the  correction  for  the 
attenuation  coefficient  is  given  by: 

&  =  /?o(To/TtocaJ) 

where  0O  —  derived  attenuation  coefficient  at  20°C  where  we  assume  a  Charles’ 
Law  affect  on  the  concentration  of  absorbing  species. 

The  program  ran  up  the  beam  power  in  increments  of  0.1  watts  up  to  the 
powers  listed  above.  This  allowed  the  continuous  correction  of  the  affects  of 
the  rise  in  temperature  and  the  determination  of  the  temperature  at  the  focus 
as  a  function  of  input  power.  A  plot  of  the  calculations  is  shown  in  Fig.  4.  The 
results  of  our  calculations  for  the  various  experiments  performed  are  consistent 
with  an  average  decomposition  temperature  for  the  gas  equal  to  ~1130K. 

5.  DISCUSSION: 

We  believe  that  a  strictly  pyrolytic  process  is  taking  place  when  a 
sufficiently  Intense  cw  ir  laser  beam  Interacts  with  a  silane  gas  mixture  allowing 
for  a  rise  in  gas  temperature  above  HOOK.  The  direct  photolysis  of  silane  by 
UV  photolysis  is  very  difficult  because  of  the  far  UV  character  of  the 
absorption.  Photolysis  can  occur  with  an  ArF  excimer  laser  (photon  energy  8.4 
eV);  however,  the  indications  are  that  this  is  either  a  muitiphoton  process  or  a 
multistep  process  with  the  creation  of  polysilanes  and  subsequent 
decomposition  via  combination  with  silane  [5],  In  the  case  of  ir  laser 
excitation,  a  single  P(20)  ir  photon  (^s.leV)  is  far  to  low  in  energy  to  cause 
photolysis.  In  order  for  decomposition  to  occur  it  would  be  necessary  to  have 
the  rapid  absorption  of  many  photons.  Previous  studies  have  shown  that 
intramolecular  relaxation  of  vibrational  energy  to  an  excited  distribution  of 
states  is  very  rapid  [6].  Therefore  it  is  unlikely  that  sufficient  energy  can  be 
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deposited  by  a  cw  laser  and  collected  in  the  state  necessary  for  bond  fission, 
before  it  is  transferred  by  collision  to  a  bath  molecule.  Further,  it  is  not  clear 
that  the  required  resonance  absorption  condition  still  exists  for  the  absorption 
to  the  v  >  2  vibrational  states.  Purely  photolytic  processes  have  been  observed 
for  silane  using  and  ir  C02  laser  source  but  these  have  involved  collisionless 
multiphoton  dissociations  using  a  TEA  laser  [7].  In  this  case  the  intense 
photon  bombardment  creates  a  quasi  continuum  where  subsequent  photon 
absorptions  lead  to  dissociation. 

Our  results  Indicate  that  at  gas  temperatures  below  some  critical  value 
(Tg  =  HOOK)  there  is  no  evidence  of  rapid  decomposition  sufficient  of  realistic 
film  growth.  Previous  researchers  have  shown  this  to  be  true  in  the  case  of  the 
less  stable  disilane  molecule  [8j.  Under  reactor  pressures  of  between  1  and  10 
torr  the  mean  free  path  for  the  excited  silane  molecule  is  ^1  to  10  microns. 
Therefore,  the  energy  absorbed  by  the  silane  molecule  will  be  equilibrated  after 
traveling  <1  mm.  However,  if  the  substrate  cannot  dissipate  the  heat  gained 
from  the  gas  phase,  its  temperature  could  rise  substantially  allowing  for 
deposition  at  the  surface.  We  believe  that  the  relatively  volatile  and  stable 
dkilane  and  trisilane,  observed  in  the  gas  phase,  and  perhaps  higher  order 
polysilanes  are  responsible  for  the  formation  of  amorphous  silicon  films  at  a 
reduced  substrate  temperature. 

6.  CONCLUSION  AND  FUTURE  DIRECTIONS 

It  seems  most  likely  that  the  first  stage  in  the  laser  pyrolysis  is  to  create  the 
highly  reactive  SIH2  species: 

SiH4  SIH2  +  H2 

In  a  system  with  only  silane  and  inert  gas,  polysilanes  will  be  rapidly  formed 
due  to  collision  with  S1HV  If  other  gas  species  are  present,  then  difierent 
intermediates  will  be  formed  and  new  materials  can  be  deposited  on  the 
substrate  at  a  low  temperature.  In  order  to  verify  this  hypothesis  we  have 
started  work  on  the  low  temperature  growth  of  silicon  nitride.  Since  the 
thermal  conductivity  of  ammonia  is  very  close  to  <ugon,  no  additional  dilution 
is  necessary  to  reach  the  cracking  temperature  for  the  creation  of  active  silicon 
species.  Using  laser  power  and  gas  pressure  conditions  similar  to  those  for  pure 
silane  we  were  able  to  deposit  either  a  transparent  film  or  white  powder  from 
the  NH3/SiH4  mixture  when  the  ratio  Is  *siO  or  greater.  FTIR  measurements 
on  the  sample  show  that  materials  grown  at  substrate  temperatures  of  200°C  or 
less  contains  SixNr  The  Si-N  stretch  is  the  major  feature  of  the  spectrum: 
however,  there  Is  observable  hydrogen,  primarily  detected  as  a  N-H  stretch. 
This  was  expected  since  it  is  very  difficult  to  to  remove  the  last  traces  of 
hydrogen  from  nitrides. 

It  has  been  shown  that  by  exciting  gas  phase  silane  with  a  C02  laser  we  are 
able  to  create  active  gas  phase  species  which  can  prepare  low  temperature 
materials.  There  are  many  areas  of  this  research  to  be  explored  both  from  a 
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fundamental  perspective  as  well  as  a  technological  one: 

(a)  Careful  characterization  of  the  films  grown; 

(b)  In  situ  detection  of  gas  phase  species  (ideally  FTER); 

(c)  Testing  with  additional  species  either  as  dopants  or  main  reactants  with 
silane;  and 

(d)  Expansion  to  III-V  materials,  ideal  because  of  the  access  to  excitation  of 
the  group  V  hydrides  with  the  C02  laser. 
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E.  RESONANT  INTERACTIONS  IN  CRYSTALS  AT  X-RAY 
WAVELENGTH  SURFACES 

Professors  G.  Schaefer  and  E.E.  Kunhardt 

Unit  SS7-5 


1.  OBJECnVE(S) 

The  objective  of  this  project  Is  to  Investigate  concepts  for  obtaining  a  non- 
dlvergent  beam  of  x-ray  radiation  with  enhanced  coherence  from  the  interac¬ 
tion  of  Ewald  modes  and  active  oscillators  in  a  crystal. 


2.  SUMMARY  OF  RECENT  PROGRESS 

In  general,  there  are  two  ways  to  enhance  the  coherence  properties  of  radi¬ 
ation  emitted  by  atomic  oscillators:  (1)  to  alter  the  properties  of  spontaneous 
emission  through  the  interaction  of  the  oscillators  wtth-a  resonant  cavity;  and 
(2)  to  utilize  induced  emission.  We  have  developed  a  new  concept  for  obtain¬ 
ing  a  non-divergent  beam  of  x-ray  radiation  with  enhanced  coherence  proper¬ 
ties.  The  approach  utilizes  the  Interaction  of  active  oscillators  and  Ewald  modes 
in  a  crystal.  Using  the  high  Q  that  can  be  achieved  in  a  crystal,  the  x-ray  emis¬ 
sion  can  be  channeled  into  a  limited  number  of  allowed  modes  (Kossel  modes); 
this  increases  the  lifetime  of  the  oscillators  and  thus  enhances  the  coherence  of 
the  radiation.  We  have  been  Investigating  various  methods  for  enhancing  the 
Q  factor  of  the  crystal  cavity  and  for  extracting  the  radiation  from  the  crystal. 
These  Issues  will  be  discussed  In  the  next  section. 


3.  STATE  OF  THE  ART  AND  PROGRESS  DETAILS 

In  recent  years,  the  concept  of  inhibited  or  enhanced  spontaneous  emission 
in  cavities,  first  outlined  by  Purcell  [l]  at  Radio  Frequencies,  has  been 
extended  by  Stehle  [2]  and  Kleppner  [3]  to  the  near  IR  and  visible  spectral 
regime.  This  concept  has  been  verified  by  several  experiments  [4-8]. 

The  basic  quantity  describing  the  interaction  between  an  individual  atomic 
transition,  (n-m),  and  electromagnetic  radiation  of  the  frequency 
y=»(En- Ej„)/h,  is  the  Einstein  coefficient  for  absorption  or  Induced  emis¬ 
sion.  The  Einstein  coefficient  for  spontaneous  emission.  A,  in  addition 
depends  on  the  availability  of  states  for  the  electromagnetic  radiation  of  the  fre¬ 
quency  u  expressed  as  the  density  of  photon  states,  p{u): 

Anm  =  p(i/)(h^)Bmn 
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The  values  for  A,  usually  used  to  characterize  the  transition  probability  of  a 
specific  transition,  are  based  on  the  availability  of  states  in  free  space  for  dipole 
radiation  p3=8ir  u2 /c3: 

87T 1/2  ■ 

A  =  - ht'B 


The  density  of  states  in  a  cavity,  pc(v),  is  quite  different.  If  the  cavity  is  tuned 
with  its  lowest  order  mode  to  the  frequency  u  of  the  atomic  transition  and  the 
quality  factor  of  the  cavity  is  Q,  then  the  density  of  states  is  pc=8ir  v*Q/c3  and 
the  transition  probability  in  the  cavity,  Ac,  is  enhanced: 


Ac  =  Q 


hi/B  =  QAS 


If  the  cavity  dimensions  are  small  compared  to  half  the  wavelength  of  the 
atomic  radiation  then  the  transition  probability  is  decreased. 

The  same  effect  is  observed  if  the  atomic  system  is  placed  Inside  a 

waveguide  so  that  the  frequency  of  the  atomic  transition,  u,  is  close  to  the 

\ 

cutoff  frequency,  ue,  of  the  waveguide.  At  the  cutoff  frequency  the  waveguide 
acts  as  a  resonator  with  quality  factor  Q  and  the  transition  probability  in  the 
guide  Ag  becomes: 

Ag(  U—1/qq)  —  QqAj 

Below  cutoff  the  transition  probability  approaches: 

Ag(i'<i'c0)  =  "q”  As 

and  above  cutoff  the  transition  probability  decreases  with  u: 


A.c>^)  cc  ——  A, 


up  to  the  cutoff  frequency  for  the  next  higher  order  mode  i/cl,  at  which  the 
waveguide  again  acts  as  a  cavity  with  the  quality  factor  Q{  for  this  mode.  Once 
a  larger  number  of  modes  can  be  excited,  the  resonance  effect  for  the  indivi¬ 
dual  mode  becomes  less  pronounced,  and  finally,  the  density  of  modes 
approaches  the  value  for  free  space. 

The  concept  of  Inhibited  and  enhanced  spontaneous  emission  can  also  be 
applied  to  the  emission  of  x-ray  radiation  inside  a  crystal  lattice.  Such  radiation 
(Kossel  radiation)  is  only  emitted  as  specific  modes  of  propagation,  the 
Bornnann  modes.  Although  the  basic  features  of  the  Borrmann  effect  are  well 
known,  they  are  briefly  reviewed  here  [9-10],  They  are  as  follows:  a  parallel¬ 
sided,  monocrystai  slab  with  lattice  planes  perpendicular  to  the  slab  faces  is 
illuminated  by  a  beam  of  x-radlation  at  an  angle  9  with  respect  to  the  planes. 
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The  field  distribution  inside  the  crystal  excited  by  the  Incident  beam  can  be 
described  as  a  superposition  of  what  we  shall  denote  as  Ewald  modes.  These 
are  the  sets  of  modes  that  are  permitted  in  a  system  consisting  of  a  periodic 
array  of  oscillators.  (See  [11]  for  a  discussion  of  the  dispersion  surfaces 
corresponding  to  these  modes.)  In  general,  these  modes  have  a  complex  wave 
vector,  and  they  are  strongly  attenuated,  so  that  the  energy  in  the  incident 
beam  is  quickly  absorbed.  When  the  angle  of  incidence,  d,  is  set  to  the  Bragg 
angle,  one  of  the  modes  has  an  almost  real  wave  vector.  This  mode,  the 
Borrmann  mode,  may  be  considered  as  the  superposition  of  the  incident  and 
diffracted  beams,  with  corresponding  propagation  vectors,  Kq  and  KH  that  are 
directed  to  reciprocal  lattice  points  (r$ps)  O  and  H  on  the  surface  of  the  Ewald 
sphere  [11].  For  this  angle  of  Incidence,  these  beams  have  essentially  the  same 
amplitude  and  wavelength.  The  resulting  wave  (the  two  beam  Borrmann 
mode)  has  a  standing  wave  pattern,  with  the  periodicity  of  the  lattice,  in  the 
direction  perpendicular  to  the  lattice  planes,  and  propagates  in  a  direction  paral¬ 
lel  to  the  lattice  planes.  The  minimum  field-planes  are  coincident  with  the  lat¬ 
tice  planes,  which  explains  why  the  mode  suffers  little  absorption  (almost  real 
wavenumber).  Thus,  the  wave  (mode)  propagates  through  the  crystal  and 
emerges  at  a  point  parallel  to  where  the  Incident  beam  struck  the  crystal. 
Again,  note  that  this  is  a  traveling  mode  in  the  direction  parallel  to  the  lattice 
planes. 

Increasing  the  number  of  r#ps  on  the  surface  of  the  sphere  to  n  >  2 
results  In  the  n-beam  Borrmann  mode.  In  such  cases,  the  diffracted  beams  are 
highly  collimated,  as  well  as  monochromatic,  and  their  phases  are  displayed  in 
the  form  of  asymmetries  of  the  diffracted  intensity  profiles  [13].  The  absorp¬ 
tion  coefficients  of  selected  modes  of  propagation  can  be  reduced  to  values 
much  smaller  than  those  attainable  In  the  2-beam  case  [14].  Calculations  indi¬ 
cate  [15]  that,  for  transmission  of  CuKa  radiation  through  Ge  crystals,  the 
absorption  coefficients  of  four  of  the  twelve  modes  of  propagation  in  a  8-beam 
interaction  can  be  reduced  from  the  average  of  352  cm" 1  to  values  smaller  than 
one  cm"1 .  Calculated  values  of  the  excitations  of  all  twelve  modes  of  propaga¬ 
tion  of  x-rays  and  the  corresponding  absorption  coefficients  are  listed  in  Table  I. 

Quantitative  experimental  verification  has  not  yet  been  reported,  but  qualitative 
photographic  results,  described  in  [15],  are  in  reasonable  agreement  with  the 
above  calculations.  Detailed  verification  of  the  latter  will  make  possible  x-ray 
Investigations  of  highly  absorbing  materials  using  relatively  long  wavelengths, 
as  well  as  the  more  general  class  of  experiments  which  Involve  long  paths  of  x- 
ray  beams  within  crystals. 

In  the  analysis  of  Ewald  modes,  a  plane  wave  description  is  normally  used 
and  the  oscillators  are  assumed  to  be  nassive  [11].  The  plane  wave  assumption 
was  removed  in  the  analysis  of  Kato,  [17]  who  instead  assumed  spherical  wave 
fronts  for  the  incident  radiation.  This  resulted  in  a  better  description  of  the 
experimental  conditions  which  invariably  involve  beams  of  finite  width  and 
non-zero  divergence.  A  Gaussian  incident  beam  profile  would  be  more  realis¬ 
tic. 
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Table  I 

/i's  and  excitations  at  the  exact  6-beam  pt. 
for  CuKc*!  radiation 


Mode 

ficm~ 1 

Excitation 

1 

0.02 

1/12 

2 

0.16 

1/18 

3 

0.16 

1/9 

1/3 

4 

0.48 

1/12 

5 

38.72 

1/18 

6 

38.72 

1/9 

7 

80.02 

1/18 

1/3 

8 

80.02 

1/9 

0 

85.43 

1/12 

10 

1,186.18 

1/12 

11 

1,377.79 

1/18 

1/3 

12 

1,377.79 

1/9 

=  352.12 

• 

To  our  knowledge,  none  of  the  descriptions  of  the  propagation  of  Ewald 
modes  in  a  crystal  takes  into  account  the  possibility  that  the  oscillators  may  be 
active.  That  is,  they  may  exist  in  a  number  of  different  states  whose  popula¬ 
tion,  to  a  certain  extent,  can  be  entirely  controlled  (by  excitation  with  an  elec¬ 
tron  beam,  or  x-rays  of  shorter  wavelength,  for  example).  Moreover,  some  of 
the  oscillator  transitions  may  be  resonant  with  the  Incident  radiation  (Ewald 
modes). 

From  the  above  description,  it  is  evident  that  for  Ka  radiation,  the  crystal 
lattice  acts  as  a  system  of  waveguides  with,  in  general,  a  large  number  of 
allowed  modes.  The  number  of  waveguides  for  which  the  cutoff  frequency  of 
the  lowest  order  modes  is  below  the  frequency  of  the  x-ray  radiation  of  course 
depends  on  the  lattice  constant  and  the  energy  levels  of  the  emitting  atom. 
The  frequency  of  Ka  radiation  Is,  to  the  first  order,  proportional  to  Z2  where  Z 
is  the  atomic  number.  For  crystals  composed  of  only  one  kind  of  atom  there  is 
only  little  choice  for  the  right  combination  of  wavelength  and  lattice  constant, 
which  means  that  only  one  or  a  few  waveguides  exist  with  <  u.  Some  crystals 
of  low  Z  materials,  in  fact,  do  not  have  any  allowed  Borrmann  modes.  Com¬ 
pound  crystals,  doped  crystals,  and  doped  compound  crystals,  however,  provide 
a  large  variety  of  combinations  of  lattice  constants  and  wavelength  so  that  a 
combination  can  be  found  for  which  the  x-ray  frequency  is  little  above  the 
cutoff  frequency  of  a  low  order  allowed  mode.  Of  special  interest,  of  course, 
are  modes  with  extremely  low  absorption  coefficients  as  shown  in  Table  I 
(modes  11  and  12  ).  These  are  the  modes  which,  at  cutoff,  correspond  to  a 


SECTION  II:  SOLID  STATE 

i 

very  high  quality  factor  Q. 

There  are  also  methods  to  provide  some  fine  adjustment  of  the  lattice  con¬ 
stant.  A  method  used  during  fabrication  would  be  to  adjust  the  doping  density 
of  some  doping  material  (not  necessarily  the  atoms  used  for  x-ray  emission). 
Fine-tuning  during  the  experiment  may  be  accomplished  through  temperature 
or  pressure  changes. 

Another  method  to  obtain  high  Q  factors  is  to  generate  a  cavity  which 
means  to  provide  a  coupling  mechanism  between  a  Borrmann  mode  and  the 
equivalent  mode  with  opposite  direction  of  propagation.  Assuming  that  one 
finds  a  combination  of  emitter  and  lattice  so  that  the  frequency  u  >  uc  is  only 
fairly  close  to  the  cutoff  frequency,  uc,  of  a  low  loss  Borrmann  mode.  In  this 
case,  the  wavelength  of  the  guided  wave,  \v  is  much  larger  than  the 
wavelength,  X,  of  the  x-ray  radiation.  A  periodic  structure  with  periodicity 
d=\g/2  and  planes  perpendicular  to  the  propagation  direction  of  the  Borrmann 
mode  would  then  create  a  cavity  for  the  x-ray  radiation.  Such  a  periodic  struc¬ 
ture  could,  for  example,  be  generated  through  a  super  lattice. 

Another  way  to  provide  feedback  between  the  Borrmann  modes  with 
opposite  direction  of  propagation  Is  x-ray-phonon  scattering.  If  we  select  a 
combination  of  lattice  constant  and  x-ray  wavelength  such  that  the  lattice  con¬ 
stant  is  only  slightly  larger  than  half  the  wavelength  a=(i+e)~  with  (e<<i) 

2 

then  the  propagation  constant  of  the  propagating  Borrmann  mode  becomes  very 
small  and  can  be  the  same  order  of  magnitude  as  the  propagation  constant  of 
phonons.  In  this  case,  it  is  possible  to  launch  an  acoustic  wave  with  the  right 
frequency  in  the  opposite  direction  to  the  Borrmann  mode  such  that  180°  back 
scattering  of  the  propagating  Borrmann  mode  occurs.  A  standing  wave  pattern 
(two  acoustic  waves  in  opposite  directions  in  combination  with  a  low  loss 
Borrmann  mode  (see  Table  I)  in  a  crystal  close  to  cut-off,  therefore,  provides 
an  x-ray  cavity  with  a  very  high  quality  factor  Q.  It  should  be  mentioned  here 
that  Rivlin  [18]  estimates  the  quality  factor  of  Bragg  modes  in  a  perfect  crystal 
to  approach  10s. 

There  are  two  important  consequences  resulting  from  the  interaction  of  x- 
ray  emitters  with  resonant  cavities  which  allow  oscillation  of  only  a  few  modes. 

1.  First,  the  cavity  provides  a  channeling  of  the  x-ray  radiation  into  the  mode 
of  the  cavity,  and 

2.  Second,  the  external  observable  lifetime  of  the  radiation  becomes  dom¬ 
inated  by  the  radiating  decay  of  the  cavity  rather  than  the  lifetime  of  the 
atomic  state.  The  interaction  of  Borrmann  mode  based  cavities  in  crystals 
with  the  active  oscillators,  therefore,  provides  an  efficient  method  to  gen¬ 
erate  beams  of  characteristic  x-ray  radiation  with  improved  spatial  and  tem¬ 
poral  coherence. 
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1.  OBJECTTVE(S) 

In  recent  years  the  general  area  of  robust  detection  and  estimation  has 
received  considerable  attention.  In  this  context,  "robustness"  was  meant  to 
signify  insensitivity  to  small  deviations  from  the  assumptions  about  the  noise 
model.  Freque  tly,  a  robust  method  was  described  as  being  a  procedure  with 
uniformly  good  behavior  within  a  family  of  contaminated  distributions  [1-9], 
However,  the  exact  meaning  of  "uniformly  good  behavior"  and  "contaminated 
family,”  though  defined  in  distinct  mathematical  terms  of  a  give  situation, 
becomes  in  effect  ambiguous  for  the  data  processor  to  which  the  input  is  a 
stream  of  observables  with  unknown  and/or  slowly  time-varying  statistics. 

To  alleviate  this  vagueness,  further  research  efforts  were  concentrated  on 
the  stability  of  various  stochastic  distances  {10-15]  to  give  these  properties  a 
qualitative  as  well  as  quantitative  measure.  Nevertheless,  the  stability  does  not 
always  Imply  good  performance  in  terms  wf  the  expected  value  of  a  cost 
functional;  moreover,  these  stochastic  distances  are  usually  too  abstract  to  be 
easily  related  to  the  system  design. 

To  bypass  the  problems  associated  with  the  differences  generated  between 
theoretical  and  practical  conceptualizations  of  robustness,  one  of  the  major 
objectives  of  the  present  research  effort  was  to  introduce  descriptors  of 
robustness  which  leas  naturally  to  the  design  of  practical  detectors  and 
estimators.  A  major  breakthrough  In  the  theory  of  robust  detection  and 
estimation  was  also  achieved  [31,32].  Specifically,  two  approaches  to  robust 
procedures  which  Include  dependent  sampling  were  developed. 

It  was  recognized  by  Chadwick  and  Kurz  [10]  that  in  low  signal- to- noise 
ratio  (SNR)  environments,  nonparametric  detectors  should  operate  in  a 
sequential  mode  to  insure  an  Increase  In  the  information  rate  while  maintaining 
a  constant  error  rate.  Brownstein  and  Kurz  [  17]  demonstrated  that  group 
rather  than  sample-by-sample  updating  of  the  data  is  necessary  for  efficient 
operation  of  detectors  in  a  sequential  mode.  They  also  demonstrated  that 
stochastic  approximation  algorithms  represent  a  useful  tool  for  an  efficient 
operation  of  adaptive  detectors  [  13,18-221.  The  natural  objective  of  the 
ongoing  research  is  to  investigate  the  use  of  variable  threshold  sequential 
detectors  as  practical  and  efficient  means  of  truncation  of  the  decision  process. 
Also,  a  class  of  sequential  partition  detectors  that  include  stationary  and  quasi¬ 
stationary  Markov  processes  was  studied  In  detail.  Applications  to  detection  of 
underwater  sounding  data  were  stressed  (3ij. 
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With  the  advent  of  satellite  communications,  the  problem  of  data 
transmission  through  satellite  channels  became  particularly  important  [23-20]. 
Some  of  the  techniques  used  In  robust  detection  lead  naturally  to  procedures 
which  suppress  interference  owing  to  the  nonlinear  nature  ol  satellite 
transponders.  In  general,  optimization  and  evaluation  of  various  modulation 
schemes  used  In  satellite  data  transmission  were  among  the  objectives  of  the 
present  research  effort.  Various  form  of  nonlinear  filtering  and  equalization 
were  also  investigated  [29,33]. 

2.  SUMMARY  OF  RECENT  PROGRESS 

In  the  last  year,  several  significant  contributions  to  the  theory  of  partition 
detection  and  estimation  were  introduced.  Specifically,  a  class  of  8  sequential 
adaptive  detectors  for  underwater  non-gaussian  noise  environments  were 
considered.34  A  new  approach  was  developed  which  permits  the  design  of  non¬ 
linear  robust  detectors  and  estimators  for  efficient  operation  in  dependent 
noise.35 

In  addition,  further  contributions  were  made  to  the  problem  of 
optimization  of  satellite  channels  to  insure  Improved  data  transmission.36,37 

Finally,  detailed  study  of  performance  of  differentially  coherent  QPSK 
systems  in  slowly  fading  Rician  channels  was  completed.33 


3.  STATE  OF  THE  ART  AND  PROGRESS  DETAILS 

A.  Sequential  Adaptive  Detectors  for  Underwater 
Non-Gaussian  Noise  Environments  [34] 

The  problem  of  sequential  signal  detection  in  the  presence  of  slowly 
varying  not  necessarily  Gaussian  noise  was  treated  by  introducing  the  sequential 
generalized  linear  tests  (SGLT)  and  the  sequential  generalized  quantiles  statistic 
Tl*.  Detectors  based  on  these  tests  combine  the  good  qualities  of  order  and 
partitioning  statistics:  robustness,  increased  efficiency,  and  ease  of 
implementation.  In  particular,  the  sequential  Ti*  detector  (STI*D)  was  analyzed 
in  detail.  It  was  shown  that  the  same  detector  structure  is  needed  for  fixed- 
sample  and  sequential  optimum  operation,  that  the  STi*D  is  asymptotically 
optimum  and  more  efficient  than  the  corresponding  sequential  partition 
detector,  and  that  its  transmission  rate  is  up  to  four  times  that  of  its  fixed* 
sample-size  counterpart. 

Furthermore,  some  truncated  and  curved  boundary  decision  rules  in 
sequential  operation  of  the  detector  were  introduced  to  improve  the  efficiency 
even  in  extremely  low  signal- to- noise- ratio  (SNR)  environments,  which  are 
often  encountered  in  underwater  acoustics.  These  decision  rules  also  eliminate 
the  influence  of  occasionally  unbounded  sample  sequences  which  are  an  integral 
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part  of  sequential  detectors  operating  in  severe  noise.  On  the  other  hand, 
asymptotic  measures  of  performance  were  evaluated  under  relatively  high  SNR 
operation  of  the  detector,  encountered  in  some  practical  situations,  such  as 
active  sonar,  and  their  performance  was  compared  to  the  exact  and  computer 
simulated  results. 

B.  Non-Linear  Robust  Detection  and  Estimation  in 
Dependent  Noise  [35] 

Adaptive  optimal  detection  and  estimation  algorithms  were  considered. 
The  algorithms  are  generaily  non-linear.  The  development  is  more  general 
than  most  previous  treatments  of  non-linear  detection  and  estimation  in  that  it 
is  permitted  for  the  noise  to  be  dependent,  subject  only  to  a  mixing  condition. 

Optimal  detection  procedures  depend  on  knowledge  of  the  noise 
probability  density  function  (pdf);  therefore,  a  methodology  for  estimation  of 
pdf’s.  The  approach  permits  estimation  based  on  dependent  samples,  and  was 
extended  to  estimation  of  conditional  pdf’s.  The  latter  are  required  for  optimal 
detection  in  dependent  noise. 

The  methodology  adopted  to  the  problem  of  estimation  of  pdf's—piecewise 
polynomial  approx  lm  ation—  re  quires  knowledge  of  quantiles  and  truncated 
moments  of  the  density  to  be  estimated.  This  led  to  the  consideration  of 
estimation  of  these  quantities;  the  Robbins-Monro  Stochastic  Approximation 
(RMSA)  algorithm  was  employed  for  this.  A  proof  of  convergence  of  the 
RMSA  algorithm  with  a  dependent  input  sequence  is  given  under  the 
assumption  that  the  sequence  is  mixing.  This  is  an  Important  generalization  of 
the  available  theory. 

Adaptive  non-linear  filtering  and  prediction  algorithms  were  also  presented. 
These  are  based  on  piecewise  approximation  to  a  regression  function  involving 
only  estimated  partition  moments  and  quantiles.  All  of  the  algorithms  are 
simple  enough  to  be  adaptively  implemented  in  real  time,  and  several 
illustrations  of  performance  were  presented. 

C.  Nonlinear  Equalization  in  Satellite  Data 

Transmission  Systems  [38,37] 

The  problem  of  designing  and  evaluating  the  performance  of  a  minimum 
mean-square  error  equalizer  (MMSEE)  for  binary  PSK  transmission  over 
band-limited  nonlinear  satellite  channels  was  considered.  The  effect  of  inter¬ 
symbol  interference  followed  by  AM/AM  and  AM/PM  conversions  were  taken 
into  account  while  optimizing  the  performance  in  the  presence  of  the  downlink 
white  Gaussian  noise.  In  analyzing  the  problem,  the  decision  was  made  on  a 
typical  signal  in  a  received  sequence  taking  into  account  past  and  future 
interfering  signals,  i.e.,  ISI. 

As  an  illustrative  example  of  the  ?  typical  channel  model  was 

considered  in  detail.  Based  on  the  analysis,  an  alternative  receiver  structure 
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which  is  more  suitable  for  implementation  was  introduced.  The  taps  gain 
coefficients  for  minimum  mean-square  error  between  the  received  sample  and 
the  actual  transmitted  bit,  were  obtained  using  numerical  methods.  The 
performance  of  the  equalizer  using  computer  simulation  techniques  and  it  was 
shown  that  significant  performance  improvement  over  the  single-sample  sign 
detector  can  be  obtained. 

The  results  of  reference  36  were  then  generalized  to  a  satellite 
transmission  systems  using  quadrature  phase  shift  keyed  (QPSK)  signals.37 

D.  Performance  of  Differentially  Coherent  QPSK  Systems 
in  Slowly  Fading  Rician  Channels  [38] 

A  detailed  analytical  approach  to  the  evaluation  of  the  error  rate 
performance  of  a  differentially  coherent  QPSK  system  with  linear  diversity  in  a 
frequency  selective  slow  fading  channels  with  Rician  distributed  signal 
envelopes  was  presented. 

First,  a  real-time  channel  parameter  estimation  technique  was  introduced, 
followed  by  the  error  probability  evaluation  using  linear  diversity  in  a  Rician 
fading  environment.  Gaussian  and  exponential  delay  autocorrelation  functions 
were  assumed,  and  the  Intersymbol  interference  for  one-  and  two-symbols 
adjacent  to  decision  intervals,  various  data  rates  and  phases,  SNR,  as  well  as 
the  ratio  of  the  direct  signal  component  power  to  the  multipath  power  were 
taken  into  account  and  evaluated. 

Extensive  simulation  results  suggest  that  if  the  normalized  data  rate,  d,  is 
less  than  .2,  the  system  is  insignificantly  affected  by  intersymbol  interference. 
If  the  diversity  order  L  =  2,  d  >  0.4  or  if  L  =  4,  d  >  0.6,  the  second  symbols 
on  either  side  of  the  detection  interval  start  generating  about  twenty  per  cent 
intersymbol  interference  as  compared  to  the  first  pair.  It  is  also  interesting  to 
note  that  below  these  rates  the  high  SNR  can  counteract  the  loss  in  a  fading 
environment.  But  for  a  high  data  rate,  the  fading  generates  the  dominant 
metric  probability  of  error  for  various  detected  phases.  This  is  a  result  not 
observed  in  linear  channels.  Typical  performance  results  are  shown  in  Figures 
1  and  2. 
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B.  A  RAISED  POWER  TEST  AND  SOME  OF  ITS  APPLICATIONS 
IN  ARRAY  SIGNAL  PROCESSING 

Professor  S.  Unnikrishna  Pillai 

Unit  IE7-2 

1.  OBJECTIVES 

The  problem  of  estimating  the  number  of  systematic  (signal)  components  in  data  ob¬ 
tained  by  a  collection  of  spatially  deployed  sensors,  simultaneously  exposed  to  spatially 
white  noise  is  considered  here.  A  method  heretofore  proposed,  based  on  the  ratio  of  arith¬ 
metic  mean  to  geometric  mean  of  the  sample  covariance  matrix  eigenvalues,  is  reviewed 
and  an  alternative  based  on  the  higher  powers  of  the  eigenvalues  is  developed.  The  lat¬ 
ter  test  is  found  to  be  applicable  in  cases  where  the  conditions  for  the  former  are  not 
met.  Specifically  this  is  used  to  analyze  the  performance  of  a  modified  spectrum  estimator 
utilizing  a  covariance  matrix  raised  to  a  higher  power. 

2.  SUMMARY  OF  RECENT  PROGRESS 

One  of  the  important  problems  in  sonar,  radar  and  geological  acoustic  data  processing 
is  that  of  resolving  an  unknown  number  of  closely  spaced  sources  from  information  collected 
at  a  set  of  sensor  elements.  Several  approaches  have  been  reported  for  solving  this  and 
related  problems  [1-6],  In  particular  determination  of  the  unknown  number  of  signals, 
when  the  additive  noise  is  element  to  element  uncorrelated  and  of  equal  variance  has  been 
formulated  as  a  sequential  hypothesis  testing  problem  [6-8].  Under  suitable  preprocessing 
such  as  spatial  smoothing,  this  test  has  also  been  extended  to  coherent  cases  [9,10],  A 
crucial  assumption  in  all  these  cases  is  that  the  number  of  unknown  sources  present  in 
the  scene  is  strictly  less  than  the  total  number  of  sensor  elements  used  to  collect  the  data. 
When  this  is  not  the  case  this  test  would  always  indicate  the  number  of  sources  present 
to  be  the  same  as  the  number  of  sensor  elements,  thus  making  the  results  meaningless.  In 
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what  follows  we  review  this  situation  and  derive  an  alternate  test  to  estimate  the  number 
of  sources  in  that  situation.  This  is  made  possible  by  deriving  the  asymptotic  distribution 
of  this  new  test  criterion. 

The  usefulness  of  this  test  to  estimate  the  number  of  unknown  arrivals,  is  illustrated 
in  section  4.  There  it  is  shown  that  if  the  source  scene  is  apriori  known  to  consist  of 
point  sources  in  uncorrelated  noise  of  equal  variance,  then  by  a  ‘raised  power  technique’ 
on  the  array  output  covariance  matrix,  one  can  achieve  the  same  type  of  performance  as 
in  the  case  of  the  so  called  Multiple  Signal  Classification  (MUSIC)  technique  developed 
by  Schmidt  [4,8].  The  advantage  here  is  that  this  ‘raised  power  technique’  avoids  the 
determination  of  the  eigenvectors  of  the  array  output  covariance  matrix  and  their  further 
partitioning  into  signal  and  noise  subspaces,  thereby  eliminating  the  associated  decision 
errors.  The  raised  power  test  developed  in  section  3.1  is  used  to  determine  the  penalty  for 
raising  the  covariance  matrix  to  a  higher  power  in  this  new  technique.  It  is  shown  that  the 
number  of  samples  used  in  estimating  the  covariances  have  to  be  increased,  as  the  power  to 
which  the  matrix  is  raised  is  increased.  Simulation  results  are  presented  to  support  these 
conclusions. 

3.  STATE  OF  THE  ART  AND  PROGRESS  DETAILS 

We  consider  the  problem  of  estimating  the  number  of  “systematic"  (signal)  compo¬ 
nents  present  in  the  data  which  is  measured  as  an  M-component  random  vector  x(tn);  n  = 
1,2,...  N.  In  general  the  observed  data  2(f„ )  have  two  unobservable  random  parts. 

$(**)  =  St(**)  +  Sl(tn);  n  =  l,2,...iV  (1) 

Here  a(t„)  and  Ji(*»)  are  random  vectors  with  zero  mean  and  covariances  given  by 
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E[vm  (t)vT  (t)]  =  a3 1  (3) 

Usually  t/(f)  can  be  interpreted  as  composed  of  errors  of  measurement:  uncorrelated  com¬ 
ponents  with,  equal  variances.  In  those  cases  where  data  is  received  at  a  set  of  M  sensors, 
this  would  correspond  to  spatially  white  noise  inducing  an  uncorrelated  random  compo¬ 
nent  in  sensor  of  variance  a3 .  Similarly  u(t)  can  be  interpreted  as  made  up  of  useful 
signal  parts  and  is  supposed  to  be  in  a  space  of  dimension  K  less  than  M.  As  a  result  £  is 
of  rank  K,  and  the  observable  covariance  matrix  R  —  £  +  a3 1,  has  an  eigenvalue  a2  with 
multiplicity  (M-K). 

Many  practical  situations  can  be  attributed  to  this  model.  One  of  its  uses  is  in 
exploratory  studies,  where  the  observer  is  interested  in  reducing  the  vast  amount  of  avail- 

i 

able  measurements.  Another  interesting  area  of  application  is  in  modem  radar  and  sonar 
problems  where  a  set  of  sensors  is  used  to  collect  data  from  spatially  distributed  targets. 
There,  the  above  model  would  correspond  to  the  situation  where  M  sensor  elements  are 
used  to  receive  signals  from  K  not  totally  coherent  point  sources  in  the  presence  of  spatially 
uncorrelated  noise  of  equal  variance  a3 . 

In  applications  where  one  makes  use  of  eigenvectors  associated  with  the  systematic 
part  (unequal  eigenvalue  part)  and  the  noise  part  (those  corresponding  to  the  equal  eigen¬ 
value  a3)  separately,  it  is  necessary  to  partition  the  eigenvectors,  and  hence  the  eigenvalues, 
into  two  groups:  those  belonging  to  the  systematic  or  unequal  part  and  those  belonging  to 
the  equal  one.  This  partitioning  is  trivial,  given  the  ensemble-averaged  covariance  matrix 
R.  However,  when  only  an  estimate  R  of  R  such  as 


R  =  **(*-)af  (*•) 


(4) 


n*  1 


is  available,  this  is  no  longer  a  trivial  problem  since  with  probability  1,  none  of  its  eigen¬ 
values  will  be  equal  [ll].  Statistical  tests,  based  on  the  asymptotic  multivariate  Gaussian 
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nature  of  the  eigenvalues,  have  been  devised  for  the  real  case  by  Anderson  [12],  Law  ley  [13], 
and  others  [14],  to  estimate  the  number  of  equal  eigenvalues  with  a  certain  confidence  level. 
Similar  results  for  the  complex  case  are  reported  in  [15,  16].  Specifically  with  u(t)  and 
v(t)  as  complex  Gaussian  random  processes  and  R  as  in  (4),  its  eigenvalue-decomposition 
would  give 

R  =  UllP  (5) 

where  (f  stands  for  complex  conjugate  transpose) 

U  —  [ttt ,  u3 , . . .  Ujy  ]  (6) 


and 


’ii  0  ...  0  * 

0  £q  ...  0 

•  •  •  • 

•  •  •  • 

•  •  •  • 

.0  0  ...  £u . 


(7) 


Here  41,a3,...itu,  are  the  M  orthonormal  eigenvectors  and  ,  are  the  cor¬ 

responding  real  eigenvalues  of  R.  A  test  criterion  for  the  equality  of  the  (M-K)  lowest 
eigenvalues  in  the  complex  case  has  been  shown  to  be  [16] 


1 


which  is  the  ratio  of  the  arithmetic  mean  to  the  geometric  mean  of  the  corresponding 
sample  eigenvalues.  Using  an  asymptotic  analysis,  it  can  be  shown  that 

N(M  -  K)lnSK  (£)  -  *3(<i)  (9) 

where 

d  =  (M  -  K)2  -  1  (10) 

represents  the  degrees  of  freedom  associated  with  the  x2  random  variable  in  (9).  This 
will  allow  one  to  examine  the  equality  of  the  lowest  (M-K)  eigenvalues  by  computing 
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2 N(M —K)tnSK  (£)  and  accepting  the  hypothesis  if  it  is  less  than  the  upper-tail  significance 
point  of  the  appropriate  x3  distribution  for  some  suitable  significance  level  a.  The  test 
will  proceed  with  K=l,2,. . .  and  the  first  value  of  K  for  which  the  test  succeeds  is  taken 
as  the  estimate  for  K[7,8,9j.  Clearly  the  test  is  meaningful  only  if  K  <  M ,  i,e;  so  long  as 
the  number  of  targets  is  less  than  the  number  of  sensor  elements. 

In  what  follows  we  derive  an  alternate  test  based  on  higher  powers  of  the  eigenvalues 
and  derive  its  asymptotic  distribution.  The  usefulness  of  this  test  is  briefly  discussed  in 
section  4  for  a  specific  application  where  the  conventional  test  based  on  (9)  turns  out  to 
be  inappropriate. 

A.  A  Raised  Power  Test  and  Its  Statistics 

With  symbols  as  defined  above  and  A^Aj,...  Aw  representing  the  roots  of  R,  the 
hypothesis  Hx  under  test  is  XK+1  —  +2  =  ...  =  AM.  The  required  criterion  for  the 

equality  of  the  (M-K)  lowest  eigenvalues  of  based  on  their  p"*  powers  would  be 


t  <r)/(n  <) 

L  '  imK  +  l  '  ’  \»JC+1  7 


1  1 
P 


M-K 
1 


-(ra  t  <)7(  n  *)M-* 

'  i~JC  +  l  7  7  \mK  +  l  ' 


We  will  show  that 


21V  (M  -  K ) 


f»AjcU,p)  ~  x3(<0 


where  d  is  as  defined  in  (10). 
Proof 

Following  Anderson  [17],  define 


in) 


(12) 


dl=VN(ti-XKtl)  ;  i  =  K  +  l,...M 


(13) 
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With  (11)  and  (14)  in  (12)  we  have 


2JV(M  -  K) 


tnAK{£,p)  =  — 


N 
P  L 


(M-K) 


1,  ('^n 


=  — 
p 


t  (.. 

v  i-X  +  i  x 


pd< 


p(p  ~  l)g  , 

l  «  e  i  •  •  • 


_  V  _ «_  +  \ 

-2[ii£=*lfc{1+( - i£*. 


+  _p(pi^1E£_ 

1  a  «  r  /  <  r  i  l  ■ 


,  tf1/3  (M  -  JO  A,  +  l  2JV(Af  -  K)  \*K  + x 


E4  .  E<?  ‘ 

+  rrrr; - H  ... 


N 

P 


*l/aA*  +  i  2WAj> +l 
{M-K)f  p£4 


_  _ i  p(p  - 1)  E 

P  \N'/*{M-K)\k  +  1  2N[M‘-K)\*k  +  1 

^(EO3  ,  ^  E*  .  E* 


(E4)3  \ 


2JV(Jl/-JT)ni  +  1  •  -J  JV/’A,*,  ■  2iVA>+, 


+  .  .  . 


(15) 


It  can  be  shown  that  the  limiting  distribution  of  dK  +  l,dK  +  3,...dM ,  is  the  same  as  the 
distribution  of  the  eigenvalues  of  a  hennitian  matrix  C/33  =  (tt0 ),  *,  j  =  A  + 1,  A +2, . . .  M. 
Further,  for  circular  Gaussian  data  [18]  elements  of  C/33  are  asymptotically  independent, 
identically  distributed,  zero  mean  random  variables  with  common  variance  Xi.  , .  Then  it 
follows  that  (15)  has  the  limiting  distribution  of 


SECTION  HI:  INFORMATION  ELECTRONICS 


(£-)*] 

JT  +  l  <■'<*/  i-JT  +  1  '  '  '  JC  + 1  7  J 

=  477 ^ 'u*y'a  +  ^77 [  27  ““ “  (m- a-) (  ^  u<<)  (16) 

AT*  ,'|4y  +  l  '  '  \«JC  +  1  '  J 

Thus  (53  Kf)/4  +  1  is  asymptotically  x3  with  (M-K)(M-K-l)  degrees  of  freedom;  inde- 
*#y 

pendently,  [53  u2.  -  (53  ti<i  )3  /  (A/  -  if )  ]/ +  x  is  asymptotically  x3  with  (M-K-l)  degrees  of 

«  i 

freedom.  Thus  (15)  has  a  limiting  x3  distribution  with  [(M  —  if)2  —  1]  degrees  of  freedom 


and  hence 


2iy(Af  -  K) 


tn&K  (Lp)  ~  X2  (d) 


where  d  is  as  in  (10).  This  completes  the  proof.  On  comparing  (8)  and  (17)  it  now  follows 
that  the  effect  of  raising  the  eigenvalues  to  a  higher  power  p  is  to  reduce  the  effective 
number  of  samples  by  a  factor  of  p. 

4.  AN  APPLICATION  FROM  ARRAY  SIGNAL  PROCESSING 

We  illustrate  the  usefulness  of  the  above  test  with  an  example  from  sensor  array  signal 
processing  [19].  The  objective  in  array  signal  processing  is  to  estimate  the  number  of  tar¬ 
gets  and  their  associated  parameters  such  as  angles  of  arrival,  power  levels  etc.  by  making 
use  of  the  data  available  at  the  sensor  element  outputs.  The  Maximum  Likelihood  Method 
(MLM)  of  spatial  spectrum  estimation  is  one  of  the  earliest  superresolution  techniques 
proposed  [20],  and  is  still  today  an  important  tool  for  estimating  directions  of  multiple  sig¬ 
nal  arrivals.  The  MLM  requires  the  inverse  of  the  matrix  of  covariances  obtained  between 
all  array  element  pairs.  This  is  then  scanned  by  a  pointing  vector  to  give  an  estimate 
of  power  vs.  pointing  angle.  The  maxima  of  this  function  are  then  used  as  estimates  of 
arrival  angles.  More  recently,  elegant  techniques  based  on  a  vector  space  decomposition  of 
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the  covariance  matrix  have  been  developed  [4,8]  which  give  results  superior  to  those  of  the 
MLM.  In  particular,  the  resolution  capability  of  the  MLM  deteriorates  with  decreasing 
signal  to  noise  ratios  (SNR).  Furthermore,  estimates  of  arrival  angle  for  sources  closely 
spaced  in  angle  tend  to  be  more  biased  under  such  conditions,  even  when  obtained  with 
a  perfectly  estimated  covariance  matrix.  On  the  other  hand,  with  a  perfectly  estimated 
covariance  matrix,  the  vector  space  methods  always  give  unbiased  estimates  even  when  the 
source  directions  are  only  infinitesimally  apart  and  have  superior  resolvability  compared 
to  the  MLM-estimator. 

To  understand  the  superior  performance  of  the  vector  space  spectral  estimators,  con¬ 
sider  the  array  processing  problem  described  above  where  an  M  element  array  located  at 
&i  t  >  •  •  •  &u »  w'th  respect  to  some  frame  of  reference,  receives  signals  from  K  narrowband 
sources  which  arrive  as  plane  waves  with  wave  vectors  , . . .  ^  in  that  frame.  Let 

the  k*K  arrival  at  the  origin  be  denoted  (t);  Jb  =  1,2,...  K.  Then  the  element  output 
Xi  (t)  at  time  t  can  be  written  as  [8] 

*<(*)  =  £  MO  «-'(***/*)  +  MO  (18) 

km  l 

where  e  represents  the  speed  of  propagation  in  the  medium  and  n<  ( t )  stands  for  additive 
noise  which  we  will  take  to  be  uncorrelated  from  element  to  element  and  of  equal  variance, 
No  .  In  vector  form  (18)  becomes 


*(0  *  AT  i(t)  +  a(t) 


(19) 


i(0  =  (M0»  M0»---  *m(01T 

*(0  =  MO.  MO.---  M01T 


where 


(20.a) 

(20.b) 


SECTION  HI:  INFORMATION  ELECTRONICS 
Here  T  stands  for  the  transpose  and  *  represents  complex  conjugation.  Further, 

A  =  [£(*),  six,),...  £(**)]  (21) 

is  an  MxK  matrix  consisting  of  signal  direction  vectors  given  by 

£($/*)  =  e,(i-  K*/a)]r  (22) 

The  correlation  matrix  of  the  received  signal  vector  x(t)  is  then  an  MxM  matrix, 

R  =  E[z‘  (i)£r  (*)]  =  ARs  s  Af  +  No  I  (23) 

where  R33  =  i?(i*(f)er  (f)]  represents  the  KxK  source  correlation  matrix  and  t  denotes 
complex  conjugate  transpose. 

When  the  sources  are  not  completely  correlated,  ARS  3  A1  will  be  of  rank  K  and 
hence  the  lowest  eigenvalue  of  R,  which  will  equal  JV0  will  repeat  (M-K)  times.  Letting 
ii,  i  —  K  +  1,  K  +  2,  . ..M,  represent  the  normalized  eigenvectors,  corresponding  to  the 
(M-K)  repeated  eigenvalue,  we  get 

R§4  »  (ARSsAf  +  N0I)Si  —  NoZi  (24.a) 

which  implies 

A'&i  =  0  (24.b) 

or  equivalently 

£?(&)*  =0;  foraU  i  =  K  +  l,if  +  2,...A/  and  Jfe  =  1,2,... K  (25) 

Thus  the  K  signal  direction  vectors  are  orthogonal  to  the  subspace  spanned  by  the 
noise  eigenvectors.  This  is  the  basis  for  the  MUSIC  estimator  Pu  v  (u)  proposed  by  Schmidt 
[21j  where 

_ 1 _ 

Itex  +  U  Sjc+a**”  4*1*  £(ii)la 


Puvik.)  = 


(26) 
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This  however  can  be  rewritten  as 


Pun  (id) 


=  (tf’SO'))' 


(27) 


where  U  is  the  modal  matrix  of  R.  Using  (24)  and  (25)  we  see  that  W  S_{y)  results  in  a  set 
of  new  direction  vectors,  which  for  true  directions  of  arrival  y^,;k  =  1,2...  K,  would  span 
the  wAtTMi  subspace  as  the  one  spanned  by  the  eigenvectors  corresponding  to  the  K  largest 
eigenvalues  (signal  subspace  eigenvectors)  and  from  this  it  follows  that  1/Pmu  (id)  given  by 
(27)  equals  zero  iff  u  €  {idi  jidj  •  •  •  vK  }. 

Within  this  framework  the  Maximum  Likelihood  spectrum  estimator  of  Capon  can  be 
written  as  [20] 


P(u)  - - = - 

£f(id)P-l£(id) 

Using  similar  manipulations  as  in  (27)  we  have 


(28) 


1 

P(*) 


(fr^Od))* 

No 


Pi 


Pi 


Pk 


VIM 


(29) 


where 


P»  - 


1  +  A*/JV0 


>0;  k—  1,2,...K 


(30) 


As  before  U  is  the  modal  matrix  of  R  (or  ARssA*)  and  A,  is  the  Ith  non  zero- 
eigenvalue  of  the  noiseless  part  {ARS  s  Af)  of  R.  Note  that  positive  definiteness  of  the 
signal  correlation  matrix  implies  that  A<  >  0  and  hence  Hi  <  1  for  all  i  =  1, 2, . . .  K. 
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The  improvement  in  performance  in  terms  of  resolution  for  the  MLM  spectrum  es¬ 
timator  will  increase  in  SNR  is  well  known  [21,22].  To  see  this,  note  that  the  reason  for 
the  inferior  resolvability  for  MLM(i/),  as  compared  to  Puv  (u),  is  because  of  the  full  rank 
diagonal  matrix  in  (29).  The  quantity  of  interest  there  (as  a  function  of  SNR)  is  nh  as 
given  by  (30).  As  increases  in  the  SNR  would  result  in  increased  values  of  the  Xks,  the 
Mb's  would,  in  turn,  decrease.  For  large  SNR’s,  fik  — ►  0,  k  —  1,2, ...  K,  and  from  (27)  and 
(29)  it  follows  that  MLM  spectrum  approaches  the  one  obtained  by  MUSIC  method. 

This,  however,  raises  the  following  question: 

In  the  presence  of  a  multiple-source  scenario,  can  the  MLM  spectral  estimator  be 
modified,  without  sacrificing  its  simplicity,  such  that  the  modified  spectral  estimator  would 
tend  to  attain  the  superior  properties  of  the  vector  space  spectral  estimators  such  as 
MUSIC?  A  comparison  between  (27)  and  (29)  shows  that,  they  differ  (barring  N0)  only 
because  of  the  nonzero  values  of  /x*,  k  =  1,2,. . .if,  appearing  in  (29).  Thus  for  true 
direction  vectors,  whereas  MUSIC  is  able  to  satisfy  the  requirement  of  zero  norm  in  the 
quadratic  product  (27)  exactly,  the  MLM  is  only  able  to  do  this  approximately.  Hence  the 
degradation  in  performance  for  MLM  due  to  the  nonzero  values  of  /xfc  can  be  eliminated 
by  suitably  modifying  it  so  as  to  drive  the  equivalent  p'k3  to  zero.  This  can,  in  fact  be 
achieved  by  modifying  the  MLM  estimator  as 


=  (31) 
where  p  is  some  positive  number  greater  than  one.  As  in  (29)  it  can  be  shown  that 

A  Q 


PM 


Ac 


(32) 


nBnafHBnonoi  ummiurujiv  ifH  ranannv  gOTpqciwwwmAwrjtnK7ucy>  njr  njtnin  v*  v*  re  ^  lh.  l-.  wv 
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Since  0  <  nk  <  1  for  all  k  =  1, 2, ...  K,  — ►  0  as  p  — ♦  oo  for  ail  k.  Thus  asymptoti¬ 
cally,  the  modified  spectral  estimator  given  by  (31)  would  approach  (26)  and  hence  would 
have  all  the  desirable  properties  of  the  vector  space  spectral  estimators.  We  note  that  this 
result  is  independent  of  SNR  and  hence,  in  particular,  is  true  at  low  SNRs  also. 

We  remark  that  the  modification  suggested  here  is  different  from  that  introduced  by 
Pisarenko  in  connection  with  estimating  spectra  by  means  of  nonlinear  functions  of  the 
covariance  matrix  [23],  He  proves  that  for  a  spatially  stationary  field,  P*r  (v),  given  by 

(33) 

convciges  in  probability  to  the  true  spectra  under  some  mild  conditions  (See  corollary  after 
remark  4  in  [23]).  This  also  follows  from  a  general  result  known  in  spectral  theory  (See 
Th.3.3.3  in  Kadison  and  Ringroee  [24]),  which  states  that  for  a  bounded  operator  A 

As  N  — *  oo,  llA^]]1^  has  a  limit.  (34.a) 

This  limit  is  known  as  the  spectral  radius.  Note  that  (33)  follows  from  (34.a)  by 
defining  the  norm  of  A  to  be 

IW|-fi»(fcM£(|f)  (34.b) 

and  substituting  R  for  A  in  (34.b). 

However  the  modification  suggested  in  (31)  discretizes  the  spectrum  and  is  valid  even 
for  spatially  nonstationaxy  situations  such  as  correlated  arrivals.  To  see  this  note  that 
as  p  — ♦  oo,  equation  (32)  represents  the  length  of  the  direction  vector  £(z/)  projected 
onto  the  subspace  spanned  by  the  minimal  eigenvector (s)  (corresponding  to  the  lowest 
eigenvalue).  Using  (24)  and  (25),  this  length  (after  projection)  will  be  zero  only  for  those 
distinct  direction  vectors  corresponding  to  actual  targets  whose  total  number  is  less  than 
the  number  of  array  elements.  As  a  result  the  limiting  spectrum  in  (31)  will  consist  of  a 
set  of  impulses  representing  those  distinct  direction  vectors.  Further  from  the  arguments 
following  (32),  it  is  clear  that  the  spectrum  so  obtained  corresponds  to  the  true  spectrum 
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only  when  the  data  itself  represents  point  sources  and  thus  the  modified  MLM  estimator 
in  (31)  is  asymptotically  optimal  for  a  point  source  scenario. 

When  the- ensemble  averaged  values  of  covariances  are  known,  the  value  of  p  to  which 
the  covariance  matrix  R  is  raised  can  be  arbitrarily  large.  In  that  case  p  has  the  same 
effect  on  all  (M-K)  diagonal  entries  in  (32)  corresponding  to  the  noise  subspace  part  and 
hence  can  take  any  value.  However,  when  only  an  estimate  R  of  the  covariance  matrix  is 
known,  as  mentioned  before,  the  probability  of  any  two  of  its  eigenvalues  being  equal  is 
zero  111],  and  as  a  result  p  win  have  a  different  effect  on  those  entries  corresponding  to  the 
noise  subspace  part. 

The  question  about  the  choice  for  the  value  of  p  to  which  the  matrix  R  may  be 
raised  can  however  be  readily  answered  by  making  use  of  the  asymptotic  result  developed 
in  section  3.  It  is  shown  there  that,  when  R  is  estimated  using  N  samples,  the  result 
of  raising  its  eigenvalues  to  a  higher  power  p  in  an  algorithm  is  to  reduce  the  effective 
number  of  samples  in  averaging  to  N/p.  Since  R  is  hermitian,  the  effect  of  raising  it  to 
a  higher  power  p  is  identical  to  replacing  its  eigenvalues  li,l9t.  ..lu ,  by  /^,Z!|,. .  .Pu ,  and 
hence  from  the  above  discussion  it  now  follows  that  when  R  is  made  use  of  in  a  spectrum 
estimation  algorithm,  as  in  (31),  it  is  necessary  to  increase  the  number  of  samples  used  in 
averaging  by  a  factor  of  p  . 

This  observation  is  also  supported  by  simulation  results  shown  in  Fig.  1  and  Fig.  2. 
There  a  four  element  aperiodic  array  receives  signals  from  two  sources  located  along  85° 
and  94*  with  respect  to  the  line  of  the  array.  The  SNR  in  each  case  is  taken  to  be  7.5  dB 
and  6  dB  respectively.  Fig.  l.a  shows  the  spectrum  obtained  by  the  MUSIC  method  where 
25  samples  are  used  in  averaging  R.  The  corresponding  MLM  spectrum  is  shown  in  Fig. 
l.b.  The  same  covariance  matrix  is  raised  to  the  fourth  power  to  compute  the  modified 
MLM  spectrum  given  by  (31)  and  this  is  shown  in  Fig.  2.b.  As  before  the  targets  are 
barely  resolvable  here.  As  required  by  the  asymptotic  result  in  section  3,  the  number  of 
samples  used  in  averaging  R  is  increased  to  100  for  performance  comparable  to  the  MUSIC 


element  aperiodic  array  with  element  locations  0,2.1, 4.3  and  5.8  (half  wavelenghts) 
receives  signals  from  two  sources  located  along  85*  and  94*  with  SNRs  7.5  dB,  6dB 
respectively,  (a)  The  MUSIC  estimator  using  25  samples  in  averaging  R.  (b)  The 
Capon  estimator  with  25  samples,  (c)  The  Capon  estimator  with  100  samples. 
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Fig.2  Effect  of  raising  the  sample  covariance  matrix  to  a  higher  power  on  the  spatial  spec¬ 
trum.  Details  are  as  in  Fig.l.  (a)  The  MUSIC  estimator  with  25  samples,  (b)  The 
Modified  Capon  estimator  with  25  samples  and  p  a  4  (See  (31)).  (c)  The  Modified 
Capon  estimator  with  100  samples  and  p  =  4. 
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estimator  and  the  spectrum  recomputed  using  (31)  is  shown  in  Fig.  2.c.  The  improvement 
in  performance,  obtained  with  increasing  the  number  of  samples  used  in  averaging  and 
then  raising  the  covariance  matrix  to  a  higher  power  prior  to  computing  the  spectrum  can 
be  seen  by  comparing  Fig.  l.c  and  Fig.  2.c. 

We  have  seen  that  by  raising  the  covariance  matrix  to  a  higher  power  it  is  possible  to 
realize  spectrum  estimators  with  performance  levels  comparable  to  those  that  make  use  of 
eigenvectors  of  the  covariance  matrix  in  their  algorithm.  Although  there  is  an  advantage 
in  not  performing  such  an  eigenvalue-eigenvector  analysis  here,  the  disadvantage  lies  in  the 
larger  number  of  actual  samples  required  in  estimating  the  covariances. 

In  the  context  of  estimating  the  number  of  signal  (systematic)  components  from  noisy 
observations  at  an  array  of  sensor  elements,  a  new  test,  that  is  based  on  higher  powers  of 
the  eigenvalues  of  the  array  output  sample  covariance  matrix  is  developed.  This  is  then 
used  in  n  dyzing  the  degradation  in  performance  of  a  modified  spectrum  estimator  which 
makes  use  of  higher  powers  of  the  covariance  matrix  in  estimating  the  arrival  angles  of  the 
point  sources  in  the  scene.  For  comparable  performance  as  in  the  conventional  method, 
the  modified  method  is  shown  to  require  larger  number  of  data  samples  in  estimating  the 
interelement  correlations. 
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